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SXUOPSIS 


A IdaetiG theory approach is used for investigating the structure 
of shoch ^ve in a dissociating diatomic gas. An axially elastic 
spheroxcylinder molecular model mhich permits transfer of rofetional 
and vibrational energy at a collision is proposed for a vlbrationally 
relaxing diatomic gas. (Che axial elasticity has been Introduced to 
account for the vibrational degrees of freedonu Ihe distribution 
function in such a case depends on a set of 14 variables t ^ 

v^ere B and <}> are reqxiired to describe the vibrational state of a 
molecule. 

(Hie maldiQiia'H.cal description used for vroiicing out the collision 
dynamics of the axially elastic sphere cylinder molecules, closely 
follows the one used for the rigid spherocylinders by Haight and Lundgre 
(The geometry ai^ d^mamics of a collision betweai two molecules are studi 
at first, 

(Kie frequency of binary collisions in the diatomic gas is obtained. 
Ihis is used in writing down the collision integral in the Boitjjmann 
equation. Bor the Boltzmann equation, thus formiilated, an H'-Iheorem If 
given and throxigh this theorem, the equilibrium solution of the Boltsaaan' 
equation Is obtained, 

Bor the shock structure- problem, a trimodal and a qUfsirimodal 
ansatz are proposed euad a moment method is used to Tind the wn 
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vibrational energy are used in addition to the conservation eq,nations. 
Itor liie quadramodal approach, the rotational energy, vibrational energy 

p 

and the C moment of the Boltzmann equation are used to determine the 

weight functions agai% in addition to the conservation equations, lUie 

collision integrals in these laDment equations are 22 fold out of which 

16 fold integrations are carried out analytically arai the remaining 

integrations were performed numericcjlly using the Monte Carlo technique, 
©le set of equations in the unknown weight functions have been solved 
nunerioally to get the shock profiles for a range of Isfech numbers and 

upstream pressures for nitrogen ^s* 

lurther, for the case of a dissociating gas the Boltzmann equation 
is modified by including the collision terns representing dissociative 
and rscombinative collisions * A separate Boltzmann equation is writter 
for the atom species present in the dissociating diatomic gas. dheae 
two Boltzmann equations are solved using a moment method for the shock 
boundary conditions* The trimodal ansata suggested earlier is used fo' 
the diatomic molecular distribution function and a bimodal distribution 
function for atom specres. 

In addition to the two species conservation equations, the 
rotational and vibrational energy equations have been used for deter- 
mining the wei^t functions. Again the collision integrals in these 
moment equations are also 22 fold and have been evaltiated as in the 
earlier case. The shock profiles have been obtained by solving the se' 
of equiations in the wei^t functions for various Mach number and upstre 
pressures for nitrogen gas* 



SYiropsis 


A kinetic theory approach is used for investigating the structure 
of shook wave in a dissociating diatomic gas. An axially elastic 
sphero-cylinder molecular model vdxLch permits transfer of rotational 
and vibrational energy at a collision is proposed for a vibrationally 
relaxing diatomic gas, Qlie axial elasticity has been Introduced to 
account for the vibrational degrees of freedom* The distribution 
function in such a case depends on a set of 14 variables (s » £ » 2 >23 
where B and ({' are requir>ed to describe the vibrational state of a 
molecule. 

The mathematical description used for working out the collision 
dynamics of the axially elastic spherocyllnder molecules, closely 
follows the one used for the rigid sphero cylinders by Haight and Lundgren. 
The geometry and dynamics of a collision betweai two molecules are studied 
at first. 

The frequency of binary collisions in the diatomlo gas is obtained. 
Kiis is used in writing down the collision integral in the Boltemann 
equation. Jbr the Bolt2imann equation, thus formulated, an H-^Theorem is 
given and through this theorem, the equilibrium solution of the Boltjjmann 
equation is obtained. 

Tbr the shock structure- problem, a trimodal and a quadrimodal 
ansatz are proposed and a moment method is used to find the xuafcnown 
weight functions. To determine Ihe three weight functions in the 
trimodal formulation, Iwo moment equations rotational energy and 
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IKPROOTCPIOf. 


1 ,1 , Phis work is concerned with the problem of shock structure 
in a dissociating diatomic gas. In gas dynamics, it is customary 
to assume the shock to be a surface of discontinuity across which 
two uniform states eocistn. {these two uniform states are related 
by the Eankine - Hrgoniot relations. Ebwever, to describe the 
narrow transition aone in the shock vra.ve it is necessary to consider 
the effects of wLscosity and thermal conductivity. Phis problem 
was solved within the frame work of the Havier-Stokes equations 
for a monatomic gas first by {Taylor (1910) and later by Becker 
( 1921 ), {Thomas (1944) and others. Ehe thickness of shock vave 
v/as found to be of the order of a few mean free paths of the 
upstream. But the whlidity of the IfFavier “Stokes equations for 
the shock phenomenon is questionable as the characteristic length 
for the shock is so small as to break the continuum hypothesis 
made in the formulation of these equations. {Ore results obtained 
by continuum and kinetic theory approaches for the shock problem, 
as given by Gilbarg and laolucci (l953)> show roughly equivaloit 
shock thicknesses. However, the merit of a theory should be 
judged by comparison of its results with those given by experiments. 
Observations made for monatomic gases by {Talbot (1962 ) and 
Sherman and Talbot (1960), after a careful study of the available 
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Kon>-diKiensional density 


Subscripts 

r j Rotational (for T, f and t ) 

t ; (Eraaslational (for T, f and t ) 

V : Vibrational (for T, T and j ) 

R : Rotational (for & ) 

T ! iEranslational (for A ) 

V t Vibrational (for A ) 

eU j Monte Carlo estimate for the integral vdth U trials 



data, bo til experimfintal and theoretical, ehow that the ffavaer- 
Stakes equatioas yield suffici^tly good results for weak shocks 
(M < 2 )j but for strong shocks, (M > z), these results are not 
satisfactory. 

{Ib oYercome this, Bfott-Stoith (1951 ) treated the problem for 
the strong shock case, from the molecular point of view by using 
moments of the Bol t^nann equation. He proposed an ansats for the 
distribution function considering the physics of the problem. Ilie 
shock region being a thin one, be suggested that molecules from 
both the hot and the cold side might freely interpenetrate into 
this region. He assumed the distribution function to be a weight 
sum of two Marwellians, one corresponding to the cold side and the 
other to the hot* (Sae two wei^t functions were determined 
approximately by satisfying the Boltzmann equation in an average 
sense (mommt method). Ihe resiiLts obtained thus were found to b 
satisfactory for a strong shock* Mott-Smith^s analysis was, 
however# restricted to rigid sphere and Sutherland molecules. It 
was further extended to a general inter-molecular inverse»power 
repulsive force law by Muckenfuss (i960). 

Ibr the ease of a diatomic gas with rotational degrees of 
freedom the shock problem was solved by !Ealbot and Seala (l96l) 
from a continuum point of view (loodified Navier*Stokes equatSnns), 
They defined two temperatures separately one for translational 
degrees and the other for rotational ones and incorporated them 
in the Navier*^tokes equations. For small relaxation times they 
proposed a perturbation technique based on a modified bulk viscosi 
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INTRODUCTION 


1,1. This work is concerned with the probleta of shock structure 
in a dissociating diatomic gas. In gas dynamics, it is customary 
to assume the shock to be a surface of discontinuity across which 
two uniform states exist., Ihese two unifoim states are related 
by the Ranlcine - Higoniot relations. Hjwever, to describe the 
narrow transition zoiie in the shock wave it is necessary to consider 
the effects of viscoaity and thermal conductivity. This problem 
was solved within the frame work of the Navier-Stokes equations 
for a monatomic gas first by Taylor (l910) and later by Becker 
( 1921 ), Thomas (1944) and others. The thickness of shock v/ave 
was found to be of the order of a few mean free paths of the ^s 
upstream. But the validity of the Navler~Stokes equations for 
the shock phenomenon is questionable as the characteristic length 
for the shook is so small as to brealc the continuum hypothesis 
made in the formulation of these equations. The results obtained 
by continum and kinetic theory approaches for the shook problem, 
as given by Gilbarg and Daolucci (1953), show roughly equivalent 
shock thicknesses. Ifowever, the merit of a theory ahoxild be 
judged by comparison of its reeiilts with those given by experiments. 
Observations made for monatomic gases by Talbot (1962) and 
Sherman and Talbot (i960), after a careful study of the available 



fldiefreae for lon^ relaxation tdaen they used a eubsidiary relaxatior 
equation to complete the set of equations, ®ie resiats obtained b; 
the bulk viscosity approach were found to be unsatisfactory* 
nevertheless, the relaxation loodel was found to give phenomenologi- 
cally acceptable results which depend on the relaxation time in 
addition to the Mach number* Just as in the monatomic case, in 
this case also, it is reasonable to expect the continuum analysis 
yield results valid for weak shocks only* IJhus Ihe results given 
by the relaxation model are limifed to weak shock cases only* 

Seals and l^llxjt (1965) further extended ihe two temperatur 
relaxation model to a vibrationally relaxing diatomic gas by 
defining one more temperature to include vibrational relaxation 
effects* Bley also used in this case, two subsidiary relaxation 
equations, instead of one, to complete the set of equations. With 
these modifications the Havter-Stokes equations were solved ^ong 
with the two relaxation equations for the shock problem. Out of tl 
two relaxation times, they fixed one, the rotational relaxation 
time, and obtained results for different vibrational relaxation 
times and for two Mach numbers* In this analysis also physically 
acceptable solutions were obtained but the limitation of such a 
(continuum) theory, as mentioned earlier, lies in the fact that it 
IS capable of predicting the structure of only weak shock waves 
satisfactorily. 

Turcotte and Scholniok (1969) proposed a Beam Continuum Mod. 
(b*G.M. ) for monatomic strong shocks in hypersonic flows. In thi£ 
the upstream hypersonic flow is taken to be a zero temperature 
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data, both experimental and theoretical, show that the Uavier- 
Stokes equations yield sufficiently good results for weak shocks 
(M < 2)} but for strong shocks, (m > 2 ), these resvCLts are not 
satisfactory* 

2b overcome this, Mott-anith (l95l) treated the problem for 
the strong shook case, from the molecular point of view by using 
moments of the Boltzmann equation. He proposed an ansate for the 
distribution function considering the physics of the problem. 33ie 
shock region being a thin one, he suggested that molecules from 
both the hot and the cold side might freely interpenetrate into 
this region. He assumed the distribution function to, be a weighted 
sum of two Maxwellians, one corresponding to the cold side and the 
other to the hot. The two weight functions were determined 
approximately by satisfying the Boltzmann equation in an average 
sense (moment method). The results obtained thus were found to be 
satisfactory for a strong shook. Mott-BImith's analysis was, 
however, restricted to rigid sphere and Sutherland molecules. It 
was further extended to a general inter-molecular inverse-power 
repulsive force law by iluckenfuss (i960). 

Bar the case of a diatomic gas with rotational degrees of 
freedom the shock problem was solved by Talbot and Scala (l96l) 
from a continuum point of view (modified Navier-Stokes equations). 
They defined two temperatures separately one for translational 
degrees and the other for rotational ones and incorporated them 
in the Havi er-Stofces equations, Bor small relaxation times they 
proposed a perturbation technique based on a modified bulk viscosity 
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:ontin,um flow is 

Itother, this change is supposed to be bro^ight about 
collision per molecule of the beam, 5he Ifavier-Stokes equations 
are v^ed to treat the continuum particles with additional terns to 
account for liie addition of mass^ momentum and energy from beam 
particles* 53h.e shocS; profile given by this theory is in reasonable 
agreement with some of the experimental profiles given by Oamao 
(1965) and HoWiQa and iCalbot (1966);- -* 

Ohis model (b.C.M.) has been later extended by furcotte (1969) 
to the diatomic case (to include the effects of rotational degrees 
of freedom). “Sie modified JTavier-Stokes equations proposed by 
lEalbot and Scala (1961 ) have been used to treat the continuum flow 
in this case. An equivalent hard sphere model was used to specify 
the interaction cross section vdiich appears in the expressions for 
the transport coefficients. Ihe results were obtained by varying 
the number of collisions required for rotational relaxation, as a 
parameter for two differ^t Mach numbers. They are in good agreemen 
with the experimental results given for strong shocks in nitrogen, 
by Bobben and Talbot (1966). 

1,2. Before dealing with the shock problaa in a diatomic gas from 
a molecular point of view, one has to develop a suitable kinetic 
theory which takes into account the internal degrees of freedom. 
There are two available approaches to the kinetic theory of diatomic 
gases >- the semi‘*quantum mechanical and the classical. 




3 


whereas for long relaxation tijues they used a subsidiary relaxation 
equation to complete the set of equations. The results obtained by 
the bulk -viscosity approach were found to be unsatisfactory. 
Nevertheless^ the relaxation model was found to give phenomenologi- 
cally acceptable results which depend on the relaxation time in 
addition to the Mach number. J-uat as in the monatomic case, in 
this case also, it is reasonable to expect the continuum analysis to 
yield results valid for weak shocks only. Thus 1he results given 
by the relaxation model are limi-ted to weak shock cases only. 

Scala and Talbot (1963) further extended ihe -two temperature 
relaxation model to a vibrationally relaxing diatomic gas by 
defining one more temperature to include vibrational relaxation 
effects. They also used in this case, two subsidiary relaxation 
equations, Instead of one, to complete the set of equations, with 
these modifications the Navier-Stokes equations were solved along 
with the two relaxation equations for the shock problem. Out of the 
two relaxation times, they fixed one, the rotational relaxation 
time, and obtained results for different vibrational relaxation 
times and for two Mach numbers* In this analysis also physically 
acceptable solutions were obtained but the limitation of such a 
(continuum) theory, as mentioned earlier, lies in the fact that It 
is capable of predicting the structure of only weak shook waves 
satisfactorily. 

Turcotte and Schoinick (1969) proposed a Beam Continuum Model 
(b.O.M. ) for monatomic strong shocks in hypersonic flows. In this, 
the upstream hypersonic flow is taken to be a zero temperature 



jLn tile sesai-quantuQi mechanical approach , the classical descity 
■Hon fbr 'the translational dogreos of freedom Is retained, tjut the 
internal degrees of freedom are handled qiiantum mechaaioaliy* The 
difficulties in such a trea-tment lie in the qmntum mechanical 
evaluation of angular dependent in*elastio cross sections, This 
theory was given by Chang and Uhlenbeck (l951 ) end later was 

used by Monchik, Yun and Ifeson (1963) for evaluating the transport 
properties of gases, To enable a treatment of a wider class of 
problems in gas dynaaios, physically plausible kinetic models were 
proposed for a modified Bolt2mam equation, by Morse (1964 ), Brau 
(1967) and Hanson and Morse (1967). These are similar in spirit to 
the B,G,K, mod^ (l954) for the monatamic gas case* 

Further, an analysis close to Grad's thirteen Bionient approxi- 
mation (1949) for the monatomic gas case has bean given by Mb Oormio. 
(196s) to get a set of kinetic equations for poly-atomic gases. He 
used Wang Chang, TIfalenbeck and de Baer equation (1964) (inelastic 
Boltamann equation) and four addi'tional moments to take care of the 
internal temperature and the internal heat-flux vector to obtain the 
seventeen moment approcimation, dhese equations are used for 
studying the near equilibrium and relaxation type of flow. The 
restJlts obtained are in very close agreement with those of Wang Chan 
et al (1964) analysis. 

In the classical approach a suitable model is chosen for the 
gas molecule possessing internal degrees of freedom. The tsro 
body collision problem is solved classically and a modified Boltzman 
equation formulated. Such a theory was developed by Curtiss (l956, 
1957) for diatomic gases having rotational degrees of freedom, and 
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molecular beam and the conversion of this beam into a suhgoale • 
continuum flow is assumed to take place throu^ collisions, 

]rurther, this change is supposed to be brou^t about by a single 
collision per molecule of the beam. The Navier-Stokes equations 
are used to treat the continuum particles with additional terms to 
account for the addition of aass| momentum and energy from beam 
particles* !Qie shock profile given by this theory is in reasonable 
agreement with some of the experimental profiles give4i;by Oamao 
(1965) and Robben and HJalbot (1966);- - . . 

This model (b.O.M, ) has been later extended by Turcotte (1969) 
to the diatomic case (to include the effects of rotational degrees 
of freedom). The modified Navier-Stokes equations proposed by 
Talbot and Scala (1961 ) have been used to treat the continuum flow- 
in this case. An eqxxtvalent ha3:d sphere model was used to specify 
the interaction cross section which appears in the expressions for 
the -feransport coefficients. The resxilts were ob-tained by varying 
the number of collisions required for rotational relaxation, as a 
parameter for two different Mach numbers. They are in good agreement 
with the experimental results given for strong shocks in ni-fcrogen, 
by Robben and Talbot (1966), 

1,2. Before dealing with the shock problem in a diatomic gas from 
a molecular point of view, one has to develop a suitable kinetic 
theory which takes into account the internal degrees of freedom. 

There are two avai4.able approaches to the kinetic theory of diatomic 
gases - the sejai-quantuci mechanical and the olassioal. 



ms fiirthor xised for evaluating the transport pro parti os fer gases 
with loaded sphere and sphero-cylinder laolecular laodels by I>alher 
and father (l965) mid Sandler (l966, 1968). 

1,5* first attaa.pt "to obtain the shock structure in a 

rotationally relaxing gas from the molecular point of view was made 
by Haight and iundgren (l965t 1967^ 1968), They used the kinetic 
theory developed by a classical approach, for the strong shock. A 
tximodai distribution function, built up using three Maxwellians, 
describing the te?o uniform states - the upstream and the down stre£ 
and an intermediate state having no rotational excitation - was 
proposed. A moment method was used to determine the three weight 
functions. This analysis is an extension of Ifott-Smith’s bimodal 
theory (l95l). Shock profiles for three different cases of 
molecular models (rough sphere, loaded sphere and sphero-cylinder) 
were obtained. As a special case of this theory, it ves possible 
to recover Mctt-Snitfa monatomic results, 

later the trimodal ansats was used by Brau, Simons and Macon 
( 1969 ) in solving the kinetic model equation proposed by Brau (196'" 
vath the shock boundary conditions. They also tried an, alternate 
ansatz based on a model developed by toderson and Efecomber (l965) 
for monatomic gases. While the trimodal theory gave good results 
for both low and hi^ upstream Ifech numbers agreeing with the 
experimental data of Eobben and Talbot (1966) for shock profiles ar 
shock thickness, the second ansatz yielded reasonably good results 
only for shock thickness. Uie shock profiles predicted by the ser 
ansatz were not in agreement with the experimental shock profiles : 
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In the semi-quantum mechanical approach, the classical descrip- 
tloa for the translational degrees of freedom is retained; but the 
internal degrees of freedom are handled quantum mechanically. The 
difficulties in such a treatment lie in the quantum mectenical 
evaluation of angular dependent in-elastic cross sections. This 
theoiy was given by Wang Chang and tJhlenbeok (1951 ) and later was 
used by Monchik, Yun and Mason (1965) for evaluating the transport 
properties of gases. To enable a treatment of a wider class of 
problems in gas dynamics, physically plausible kinetic models were 
proposed for a modified Boltzmann equation, by Morse (1964- ), Brau 
(196?) and Hanson and Morse (l967)« These are similar in spirit to 
the B,G,K. model (1954) for the monatomic gas case, 

]?urther, an analysis close to Grad’s thirteen moment approxi- 
mation (1949) for the monatomic gas case has been given by Mo Oormlok 
(196s) to get a set of kinetic equations for poly-atomic gases. He 
used Wang Chang, Uhlenbeck and de Boer equation (l964) (inelastic 
Boltzmann equation) and four additional moments to take care of the 
internal temperature and the internal heat-flux vector to obtain the 
seventeen moment approximation. !Ihese equations are used for 
studying the near equilibrium and relaxation type of flow. The 
results obtained are in very close agreement with those of Wang Chang 
et al (1964) analysis. 

In the olasaical approach a suitable model is chosen for the 
gas molecule possessing internal degrees of freedom. The two 
body collision problem is solved classically and a modified Boltzmann 
equation formulated. Such a theory was developed by Curtiss (l956, 
1957) for diatomic gases having rotational degrees of freedom, and 



high upatreaa Mach numbers. 

Venkatramon and Morse (1969) used the kinetic model gavan by 
Moree (l964)y to study the structure of a shock wave in diatomic 
gases with rota-y.onal degrees of freedom. Ihey assumed a suitably 
modified Mott«3mith bimodal distribution. In this ansatz, in additic, 
to the two weight fuiiffl tions , the parameters of the second Jiaxwellian 
were also taken to be functions of the distance, {Qais was necessary 
in order to satisfy the energy relaxation equation. Hie three 
conservation equations and two moment equations ( and 1/2 me - Ej^) 
were used to determine the five unknowns in the problem. Comparison 
of density and temperature profiles with the experimental data of 
Robben and Talbot (1966) is quite ©eod. 

Recently ]\fe.chperson (l97l) followed the classical approach 
and a ifonte Carlo scheme for obtaining the rotational taaperature 
profiles for shock waves in oxygen and nitrogen* Ife also calculated 
the variation of translational temperature and density inside the 
shock. The intermolecular potential used in this analysis is the 
one given by Parker (l959)« Hie shock thicknesses predicted thus, 
compare well with the experimental values obtained by Mairone (l967). 

Venkataraman and Morse (l970) extended their earlier analysis 
(1969) to F^dict the effects of vibrational relaxation on the 
shock 'Structure In dis'fcomtc gas, A kinetio model, sijallar to the 
one given by Morse (1964), with additional terms to include ■vibration 
effects, was also proposed by the®. In addition to the three 

2 1 2 N 

conservation equations, they used c^,( sQC *“ ®vi' ^^oments of 

the kinetic equation. Shock structures in oxygen, nitrogen and 
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was further used for e-vaiuating the transport properties for gases 
with loaded sphere and sphere *oylinder molecular models by Dalher 
and fiather (l963) and Sandler (1966, 1968). 

1 ,3. The first attoapt to obtain the shock structure in a 
rotationally relaxing gas from the molecular point of view was made 
by Haight and Lundgren (1965» 1967, 1968), They used the kinetic 
theory developed by a classical approach, for the strong shock, A 
trimodai distribution function, built up using three Maxwellians, 
describing the two xaniform states - the upstream and the down stream 
and an intermediate state having no rotational excitation - was 
proposed, A moment method was used to determine the three weight 
functions. This analysis is an extension of Mott-Saith*s blniodal 
theory ( 1951 ), Shock profiles for throe different oases of 
molecular models (rough sphere, loaded sphere and sphero-oylinder) 
ware obtained. As a special case of this theory, it «as possible 
to recover Mott-Smith monatomic results, 

later the trimodai ansatz was used by Brau, Simons and Maoomber 
(1969) in solving the kinetic model equation proposed by Brau (l967) 
Vvrith the shock boundary conditions. They also tried an alternate 
ansata based on a model developed by Anderson and Maoomber (1965) 
for nonatomie gases. While the trimodai theory gave good results 
for both low and high upstream ASach numbers agreeing with the 
experimental data of Ho bben and Talbot (l 966) for shock profiles and 
shook thickness, the second ansatz yielded reasonably good results 
only for shook thickness. The shook profiles predicted by the second 
ansatz were not in agreement with the experimental shock profiles for 



carbon di-oxxde were obtained, for high Mach ntaabers using this 
nod el. 

1*4» In the present woik, an axially elastic sphero --cylinder model 
IS proposed for a diatomic gae molecule. Ajcial elasticity has been 
introduced to account for the vibrational degrees of freedom. Ihe 
distribution function in such a case depends on x, 
and <f» . ffiie last two of the variables are required to describe the 
vibrational state of a molecule. A mathematical description of the 
molecijle is required for working out the dynamics of a collision 
between two laolecules. Ihe mathematical representation followed for 
the axially elastic vibrating molecule follows closely, the one used 
by Kaight and Iiundgren {l965, 1967, 1968) for the rigid non-central 
particles. 

21ie geometry and dynamics of a collision between two molecxfLes 
are studied at first. IJhe collision integrals In the modified Boltaf 
equation are formulated. The moment equation corresponding to the 
Maxwell transfer equation la obtained. As special cases of this, tt 
three conservation equations are written down. iChrough an H-theoren 
the corresponding equilibrium distribution is obtained. 

Purther for investigating the shock structure in a vibrational! 
and rotationaUy relaxing diatomic gas, the moment method is used, 
ibllovdng Efott-Smith (1951 ) and ihiight (l965) two multi-snodAl ansatz. 
for the distribution function have be^x proposed - the first one a tri 
using three Maxwellians, the second one, a qtiadrimodal using fotir 
Maxwellians. In the trlmodal case a molecule inside the shock is 
assumed to belong either to the upstream eqxiilibrium state, the 
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high upstream Mach numbers. 

Venkatraman and Morse (1969) used the kinetic model given by 
Morse (l964), to study the structure of a shock wave in diatomic 
gases with rotational degrees of freedom. They assumed a suitably 
modified Mott-Smith bimodal distribution. In this ansatz, in addition 
to the two weight functions, the parameters of the second Maxwellian 
were also taken to be functions of the distance. Ihis was necessary 
in order to satisfy the energy relaxation equation. !Phe three 
conservation equations and two moment equations ( and 1/2 mcJ ~ Ej^) 
7/ere vised to determine the five unknovwia in the problem. Comparison 
of density and temperature profiles vdth the experimental data of 
Robben and Talbot (1966) is quite good. 

Recently Machperson (l97l) followed the classical approach 
and a Mbnte Carlo schtme for obtaining the rotational temperature 
profiles for shock waves in oxygen and nitrogen. He also oaloulated 
the variation of translational temperature and density inside the 
shock. The intermolecular potential used in this analysis is the 
one given by Parker (1959)» 53ie shock thicknesses predicted thus, 
compare well with the ejqjerimental values obtained by Marrone (1967). 

Venkataraman and Morse ( 1970 ) extended their earlier analysis 
( 1969 ) to pafiedict the effects of vibrational relaxation on the 
shock •structure in diatonic gas. Akinetic model, Blrullar •fco •the 
one given by Morse (1964), with additional terms •to include vibrational 
effects, was also proposed by them. In addition •to the three 
conservation equations, th^ used "I: ) moms of 

the kinetic equation. Shock S true tures in oxygen, nitrogen and 
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iownaia^oaia . M^wellian state or an intermediate etjuilibrivnn state 

without vibrational excitations. In the quadrimodal case, in 

addition to these three Maxwellians, atKither ISaxwelliaa^ which 

corresponds to an intermediate equilibrium state with only translatior 

degrees of freedom fully excited and rotational and vibrational 

degrees of freedom frozen^ is taken. !Ilie three conservation equatior 

effectively give one algebraic equation in the unknown weight 

f^inction^, To determine three weight functions in the trimodai 

formulation, Iwo moment equations - rotational energy and vibrational 

energy - are used. fbr the quadrimodal approach the fourth equation 

2 

used is 0 moment of the Boltzmann equation. The collision integral? 
in these moment equations are 22 fold, out of which 16 fold integra- 
tions are carried out analytically and the remaining Integrations 
were done nuiaeriG£0.1y using the &bnte Carlo technique. Effectively 
we obtssin a set of two non-linear, first order, coupled total 
differential eqmtionsand an algebraic equation in the three unknown 

weight functions for the txiEio.dal case vdiereas for the quadrimodal 

2 

case an .additional differential equation is provided by the cr 
moment equation. These equations are solved nxmierically on a 
computer and shock profiles are obtained for various Hach numbers and 
upstream pressures for nitrogen. 

1,5» As the Mach number increases the degree of dissociation in a 
diatomic gas behind the shock wave becomes increasingly larger* 

Thus 1he theory proposed above for the shock structure is valid at or. 
moderately hi^ Mach numbers and it is necessary to include the effec 
of dissociation in the analysis of the shock structure at higher Ifech 


numbers 
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oarbon di-oxide were obtained for high Maoh numbers using this 
model 4 

1.4* In the present work, an axially elastic sphero-Kiyliader model 
is proposed for a diatomic gas molecule. Axial elasticity has been 
introduced to account for the vibrational degrees of freedom. !Dhe 
distribution function in such a case depends on x , c . e , w , B 
and <j) . The last two of the variables are required to describe the 
vibrational state of a molecule. A mathematical description of the 
molecule is required for working out the dynamics of a collision 
between two raolecules. The mtheraatical representation followed for 
the axially elastic vibrating molecule follows closely, the one used 
by ftaight and Iiundgren (1965, 1967, 1968) for the rigid non-central 
particles. 

Tlie geometry and dynamics of a collision between two molecules 
are studied at first. The collision integrals in the modified Boltzmann 
equation are formulated. The moment equation corresponding to the 
Maxwell transfer equation is obtained. As special cases, of this, the 
three cons oiwation equations are writtexi down. Through an H-theorem, 
the corresponding equilibrium distribution is obtained. 

further for investigating the shock structure in a vlbrationaliy 
and rotationally relaxing diatomic gas, the moment method is used, 
following Mb tt~Smith (1951 ) and Hhight (l965) two multi-miodftl ansatzt 
for the distribution function have been proposed - the first one a trimodal, 
using three Maxwellians, the second, one, a quadrimodal using four 
Maxwellians. In the trimodal case a molecule inside the shock Is 
assumed to belong either to the .upstream eqtiilibrium state, the 
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2!he analysis outlined, above has to be modified in two ways; 
firstly because of the appearance of a second species, the atom, due 
to the onset of molecular dissociation, an additional Boltzmann. 
eq^^ation has to be introduced governing the distribution function of 
the atom species j secondly, additional collisional terms appear on 
the right hand side of iiie Boltamann equation for the niatomic species 
representing dissociative and reoombinatlve collisions corresponding 
to atom - E»lecule collisions# dihus the problem reduces to the 
calcxilation of the stoclc structure in a ohOTdoally reacting ntbcture 
of two Species, atoms and molecules# In this case also, a inoment 
method is used to solve the Boltemann equations for two species. 32 iq 
trinodal ansats suggested earlier is adopted for the diatomic species 
and the Jibtt-Sinith bimodal distribution is used for the atom species. 
Two species conservation equations and tm moments of the Soltmana 
equation for diatomic species (rotational energy and vibrational 
energy) are lised to determine the unknown weight functions, She 
shock pix) files are obtained for a range of upstream Ifeoh numbers and 
apetream pressures in nitrogen# 

In chapter 2 the collision dynamics of two vibrating noleeules 
is solved. In Chapter 3 is given the fomulatlon of Boltaaam equati 
for a diatomic and also the equilibrium solution for this equatic 
fhe tcimodal and quadrimodal analysis for the structure of a shock 
wave in a vlbrationaHy relaxing diatomic ^s has been dealt in 
Chapter 4. fhe dissociation and recombination rates are written foi 
intemaol ecular and molecwle-etoa collisions in Chapter 5# In this 
chapter, the formiilatioa for the shock structure problem con^dering 
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downato'oam . Btewellian state or an intermediate equilibrium state 

without vibrational excitatione. In the quadrimodal case, in 

addition to these three Max wellians, another Maxwellian, which 

corresponds to an intermediate equilibrium state with only translational 

degrees of freedom fully excited and rotational and vibrational 

degrees of freedom frozen, is taken* Ihe three conservation equations 

effectively give one algebraic equation in the unknown weight 

functions. Ib determine three weight functions in the trimodal 

formulation, two moment equations - rotational energy and vibrational 

energy - are used* Bbr the quadrimodal approach the fourth equation 
2 

used is 0 ^ moment of the Boltzmann equation. The collision integrals 
in these moment equations are 22 fold, out of which 16 fold integra- 
tions are carried out analytically and the remaining Integrations 
were done numerically using the Monte Carlo technique. Effectively 
we obtain a set of two non-linear, first order, coupled total 
differential equations and an algebraic equation in the three unknown 

weight functions for the trimodal case whereas for the quadrimodal 

2 

case an additional differential equation is provided by the c 
laoment equation. These equations are solved numerically on a 
computer and shock profiles are obtained for various Mach numbers and 
upstream pressures for nitrogen. 

1*5, As the Mach number increases the, degree of dissociation in a 
diatomic gas behind, the shock wave becomes increasingly larger. 

Thus Ihe theory proposed above for the shook structure is valid at only 
moderately high Mach numbers and it is necessary to include the effects 
of dissociation in the analysis of the shock structure at higher Maoh 
numbers . 
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the dissoQiation effects is also given. In Chapter 6 is discussed 
the nunerical teciunique (Monte Carlo method) used for evaluating the 
multiple integrals in the moment equations along with the error 
estimates. linally in, Chapter 7 are given results and conclusions 
of the present investigations. 
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Sie analysis outlined above has to be modified in two ways: 
firstly because of the appearance of a second species, the atom, due 
to the onset of molecular dissociation, an additional Boltemann 
equation has to be introduced governing the distribution function of 
the atom species; secondly, additional collisional terms appear on 
the right hand side of the Bolt 2 a 3 iann equation for the clatomic species 
representing dissociative and recombinative collisions corresponding 
to atom - molecvile collisions. Thus the problem reduces to the 
calculation of the shock structure in a chemically reacting mixture 
of tv/o species, atoms and molecules. In this case also, a moment 
method is used, to solve the Boltzmann equations for tvAO species. The 
trimodal ansatz suggested earlier is adopted for the diatomic species 
and the Mott~Smlth blmodal distribution is used for the atom species. 
Two species conservation equations and two moments of the Boltzmann 
eqmtion for diatomic species (i^otational energy and vibrational 
energy) are used to determine the unknown weight functions. The 
shock profiles are obtained for a range of upstream JkEioh. numbers and 
tipstream pressures in nitrogen. 

In Chapter 2 the collision dynamics of two vibrating .nolocules 
is solved. In Chapter 3 is given the formulation of Boltzmann equation 
fora diatomic gas and also the equilibrium solution for this equation. 
The trimodal and quadrimodal analysis for the structure of a shock 
wave in a vibrationally relaxing diatomic ffiB has been dealt in 
Chapter 4. The dissociation and recombination rates are written for 
intermolecular and molocule'-a tom collisions in Chapter 5* Ih this 
chapter, the formulation for the shock structure problem con aider jng 
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coiusioN nmsmoB 


2.1. MoleciJlar Mjdel t 

Before developing a kinetio theory for a vibrationally, 
rotationally and translationally relaxing gas, a suitable molecvaar 
model has to be postulated* Jeans (1925), proposed a rigid loaded 
sphere nodel for rotationally relaxing diatomic gas molecules. later, 
Hirsehfelder, Ciirtiss and Bird (1954) introduced a rigid spherocylinder 
model in the kinetic theory of diatomic gases. However, one of the 
first few attempts in developing a kinetic theory of non-spherical 
molecules, with rotational degrees of frtiedom by Gxirtiss (l956), used 
a general class of molecules which are rigid convex bodies of revolution. 
The loaded sphere model is a special case of this general model and the 
spherocylinder, in some sense, also belongs to the general class of 
convex bodies of revolution* 

With the available models the effects of rotational relaxation 
ofl transport properties in a diatonic gas, have been estimated by 
Gather (1963) and &ndler (1966,1968) and on the shock structure by 
feight (1965) in such a gas, A classica.1 theory of vibrationally 
relaxing diatonic gas is not available, so far in the literature. 

Por this purpose, a molecular model which permits the transfer of 
vibrational energy in addition to the rotational and translational 
energies at a collision has to be chosen at first* 
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the dissociation effects is also given. In Chapter 6 is discussed 
the numerical technique (Monte Carlo method) used for evaluating the 
multiple integrals in the moment equations along with the error 
estimates. Elnally in Chapter 7 are given results and conclusions 
of the present investigations. 
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In tills thesis an axially elastic sphere-cylinder model is 
proposed for such a gas. The axial elasticity is introduced for 
the central cylindrical portion to allow the molecules vibrational 
degrees of freedom. This model can be regarded as composed of an 
assembly of tm rigid spherical atoms connected by a linear elastic 
spring guided in a rigid hollow cylindrical shell as shown in 
Iig.(2,1 *1), 

The vibrational degrees, as indicated earlier are treated class! 
The justification for this is given in Chapter 3. The molecule Is 
assmed to vibrate from an initial equilibrium state. The amplitude 
of the exial vibratdnn is taken to be small compared to the eqvolibriUi 
separation distance between the two fitoms constituting the molecule. 
The molecule is assumed to be rigid in the transverse direction, so as 
not to excite any bending vibrations in that direction, Jbr a given 
homogeneous species of molecules the frequency of vibrations of a 
molecule is the sane for all the molecules. However, the two varlabl 
used to describe the vibrational state of a molecule are; the maximum 
amplitude, B, of vibration measured from the equilibrium positions of 
two atoms and the phase angle, ^ , of the simple harmonic vibrations 

of an individual molecule. The amplitude B is a measure of the vibra 
tional energy of a molecule. 

In addition 1 d the mss and the two geometrical parameters ( a 
and b) required to specify a rigid sphere -cylinder model, the axially 
elastic spherocylinder model also needs a spring constant, s, for its 
complete description. As the analysis is restricted to small amplitui 
of vibration, a linear spring is an adequate representation. The two 
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COIMSIOK PmAMICS 


2*1. Molecular Model ; 

Before developing a kinetic theory for a vibrationally, 
rotationnlly and translationally relaxing gas^ a suitable molecular 
model has to be postulated. Jeans (l925), proposed a rigid loaded 
sphere j’jodel for rotationally relaxing diatomic gas molecules. later, 
HLrschf elder, Curtiss and Bird (l954) introduced a rigid sphexo cylinder 
model in the kinetic theory of diatomic gases. However, one of the 
first few attempts in developing a kinetic theory of non-spherlcal 
molecules, with rotational degrees of freedom by Curtiss (l956), used 
a general class of molecules which are rigid convex bodies of revolution. 
The loaded sphere model is a special case of this general model and the 
spherocylinder, in some sense, also belongs to the general class of 
convex bodies of revolution. 

With the available models the effects of rotational relaxation 
oiX transport properties in a diatomic gas, have been estimated by 
Sather (l963) and Sandler (1966,1968) and on the shock structiore by 
Efeight (1965) in such a gas. A classical theory of vibrationally 
relaxing diatonic gas is not available, so far in the literature. 

Bor this purpose, a molecular model whioh permits the transfer of 
vibrational energy in addition to the axjtational'and translational 
energies at a collision has to be chosen at first. 
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T^iriables B and, d y describing the vibrational state of a noleenle enter 
not only in the analysis of collision dynanica, but also In the 
definition of a distribution function for the gas, Glie iciportant 
physical constants for some diatoniic molecules are stovii in Table 

(2.1.1 ). 


Surface Representation 


The nolecular surface is described in terns of a unit vector 
e along the longitudinal axis of the molecule and a unit vector n 
normal to the surface » such that the location vector r for a point 

•w 

on the surface from the centre of mass of the molecule may be written 
as follows (See lig. 2. 2*1.} t 


r = {b + B cos (vt + <j)) } e + a n ) for points on the 

r • n = {b + B cos (vt + <*>)} e*n+aj positive e hemisphere 

4# ftr 


r = - {b + B CDS (vt + 4t)>e + ap 1 for points on' the 
r * n = - {b + E cos (vt + <&) } e-n + a } ^aegative e hemisphere 


(2.2.1) 


r»n = a , e<>n = 0 


for points on Idle cylindrical 
surf ac e 


Bynaolcs of a Two Body Collision 


In deriving the Boltzmann equation, we shall consider a dilute 
gas and neglect the effects of ternery and higher order collisions, 
as in the case of a monatomic gas. Thus we need solve only a binaiy 
collision problem. The collisional contact time for vibrating 
molecules is finite and is appreciably large compared to the period of 
vibration (Stupochorto, 1966 ) . However ^ in this analysis we assume 
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In this thesis an axially elastic sphero-cylinder model is 
profiosed for such a gas. The axial elasticity is introduced for 
the central cylindrical portion to allov/ the molecules vibrational 
degrees of freedom. This Eiodel can be regarded as composed of an 
assembly of tvro rigid spherical atoms connected by a linear elastic 
spring guided in a rigid hollow cylindrical shell as shown in 
Sig,(2.l4l). 

The vibrational degrees, as indicated earlier are treated olassioally. 
The justification for this is given in Chapter 3. The molecule is 
assumed to vibrate from an initial equilibrium state. The amplitude 
of the axial vibration is taken to be small compared to the equilibrium 
separation distance between the two »toms constituting the molecule. 

The molecule is assumed to be rigid in the transverse direction, so as 
not to excite any bending vibrations in that direction, Ibr a given 
homogeneous species of molecules the frequency of vibrations of a 
molecule is the sane for all the molecules. However, the two variables 
used to describe the vibrational state of a molecule are: the maximum 
amplitude, B, of vibration measured from the equilibrium positions of 
two atoras and the phase angle, (|) , of the simple liarmonlc vibrations 

of an individual molecxfle. The amplitude B is a measxore of the vibraw 
tional energy of a molecule. 

In addition to the mass and the two geometrical parameters ( a 
and b) required to specify a rigid sphere "Cylinder model, the axially 
elastic sphero cylinder model also needs a spring constant, s, for its 
complete description. As the analysis is restricted to small amplitudes 
of vibration, a linear spring is an adequate repfesentatton. The t wo 
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TABLE 2 . 1 .: 

1 



Molecule 

^2 

'^2 

^2 

Atomic mass nuraber 

1 

14 

16 

d 

Inter-atomic distance re(A) 

0.7416 

1.094 

1.2074 

Dissociation Energy D(eY) 

4.476 

9.T4 

5.11 

Dissociation temperature - D/k(^K) - 

Characteristic "Vibration temperature 

115000 

59380 

If 

5910 

5340 

2250 

Characteristic rotation 

temperature (®K) 

85.4 

2.9 

2.1 

a/h 

5.8 

3.1* 


v2 

mb 

2J^ 


0 . 025 ^ 



* Galculations are made for a/l* - 5 a^ad, » l/52 ard tliese p/ce 

Zo ^ 

close to recommended values for nitrogen molecules dy Haigiit, 
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v^-iriables B and <> , describing the vibrational state of a nolecule enter 
not only in the analysis of collision dynamics, but also in the 
definition of a distribution function for the gas. The important 
physical constants for some diatomic molecules are shovai in Table 

(2.1.1). 


2.2. Surf ao e R epr es an ta tion : 


The molecular surface is described in terms of a unit vector 
e along the longitudinal axis of the molecule and a lAnit vector n 

m w 

normal to the surface, such that the location vector r for a point 

le 

on the surface from the centre of mass of the molecule may be written 
as follows (See iig, 2.2.1.) : 


r a {b + B cos (vt + <i))} e + a n ^ for points on the 

r • n =» {b + B cos (vt + «}»)} e*n+a3 positive e hemisphere 

r = - {b + B cos (vt + (|))} e + a n 1 for points oa' the (2,2.1 ) 

r .n - {b + B cos (vt + d.)} e»n + a \ negative e hemisphere 


r'n = a , e»n=0 

^ J ww ** 


for points on the cylindrical 
surface 


2.3. Dynamics of a Two Body Collision s 

In deriving the Boltzmann equation, we shall consider a dilute 
gas and neglect the effects of ternary and higher order collisions, 
as in the case of a monatomic gas. Thus we need solve only a binary 
collision problem, ahe collisionai contact tine for vibrating 
molecules is finite and is appreciably large compared to the period of 
vi.bration (stupochonto, 1966 ) * Ifcwever, in this analysis we assume 
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the collisioQS to be instantaneous. At a collision an impulse I 
is transnlttcd from molecule 1 to molecule 2 in the direction of the 
coiaoon normal n at the point of contact (lig. 2*3*1 )* line of 

«p 

action of the impulse does not always pass through the centres of 
mass of the molecules. iThusj in addition to the transfer of linear 
inpulse, an angular impulse Is also transmitted at a collision, 
Further, the simple harmonic oscillations absorb some energy due to 
the Impulse, From the conservation of linear momentum we have, 

In 


c’ 

“1 




In 

s C,2 + — 


for laolecule 1 
for molecule 2 


Ihe angular momentum conservation gives us 


' -.-I- 

m. = (d-t - J- « in X r, 

•'i 1 ** 

(2*3.2) 

w' = Wo » In X 

^4 


The energy absorbed by the vibrational degrees of freedom for molecule 
1 and 2 may be written as 


C&V,B.)j = - (1 n) • Cgj V Bj sin (}>i) 

(iV,E,)2 = 2) ' C®2 ^ ^2 ^2^ 

Finally we have from the conservation of energy, 


4”ni * ?“lUi«i ♦ isC2Bj)h ♦ {ijij * \'i2‘l2-«2 

» |s(2»2)^) ' { 7 ®! ♦ j- el -il • “1 * + (4VE) 

* 4 "2 S2 * I al • ^2> si * T + (AVEip (a.3 
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miiE 2,i»i 


... 

Molecule 

'^2 

'^2 

O2 

Atomic mass hiiralDer 

1 

14 

16 

o 

Inter-atomic distance re(A) 

0,7416 

1.094 

1,2074 

Dissociation Energy D(eV) 

4.476 

9.74 

5.11 

Dissociation Temperature » d/ 1 c(°K) 


113000 

59380 

Characteristic Vibration temperature 



= hv/k(®lO 

5910 

3340 

2230 

Characteristic rotation 

temperature ("^K) 

85.4 

2.9 

2,1 

a/b 

3.8 

3 « 1 * 

»«*• 

2 



.V 

mb 

2J^ 

Mm- 

0,025* 


* Calculations are made for a/b « 

•7 a rab^ 

. 3 oM jj- 

m 1/52 and 

tliese orce 


close to recotamended values for nitrogen molecule a by Haight, 
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Substituting fron Sqns . (2,3.1.) to (2.3.3.) in Eqn. (2.3.4* )j 
v/e have for the iiapiiLsG I : 


I - (- g.ii)/{i 




2 J 


Cr2 X 

+ } 


11 


2 J 


12 


(2.3.5) 


where 


g = ($2 + U 2 ^ ^2 - 22 ^2 ^ ‘*’ 2 ^ 

( 2 . 3 . 6 ) 

* ^Sl ^ i!)l ^ ^1 " tl ^1 ^ '^l) 

Once I is knom in terns of the pre-collision variables and the 

collision geonotry, the post collision variables can be deteminei 

fron the Eqns. (2.3*1.) to (2.3.3*). ^hus binary collision 

problea is conpletely solved for vibrating molecules. 

In the fomulation of Boltanann equation for dissociating 
diatoxrlc gas (see Chapter 5)j we need to consider the collisions 
between an aton and a molecule. Eie post collision variables for 
such a colli pion can be obtained as a special case of the above by 
putting for particle (atom ) 

1^2 = — f - Of <^2 - a)2 = 0 (2.5*7) 

2 * • 


In Chapter 3j a Boltan^nn equation for a diatomic gas consisting of 
elastic spherocylinder molecules is formulated and an equilibrium 
solution is obtained for this equation. 
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the collisions to be instantaneous. At a collision an impulse I 
is transnitted from molecule 1 to molecule 2 in the direction of the 
comon normal n at the point of contact (lig. Bio line of 

action of the impulse does not always pass through the centres of 
mass of the molecules. Thus, in addition to the transfer of linear 
impulse, an angular impulse is also transmitted at a collision, 
Shrther, the simple harmonic oscillations absorb some energy due to 
the impulse, From the conservation of linear momentum we have. 

In 


0 

«* 


! 

1 



£i 


^2 


% 

In 


for molecule I 

( 2 * 3 * 1 ) 

for molecule 2 


The angular momentura conservation gives ua 


0)1 = ui •* 0 In X r, 

~1 - -4l 

( 2 . 3 . 2 ) 

0)' = Wo + • In X r-j 

The energy absorbed by the vibrational degrees of freedom for molecule 
1 and 2 nay be written as 


= - (I n) * (©1 V sin ({)j) 
(AV.E .)2 = (1 n) • (©2 B 2 sin ^ 2 ) 

Finally we have from the conservation of energy, 


{^rajc^ + + is(2Bj)^} + {|m2 02 + |-o)2*£2*W2 

+ - 5(2^2)^} » I- !!Jl -il • I sC2Bi)^ + CAVE)^} 

^ 4^2 2^-^ I -^(2B2)^ y (2.3.4) 




^ k jffrbcLt.'l- Ct- ^ 
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3.1- Introduction - In this chapter a kinetic theory of vibratiomlly 
relaxing diatomic gas is given. As the end result of this thesis is not 
hinetio theory per se f a less general view is taken in formulating the 
Boltainann equation. This is especially so in treating the vibrational 
degrees of freedom, classically. Strictly speaking, the vibrational 
degrees of freedom exhibit quotum behaviour and the vibrational ener^ 
of H identical oscillators (diatomic molecules) is given by Planck’s 
function t 


^vib 


- »TVT f hv/kT 


-} 

1 


(3.1.1) 


vdiere u is the frequency of vibration and T is the translational 
temperature of the gas* In the limit kT >>hv the vibrational 
energy approaches its classical value KkT, Actually the energy 

is already close to its classical value, when kT -v hv * Thus, for 
example, at kT/hv = 0,5j we find 0.724 KkT; at kT/hy * 1 , we 

fxnd 0.928 HkT; at kf/hy =2, v/e find 0.979 llkT. 

The value of hv/k for oxygen is 2230®K and for nitrogen is 
3340''K, (See Table 2.1,1.). 

Ibr a shock wave in a diatomic gas there exists a narrow region 
(translational shock) in v/hich vibrational excitations are practically 
frozen and at the tall Of this region the value of kT/hv >1, 



Substituting from Eqns. (2,3*1*) to (2.3*34) in Eqn. (2.5.4. )» 
we have for the impulse I s 
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I = (- g.ii3/{i (~ 

-ci IHi 


1 , X 

+ — ) + 




'2 


2 J 


(ro X n)" 

+ " } 


11 


2 J 


12 


(2.3.5) 


where 

g = (£2 + ^2 X 52 - 22 B2 V sin <1.2) 

(2.3.6) 

“ ^Si -1 ^ 2i “ Si 

Once I is known in terns of the pre-oollision variables and the 
collision geometry, the post collision variables can be determined 
from the Eqns. (2.3.1.) to (2.3.3.)* %us binary collision 
problem is completely solved for vibrating molecules. 

In the formulation of Boltzmann equation for dissociating 
diatoiaio gas (see Chapter 5), we need to consider the collisions 
between an atom and a molecule. Ihe post collision variables for 
such a collision can be obtained as a special case of the above by 
putting for particle -2 (atom ) 

® — » B2 = 0, 4»2 = 0, and u)2 » 0 (2.3.7) 

In Chapter 3, a Boltzmann equation for a diatomic gas consisting of 
elastic spherocylinder molecules is formulated and an equilibrium 
solution is obtained for this equation. 
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^Hiere^ore* over a significant portion of the shock, following the 
translational shock, the vibrational degrees behave almost classically. 
Henco in the kinetic theory developed here, all the degrees of freedom 
inoluding the vibrational ones are treated classically , 

3.2. Bolt amann J Egnation ~ 

lo describe the dynamic state of an axially elastic sphere- 
cylinder laDlecule, a set of 15 coordinates is required, 

Ihe set consists of : 

X - the physical co-ordinates of the centre of mass of the 
molecule 

o - the translational velocity of the centre of mass of the molecxile 
with respect to the laboratory co-oi‘dinate system 
0 - the unit vector along the axis of the molecule to fix the 

orientation 

[i) - the angular velocity of the molecule with respect to the 

laboratory co-ordinate system. 

B - the maximm amplitude of vibration 

A - the phase angle of the simple harmonic oscillations of the two 
atoms constituting the BKjlecule. 

®ius for a gas mad© up of such molecules, the distribution function is a 
function varying In the hyper-space q , t where q denotes the space 
defined by ( x, e, 2he probable number of raoleculee with 

mit- ^ 

their phase coordioates between q and q at any time t , is given 

by 


f(q;t} dq 



chapter 


3 


BOLa)Zl!.lAMH EQUATION AEP) THE H-THEO RHvI 


3.1. Introduction - In this chapter a kinetic theory of vibrationally 
relaxing diatomic gae is given, Aa the end result of this thesis is not 
kinetic theory per se , a leas general view is taken in formulating the 
Boltzmann equation. This is especially so in treating the vibrational 
degrees of freedom classically, . Strictly speaking, the vibrational 
degrees of freedom exhibit quantvim behaviour and the vibrational energy 
of N identical oscillators (diatomic molecules) is given by Hanck’s 
function: 


^vib 


= ffkT } 


( 3 * 1 * 1 ) 


v/here v is the frequency of vibration and T is the translational 
temperatirre of the gas. In the limit kT >>hv the vibrational 
energy' approaches its classical value llkT, Actually the energy 

is already close to its classical value, when kT '\i hv . Thus, for 
example, at kT/hu = 0.5, we find 0.724 HkTj at kT/hv » 1 , we 

find 0,928 MkTj at kT/hv = 2, we find ®* 0.979 HkT. 

The value of hv/k for oxygen is 2230®K and for nitrogen is 
3340°K. (See Table 2.1.I.). 

ibr a shoclc wave in a diatomic gas there exists a narrow region 
(translational shock) in which vibrational ecxcitations are practloally 
frozen and at the tail of this region the value of kT/hv > 1. 
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or 


£(Xj,c,e,U),B,([>;t)dx dc de dcj dE do 

where £ is the distribution function describing the state of the 
gas and 


dx - dx dy dz j dc = dCx dcy do^ 


de = sin 3 dE dc ; dto = dtoi 

*r « X 


da)2 dw-j . 


(3.2.1 ) 


Shus the number density of the molecules at a point x In the physical 
space and at a given instant t is given by the normalisation relation- 
ship; 


nCx^t) = / f Cx,c,s,a>,B,<{);t) dx dc de dm dD d* (3*2«2) 

Hhe variation of f with time t, is now studied to obtain a 
governing equation for f . In the absence of colllsional effects, 
the tiiae derivatives q » of the I 3 phase coordinates of a molecule 
are functions of the coordinates themselves. After an interval of 
time dt from the instant t, a free molecule initially having q as 
its co-ordinate will occupy a point q + q dt in the hyper-space. Now 

« >v 

consider a set of molecules in a phase volume dq at time t around 
a point q . At time t + dt the same molecules will be present in 

•r 

a phase volume dq/ around the point q + d dt, if they suffer no 

^ « at 

collisions during the interval dt. However due to colllsional effects, 
a fraction of the tolal number of molecules in dq , fall to reach dq' 
while some molecules which were not present in dq at time t, manage 
to reach dq' , The total number of particles in dq' , nt time t+dt, 
is related to the number in dq at t by the following equation. 
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Therefore, over a significanl: porljxon of the shock, following the 
translational shock, the vibrational degrees behave almost classically. 
Hence in the kinetic theory developed here, all the degrees of freedom 
including the vibrational ones are treated classically. 


3«2. Boltzmann Equation - 

To describe the dynamic state of an axially elastic sphero- 
cylinder molecule, a set of 13 coordinates is required. 


The set consists of s 

x - the physical co-ordinates of the centre of mass of the 
molecule 

c - the translational velocity of the centre of mass of the molecule 
with respect to the laboratory co-ordinate system 
e - the unit vector along the axis of the molecule to fix the 
orientation 

to - the angvilar velocity of the molecule with respect to the 
laboratory co-orxiinate system. 

B - the maximum amplitude of vibration 

<!> - the phase angle of the simple harmonic oscillations of the two 

atoms constituting the molecule. 

Thus for a gas made up of such molecules, the distribution function is a 
function varying in the hyper-space q , t where g denotes the space 
defined by ( x, c, e, co,B,(!)), The probable number of molecules with 

y 

their phase coordinates between q and q +dq at any time t , is given 
by 


f(q;t) dq 
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fCa+o dt jt 4 -dt)dq' = £(Q;t)tia + {G-L} dq dt 


(3*2.5) 


where the phase volme dq’ at 0 + q dt, t + dt is reifiihd . to dq 
at (q , t) by 


dq ' 


% titj 


dq 


(3.2.4) 


where J is the Jacobian. G dq dt is the number of molecules reaching 

0 «r 

dq' during dt from outside dq j 03? the number of molecules 'gained 
by do' due to colliBlons; and L dqdt is the number of molecules in 
dq failing to reach dq’ during dt, or the number of molecules 'lost' 
due to collisions. Thus, we have (by expanding the left hand side 


and canoelling the phase volume term on either side ) 


3 f 

at 


+ 



G-L 


Iho following are the equations of motion of a free moleciae: 


(3.2*5) 


X = c 
■** 

c = F . the force per unit mass 

e = M X e 

^ w ««■ 

j*w“0 

Prom the last equation we have 



J = 



i— ) e e 

h ^ " 


( 3 . 2 . 6 ) 


where 
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or 


£Cx,c,e,a),B,<j)jt)dx dc de du dB d(}) 

where f is the distribution function describing the state of the 
gas and 


clx = dx dy dz ; dc a dCy do^ 

de = sin 3 dp de j du » dtoj do )2 dwj . 


( 3 . 2.1 ) 


Thus the number density of the molecules at a point x in the physical 
space and at a given instant t is given by the normalisation relation- 
ship; 

nCXjt) a / £(x,c,e,to,B,<j)jt) dx dc de dio dB d(j> (3*2. 2) 

Hie variation of f with time t, is now studied to obtain a 
governing equation for f » In the absence of collisional effects, 
the time derivatives q , of the 15 phase coordinates of a molecule 
are functions of the coordinates themselves. After an interval of 
time dt from the instant t, a free molecule initially having q as 
its co-ordinate will occupy a point q- + q dt in the hyper-space. Nov; 
consider a set of molecules in a phase volume dq at time t around 
a point q . At time t + dt the same molecules will be present in 

a phase volume dq' around the point q + d dt, if they suffer no 

'J ^ ^ . 

collisions during the interval dt. However due to collisional effects, 
a fraction of the total number of molecules in dq , fail to reach dq' 
v/hile some molecules which were not present in dq at time t, manage 
to reach dq' . The total number of particles in dq' , at time t+dt, 

is related to the number in dq at t by the following equation. . 
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so that 


J2-Ji 


(j) ^ _ — _ — (u • e) w >: e 

' Jj 

Substituting the above resiilts for in Eqn* (3.2.5) we get 

3f , . , r ^ \ ri ^ 

^+c •~+F •— + — { — ; — (tjj*ejca x e £} + — 

3t - 9x - 3c 3<;) J, - 


(3.2.7) 


{(0 X e£} 


^ Ik . G - L.« J . 

3B 3<|) 


( 3 . 2 . 8 ) 


Ihe quantities G and L are some functionals of the distribution 
function whose formulation is dealt with, in the following sections. 

3 * 3 . Collision Geometry and the Interaction Surface: 


Consider a collision between molecule 1 described by the 

coordinates q and molecule 2 vrfaose coordinates are q„ • Hie 
-1 ~2 

location vectors of the two molecules and differ negligibly 

at the instant of collision. Ibr the type of the geometry of the 

molecule assumed, there eod-sts a common tangent at the point of impact 

of the two collision partners and also a common normal to the surface 

of t\TO molecules at this point* Hhus a collision between two molecules 

1,2 is completely described by the coordinates ( o , e^, ^ 

for molecule 1 , ( t a- ) for molecule 2 and the common 

d *- 

normal h at the point of Impact, iin 'interaction surface' enclosing 
an ’excluded volume' is defined in the case of a bimry collision, as 
the locus traced by the centre of mass of molecule 2, as it moves in 
contact ?d.th molecule 1, fixing the orientation of both the moiecules, 

A point on the interaction siirface is separated by a vector R from 
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fCq+o dtjt+dt)dq’ = f(q;t)clq + {G-L} dq dt 


(3.2.3) 


where the phase volvune dq’ at q + q dt, t + dt is relpitqd, . to dq 

at ( q , t ) by 

*>• 

[q + q dt] 


dq' = 


( 3 . 2 . 4 ) 


where J is the Jacobian. G dq dt is the number of molecules reaching 

dq’ during dt ffom outside dq j or the number of molecules 'gained' 

** « 

by dq’ due to collisions? and 1 dq dt is the number of molecules in 
dq failing to reach dq’ during dt, or the number of molecules 'lost' 
due to collisions. Thus, we have (by expanding the left hand side 
and cancelling the phase volume term on either side) 


af 

^ = G - L 

3q 


{ 5 . 2 * 5 ) 


The following are the equations of motion of a free molecule » 


X = c 


c = F , 


the force per unit mass 


(3.2.6) 


e = oi X e 


(J*w) “J 


]?roni the last equation Yre have 




OD - J • J • (D 




1 1 


•) e e 


where 
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the eeatre of mass of molecule 1. 5f*om the geometry we have 

(Slg. 5.3.1 •) 


5 ' -Ij - l2 i - ^^2 


(30.1) 


go that 


n . dR = n , drj^ - B*^i2 " ^ 


(30.2) 


Hence n is saso normal to the interaction surface. Ibr l^ae case 
of a rigid spherocylinder it has already heen established by Ifeight 
(1965) that an element dS on the interaction, surface depends on 
61,62 , n for a given hind of molectsle and is related to an element 
of surface dn on a unit sphere by 


dS * K (ep & 2 , n) dn (300) 

Extending this for the axially elastic sphero— cylinder case. 

dS = KCej.e2.n,ri,B2>(l>j.<^2) *} (300) 

where 

1 * 1? ’^1 *^2 *^ 2 ) ~ ^ 1 ) 

7 2 

d (ei .n) (e.xn) + (be *^^2) ^(-©2*^) (©2^^) } 

+ {(bj+Bicos <t<p5C®i-5) + (^2*^2 (5.3.5) 

In Eqn. (5.3-I ) and r2 are functions of ( e^ , n , ) and 

(e, , n , B„, ) respectively, where e. and e, depent^ on f^r^' " 

rotating molecules. Thus the interaction surface S is variable with 

time. In obtaining the collision frequency we need the normal velocity 
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so that 


0 ) = - 




J, 


(w * e) w X 6 


(3.2,7) 


Substituting the above results for q in Eqn. (3.2.5) we get 

P 


9f 

+ C 

9t 


+ F • |i + i-* {^i^(u)*e)a) x e £} + - 
3x - 3c gw J, - 


'1 


3e 


{( 1 ) X ef} 


. M=g. l.« j. 

3B 34. 


( 3 . 2 , 8 ) 


The quantities G and L are some functionals of the distribution 
function whose formulation is dealt with, in the following sections, 

3 . 3 . Collision Geometry and the Interaction Surface ; 

Consider a collision between molecule 1 described by the 

coordinates q and molecule 2 whose coordinates are « Ohe 

-1 ~2 

location vectors of the two molecules and X 2 differ negligibly 
at the instant of collision. lor the type of the geometry of the 
molecule assmed, there exists a common tangent at the point of impact 
of the two collision partners and also a common normal to the surface 
of t-wo molecules at this point, Thus a collision between two molecules 
1,2 is completely described by the coordinates ( 0 . , e^, , 4^ ) 

for molecule 1 , ( Cgj Sg » Wg# B^, ) for molecule 2 and the common 

normal n at the point of impact. An 'interaction surface' enclosing 
an 'excluded volume' is defined in the case of a binary collision, as 
the locus traced by the centre of mass of molecule 2, as it moves in 
contact with molecule 1 , fixing the orientation of both the molecules, 

A point on the interaction surface is separated by a vector R fiom 
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of a point on. tho interaction surface. ]?roiii Eqns. (2»2«1 ) and 
we can write 

rf;: V 

g = n . — = V sin ^ ^2^ t3*3*6j 

~ dt 

where c is the normal component of the velocity of a point on the 

n,S 

interaction surface S, 

She interaction surface and the nonsal velocity of a point on 
the surface for a collision between an atom and a molecule may be 
obtained as a special case of the above results by putting 



= 0 and 1^2 ~ 


for the second collision partner vis. atom. 


3.4. Gollision Erequenc^ 


fo write down the Gain and the Loss teims, G and L in Bqn. 

(3,2.8) we need the number of collisions between type 1 and type 2 
molecules taking place during an interval t,t + dt in the region 
2; , X + dx of the physical space. fhe assumption of molecular chaos 
is made in finding an expression for collision frequency. Thus the two 


sets of molecules considered in the phase space regions q , q + dq 
and q } q + dq are assumed to be distributed at random without any 

.V 2 tv 2 2 

correlations between linear velocityj position, angular velocity^ 


orientation j maxitmjm amplitude of vibration and the phase angle of 
vibrations. In other words, all the pre-collision variables of the two 
sets of molecules, are uncorrelated. Ibr long sphero cylinder molecules, 
chattering collisions, in which the pre-collision variables are strongly 



the centre of mass of molecule 1. Itom the geometry we have 
(ELg. 50.1.) 
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5 “.ri - r2 5 dR = dri - dr2 


( 3 . 3 . 1 ) 


so that 


n . dR = n . rlr^^ - n.dr 2 = 0 (3*3.2) 

Hence n is also normal to the interaction surface. Ibr the case 
of a rigid epherooylinder it has already been established by Haight 
( 1965 ) that an element dS on the interaction svirface depends on 
®1»®2 * - given kind of molecule and ie related to an element 

of surface dn on a unit sphere by 


dS a K (ej, 62 f n) dn (3*5*3) 

Extending this for the axially elastic sphero-cyiinder case, we 

dS = KCe2,e2>n,Bj^,B2,^j>'Ji2) dn ( 3 * 3 * 4 ) 

where 

“ 4a + 2a{(b2+B2 cos tjij) 

2 2 
d(ej^.n)(ej^xn^‘‘ + (b2 + B2 cos <{>2) 6(“e2.n) (©2^^) ^ 

+ iCbj+B2cos <l'j)d(®i-2) (b2"''®2 't'2)'^ ('■®2*5) J (2l’*S2’^2)^ (3*3*5) 

In Eqn, (3,3.1 ) Tj and T 2 are functions of ( gj , n , ) and 

(e^ f n . B„. (to ) respectively, where e, and e depend, on 

•r/ ' ^ ^ ^ 

rotating molecules. iChus the interaction surface S is variable with 
time. In obtaiaing the collision frequency we need the normal velocity 
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dependentj are possilsle* The effects of such chatteriag collisions 
jna^ be reduced by choosing a small b/a. 

A nioleevile of type 2 has a velocity c. ) relative to the 

-HA -»1 

caatre of mass of molecule of type 1 . E\xrther the normal component 
of the velocity of a polntj on the interaction surface with respect to 
the centre of mass of molecule 1 is non-aero and is given by Eqn, (3.5»6)* 
Hence the velocity of the centre of mass of molecule 2 in a direction 
normal to the interaction surface and relative to the surface is given 

by 

where g is given by Eqn. (2.3.6). 

A collision between molecule 2 and molecule 1 (or collision of 
type I- 2 ) would take place when the centre of mass of molecule 2 
occupies a point on the interaction surface defined. Thus, the 
molecules of type 2 whose centres of mass are lying inside a cylinder of 
base area dS on the inteiaction surface S and length - (g n) dt 
noimal to the surface (see Eig. 5»5»1 ) would encounter a molecule of 
type 1 during a time interval dt . The probable number of type 2 
molecules lying inside this cylindrical volume (- g . n dt ds) in 

•V 

physical space and having co-ordinates between a and q + d q is 

-2 ,2 Z ?. 

^^2 

This also gives the ntraber of type 1-2 collisions taking place dijring 
the time interval dt for one molecule of type 1 present in the phase 
space q , q + d q , 


Thus the total number of 1-2 collisions taking 
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of a point on the interaction surface. Prom Eqns. (2.2.1 ) and (3*3»1 ) 
we can write 


Cyj g = n . — = V sin *^2^ 

, ~ dt “ ~ ~ 


(3.3.6) 


where c is the normal component of the velocity of a point on the 
interaction surface S. 


Ihe interaction surface and the nonnal velocity of a point on 
the surface for a collision between an atom and a molecule may be 
obtained as a special case of the above results by putting 


b 2 = 0 , = n , U 2 ” ® > "^2 ~ ^*2 “ 

for the second collision partner viz. atom. 

3.4. Oollislon Prequenoy s 

To write down the Gain and the Loss terms, G and 1 in Bqn. 

(3.2.8) we need the number of collisions between typo 1 and type 2 
molecules taking place during an interval t,t + dt in the region 
x , X + to of the physical space. The assumption of molecular chaos 
is made in finding an expression for collision frequency. Thus the two 
sets of molecules considered in the phaso space regions q » q. + ^<1 
and q^ , *^*^2 assumed to be distributed at random without any 

correlations between linear velocity, position, angular velocity, 
orientation, maximum amplitude of vibration and the phase angle of 
vibrations. In other words, all the pre-collision, variables of the tvro 
sets of molecules, are uncorrelated. Por long spherooylinder molecules, 
chattering collieions, in which the pre-collision variables are strongly 
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place ia the phase space q ^ q + d q per urdt time is given by 

•"1 1 ^ i 


{ / f 


(3.4.2) 


2Ihe nutnber of collisions taking place in a unit volume and unit time 
known as collision frequency can be widtten as 


* / fj^f2 dc2<i«2^®2 ^^2 ^^2 (3t4#3) 

Jbr the case of collision between an atom and ,a molecxile the collision 
frequency can be y«ritten from Eqn. {3»4*3) as follows* SCaking the 
particle 2 to be an atcrm* we hare 0, 0, 0)2 » 0 and ?2 ** 5 • 

Jttrther the area dS^ on the in teraotioa surface S-j# in this ease is 
given by putting 0 and ©2 ‘ ^ ^ Bqns. (3#3#4) and (3«3*5). 

Ihus the frequency of moleciil e-atom collision is given by 

- / f jF2(gj.n)dSx dC2 dcj^ de^ dwj dBj d<{ij (3, 4.. 4) 

gl.n<0 

where Pg is atomic distribution function and 

% ~ * ^Si !^i ^ ~ 2 i '^’ 1 ^ ( 3 * 4 * 5 ) 

3.5- Collieiofl Integral j 

A molecule of type 1 loses its identity when it undergoes a 
collision* Ihue the number of laolecules of typo 1 lost from the 
phase 'volume q. , q* + dq. during a time interval dt due to collisions 
Is given by 

,{ f • (^* 5 * 1 ) 
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dependent, are possible. The effects of such chattering collisions 
may be reduced by choosing a small b/a. 

A molecule of type 2 has a velocity (c - c. ) relative to the 

•Tii -ri 

centre of mass of molecule of "type 1. Thrther the normal component 
of the velocity of a point, on the Interaction surface v/ith respect to 
the centre of mass of molecule 1 is non-zero and is given by Eqn. (5*3«6). 
Hence the velocity of the centre of mass of molecule 2 in a direction 
normal to the interaction surface and relative to the surface is given 
by 

” (S • (3«4.1 ) 

where g is given by Eqn. (2 .3 •6). 

A collision between molecule 2 and molecule 1 (or collision of 
type 1-2) would take place v/hen the centre of mass of molecule 2 
occupies a point on the interaction surface defined. Eius, the 
molecules of type 2 whose centres of mass are lying Inside a cylinder of 
base area dS on the interaction surface S and length - (g n) dt 
normal to the surface (see Big. 3.3*1 ) would encounter a molecule of 
type 1 during a time interval dt . The probable number of type 2 . 
molecules lying inside this cylindrical volume (- g . n dt ds) in 
physical space and having co-ordinates between q and q + d q is 

t'-'t tiS) dc2 de2 ^^2 ^^2 ^*>*2 

This also gives the number of type 1-2 collisions taking place during 

the time interval dt for one molecule of type 1 present in the phase 

space q , q. + d q Thus the total number of 1-2 collisions taking 

V 1 i ..1 
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where Y(jO - for Y ^ ^ (3*5«2) 

= y for y T n 

The g-uantity inside the hmcket corresponds to the loss term 
If in the Egn, (3*2*8). ly cbcosing properly the pre-collision 
variahleB ( ^ ^ ^ ^ ^ ^ ^ X ^ ^^d ( x ^ ^ ^ ^ 

for the two sets of molecules, it is possible to obtain molecules with 
post collision -variables as ( x , e-j » » 4 > 3 _) and 

( X , c , {j , e„, d^). In such a collision a molecule of type 1 

is generated. Using the expression for the collision frequency 
Eqn. (3*4-.5)5 can write down the total number of 1’* 2* collisions 
occurring in time dt as 




(I f 'f^yC-f ' .n)dS dc2 ^62 dUl 04^} dq^ dt 


f . ’ 


( 3 . 5 . 3 } 


^tom Iiiouyille*s theorem of statistical mechanics we tare for the 
conservation of phase volume* 

t , t , f > ' - * . ' 

' rr* 


iCjCc2 do}^ dl ^ dC>j c-^^2 ” 

'‘ 5i ^52 ^'^1 ‘-^2 


(3.5*4) 


It can also be show that 


(g.n) = -Cg’.n) 


( 5 . 5 * 5 ) 


Using these results, the number of 1 '-* 2' collisions may be written as 
^'“S*^) dS dc2 ^^2 de2 ^(^2 (3.5*6) 

This also gLres the number of moleculce of type 1, generated in time 
dt, in space q^, q^ + d q^ * Hence the quantity in the brackets in 
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place itt the phase space q , q + d q per unit time is given by 

••1 -^1 V 1 


{ / dC2 dco dw2 dB2 d(>2} dqj 


(3.4.2) 


The number of collisions taking place in a unit volume and unit time 
knovtti as collision fiequengy can be written as 

“ / fjf2(g.n)dS dC2dm2de2 <1^2 d4i2 dcj du)]^ de^ dEj^ d<)j^ . (3.4.3) 

g.n<0 ~ - 

^ *4 

Ibr the case of collision between an atom and, a molecule the collision 
frequency can be written from Eqn. (3*4»3) as follows: (Baking the 
particle 2 to be an atomj we have BgS* 0, <;) 2 »i 0| ‘;)2 ^ ?2 * 5 ’ 

further the area dS^ on the interaction surface in this case is 
given by putting b^s 0 and ^2 ' 0 ” ^ ^ Bqns. (3.3.4 ) and (3*3*5). 

(Thus the frequency of molecule-atom collision is given by 

- / dc2 dcj dej dwi dB^ 6 ^^ 

gl.n<0 

where Pg atomic distribution function and 

li “ 32 “ (3i * <£1^ h “ 3i 

3.5. Collision Integral : 

A moleovile of type 1 loses its Identity when it undergoes a 
collision, (Blue the number of molecules of type 1 lost from the 
phase "volurae q^^ , q^ + dq^ during a time interval dt due to collisions 
is given by 

i /' fjf2Y(-g.n) dS dc2 dt02 de2 dBo d(j)2} dq^ dt, (3.5*1 ) 


(3.4.4) 


(3*4.5) 
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S^,(5*5*6) correspon<3s to the gaiii teira G in Eqn. ( 5 . 1.7), Ohus 
we have the expression for the gain and logs terms for molecules of 

lying in q^, q^ + dq^ . Krorn this we can write the collision 

integral J for the siolecules in q j q + dq by dropping the 

** P* tv 

subscript 1 and changing subscript 2 to 1 as 


J = / f' 


l("S-n)d8 citj^ (3»5.7) 


where 


^rj_ = dej_ ( 3 . 5 . 8 ) 

This completes the formulation of the Boltzmann equation for a gas 
with axially elastic sphero -cylinder molecules. 

3.6. tSoment Equations : 

'2ie corresponding Maxwell transfer equation for this case may 
be obtained as usual by multiplying both sides of tho Boltzmann 
equation (3»2.8) by a molecular property i|j( x ^ c j e , w , ^ ) 

and integrating with respect to c , m « e , B suad 0 . This 

t¥ JV 

gives us the f:>llowing equation 



c 



3f 3 
3c 


{ — - — (m.e) (a)>ro)f) 

Jj — " ^ 


3e 


Cwxe)f + 


21 " 


'3 ^J'(x,C,©,w,'3,<f>)dT 


" j J di 


( 3 * 5.1 ) 
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v;i>.ero Y(y) " ^ (3*5*2) 

= y for y 7 0 

^Che quantity inside the bracket corresponds to the loss term 
1/ in the Eqn* (3*2*8), choosing properly the pre-collision 

variables ( x , c^', e^ ’ V ’ ^ ^ ' -2 ' 22 '^2 ’ ‘^’2 ^ 

for the two sets of molecules, it is possible to obtain molecules with 
post collision variables as ( x , , 03 ^^, ) and 

( X , c , to_> e„, B_, ^ collision a molecule of type 1 

is generated. Using the expression for the collision frequency 
Eqn. (3*4*3 )» we can write down the total number of 1 2* collisions 
occurring in time dt as 


{/ fjf2V(-p».n)dS dc2 clw^ d<i2 <iql dt 


(3.5*3) 


3h?om Liouvllle‘s theorem of statistical mechanics Vi;e have for the 
conservation of phase volume. 



dci ^32 dm^ dm 2 ^<>1 ‘^<'2 


It can also be shown that 


(3*5.4) 


(g.n) = -(g’*") 


(5*5.5) 


Using these results, the number of 1 •- 2' collisions may be written as 


{/ YC-g.n)dS dc 2 dw 2 de 2 (102 dij) 2 } doj^ c?:t (3.5*6) 

This also gives the number of molecules of type 1, generated in time 
dt, in space q|, q^ t d . Hence the quantity in the brackets in 
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Using the folio vdng notation for the moment of f 


/ dc uej; cia d.3 d<^ = rL(x,t) ^(x^t) 


( 3 . 6 . 2 ) 


where n( x , t) is the nxmiber density of molecules at a point z 
and time t^ the Maxwell transfer equatioHj for the case where j|, is 
a function of only o ^ ® ^^.ay be simplified to the form 


— B * nc xj.* \ ^ 1 


Bn il> + 


- n {F . + — r- 

“ Dc d-l 


(ra^.e^) w X e. 


B y(-F.i^) K(e,e^,n, 

Bo) cl- 1-^ 

dn dcj de^ dLj 


( 5 . 6 . 3 ) 


in v/hich Eqn. (3.5.4) ie used for ds. !Ihe right hand side of this 
equation, which is the change ( in the property ,]j due to 

collisions, may be vsrltten in different forms by effecting some changes 
in the variables of integration as done by Bright (1965)* Kius one 
of the useful forms is 


Uie prined quantities ( 0 ^' , mj,, and (c*, p, B*, <!>'} 

are the pre-collision collision variables of two colliding moleoriLes 
giving rise to molecules with post collision variables ( c^, , (j^p 

and ( G , m , B , ^ }. In the present case, a sunamational invariant xji 
Is defined, which satisfies the condition 

4- B ( 5 . 6 . 5 ) 

!2hus, in the Maxwell transfer equation, if is a summational 
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Egi!i,(3»5.6) coiresponds to the gain terra G la Eqn. (3.1*7)# ^I!hvls 
we have the expression for the gain and loss terras for molecules of 
type 1 f lying in + dq^ . Erom this we can write the collision 

integral J for the molecules in q » q + dq by dropping the 
subscript 1 and changing subscript 2 to 1 as 


J = / f’ f* vCh.n) dS + f Y{-}j,n)d£ dxj (3»5.7) 


where 


'’■^1 “ ‘-’Sl ^’?i '’•^1 (3*5.8) 

Ihis completes the formulation of the Boltzmann equation for a gas 
with axially elastic sphere -cylinder molecules. 


3.6. Ivfoment Equations ; 

[Ehe corresponding Maxwell transfer equation for this case may 
be obtained as usual by multiplying both sides of the Boltzmann 
equation (3*2.8) by a molecular p3X)perty i})(x^o,e,w,B, ^ ) 
and Integrating with respect to c , m } e , ■ g ' and 6 • Siis 
gives us the follov^lng equation 


/ ( 


Df 

9t 


+ c 


M. 

ax 


+ p 


M 

9c 9(u 


J1-J2 

I — - — (w.e) (wxc)f} 



(raxe)f + 


afif ^ 

as 3^ 


■] \p(x,c,e,a),3,(l>)dT 


= / J >fi cH 


(3.6.1) 
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invariant:, the collision integral on the right hand side, vanishes. 
A oarefuL examination of the collision dynaaiics, enables us to find 
the following summational invariants 

1 




n c 

1 ^ ? It 7 



any function of e 


any function of 
o.ny function of 






0 


(5-6.6 ) 


% choosing the first three of these, the familiar conservation 
equations for mass, linear momenttaa and energy my he vsritten as 
follows: 


|a. » . nu = 0 

St 3x 


u 

at * 


;)x * L ^ ^ “ 

ite (91, „ . ii, _ !i, 

2 St " 3x 


(3.6.7) 


3x 


where u , 1, q and p are defined as the Sallowing moments of the 
distribution function f ; 


u = 


n 


/ C f dT 


and 


7 

T nk 


f {j K C * 2 ‘k 


j sCr;)“) f dT 


( 3 - 6 . 8 ) 


m c c f dx 


n =: / c r c' + i Jj \ * J ?(2n'^}Ux 
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Using the following notation for the moment of f 


/ dc do; de dB = nCx,t) 


(5.6.2) 


where n( x , t) is the number density of molecules at a point x 
and time t, the Maxwell transfer equation, for the case where ^ is 
a function of only c , 01 and B may be simplified to the form 



dn dC]^ dmj^ dej^ dBj d(|>j^ (3.6.3) 

in which Bqn. (3.3*4) is used for ds. dhe right hand side of this 
equation, which is the change ( AiIj )^^^ in the property ijj due to 
collisions, may bo written in different forms by effecting some changes 
in the variables of integration as done by Haight (1965). ®ius one 
of the useful forms is 


" I" ’*'1 ■ ^^1 ‘^■^1 (3.6.4) 

The primed quantities ( c^' , , <1)^ and (o’, m’, B', <)' ) 

are the pre~Gollision collision variables of two colliding molecules 
giving rise to molecules with post collision variables ( ) 

and ( c , CO , B , ({1 ). In the present case, a summational invariant i}/ 
is defined, which satisfies the condition 

lii^ + ijU - ^ ;|J a 0 (3.645) 

' ' ■ \ ' ■ 

Thus, in the Maxwell transfer equation, if ijj is a summational 
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Out of these moments, u and T refer, respectively, to the 
macroscopic velocity and temperaturej p and q are related to 
the pressure tensor P and the heat flux vector Q by the following 
relationships : 

- 2 

q = Q +p (3.6,9) 

and 

t! = P + p U U 

It should be noted that, though the angular momentum, is conserved 
during a collision, it is not a sujmaational invariant due to the 
tvTo distinct location vectors chosen for the tvro collisional pariacers. 

3,7* H-!Iheor&m and the Bqtrilibriuin Distribution ; 

In this section, an equilibrium solution is obtained for the 
Boltsmann Sqmtion (3.2,8), through an H-fheorem. ®ie nimber of 
molecules with dynamic variables q. In a phase volume q. , q. + dq. 
is given by f(q., , t) dq. . If the state of a g^-s changes in a time 
interval dt, this number also changes to (f^+ Af^ ) dq.j . The 
fractional change in the number of molecules accompan5dng a change of 
the state of the gas is ( Af^/f.j) dq.j or ( A Inf .j ) dq^ . Such a change 
takes place for molecules with different dynamic variables. Thus, it 
is seen that Inf is an important molecular properly associated with 
the state of a gas. Choosing » Inf in the Hfexwell transfer 
Equation (3*6 *3 we get 

_ £f 

. Cnc^lHF)= YC-g.n)K dn dr, dt (3.7.1) 
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invariant, the (xillision integral on the right hand side, vanishes. 

A careful examination of the collision dynamics, enables us to find 
the following summational invariants 

1 

m c 

- c + ■:7 J-i tOj. + “ 5(23) 

= any function of e ( 3 . 6 . 6 ) 


any function of 
any function of 





Ily choosing the first three of these, the familiar conservation 
equations for mass, linear momentum and energy my be written as 
follows s 
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( 3 . 6 . 7 ) 
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+ r) : 


ax 9x 

tv ^ 

where u , T, q and p aro defined as the fbllowdng moments of the 
distribution function f : 


u = — / c £ dx 

- T1 ~ 


and 


/ £~ r.’ '-i " J 5 (2B)^} f dr 


7 

T nk 

rt -■ f m c c f dx 


( 3 . 6 . 8 ) 
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sC2B)'^}£dx 
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Porther, it can te shown that for the case of an axially 
elastic sphere cylinder particle the right hand side of Eq.n. (3.7*1 ) 
is always negative or at laost ecLuai to zero. In the case of rigid 
spheroeylinder molecules a similar result was used by Pfoight (1965). 
Thus for a uniform case^ we can write the above equation as 

ff 

/ In C“7^) ffj yC*|.P-) K dn dr dt^ (3*7*2) 

f 

<0 

where 


H = n TnT « / £ 2n£ dr (3.7.3) 

Eqh. (3.7*2) is the statement of the H-fhoorem, following the 

argueneat given by Chapman and Cowling (l970)^ it can be concluded 

that H decreases monotoraically to a constant value corresponding to 

the case s 0* f^irthe^, it is also knewn that this corresponds 
dt 

to msKimiM entiXJpy or statistical equilibrium* ffihus for equilibrium, 
we have from Eqn. (3.7*2) 

Inf^' + Inf' = inf^ + Inf (3.7.4) 

Composing this v/lth Bqzi* (3*6*5) it can oasiljr be seen that Inf 
is a Busairiational invariant or a linear combixiation of the summational 
invariants (5.6*6)* fhus we can Inf in the form 

g^2B)^}+a4F^Ca)jj)+a^ G3^C<{i) (3*7*5) 

where » ^*3 > and are the quantities which may be 

related to physical par^erers like number density, temporature 
and velocity of the ^s. Using the normalisation relationship^ 
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Out of these momenta, u and T refer, respectively, to the 
macroscopic velocity and temperature; p and q are related to 

mr 

the pressiAr© tensor P and the heat flux vector Q by the following 
relationships; 

2 

q = Q +p (3.6.9) 

•"f M A 

and 

p S= P + p U U 

It should be noted that, though the angular momentum is conserved 
during a collision, it is not a sumraational invariant due to the 
two distinct location vectors chosen for the two collisional partners. 


3.7. H-Theorem and the Equilibrium Distribution : 


In this section, an equilibrium solution is obtained for the 
Boltzmann Equation (3.2,8), through an H-Theorem. The number of 
molecules with dynamic variables q. in a phase volume q^ , q, + dq. 
is given by ^(q- , t) dq. , If the state of a gas changes in a time 

IV ' «v I * 

interval dt, this number also changes to (f^+ Af^ ) dq^ . The 
fractional change in the number of molecules accompanying a change of 
the stats of the gas is ( A or ( A Inf^ ) dq^ ♦ Such a change 

takes place for Etolecules with different dynamic variables. Thus, it 
is seen that Inf is an important molecular property associated with 
the state of a gas. Choosing » Inf in the Maxwell transfer 
Eq\:iation (3«6.3), we get 


3n Inf ^ ^ 
Dt ax 




Cnc ln?)= i / ln(--^]£f, Y(-g.n)K dn dt, dt 
f'f' " " " ■ 


(3.7.1) 



(3.6 *2) £^d the definitions (3*6.8) of velocity and teu^eratnre 
the Eqn. (3*7.3) can he witten in a standard form as 
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B , TP . 5/2, 

4ir hnkT 4irkT 


.U«* 

■ w 


> fCuIj pc Cd) (3.7.6) 


where 


/ PCwi[) = 1 


(3*7.7) 


and &(<!>) is the distribution function for the phase an^e f 
of the vibrational state of the molecules. In the absence of any 
other infomation, the distribution G-( } for will be assumed to 

be unilorm. Ihus we shall tales 


G(d) - 


Ztt 


(3.7.8) 

^\irth ST in writing Eqn . (3.7*5) the fo 11 owin g approximation 
is employsdt 



^max^ ” Y" * 



2 


(3.7.9) 


This is justifiable as the value of for is 115000<’K 

(see Ihble 2.1.1). fhus we hare the J^wellian distribution for 
the case of a vibmtionally relaxing diatomic gas. In the next 
chapter two solutions for the plane shodc structure problem is 
discussed for a vibrationaily relaxing di atomic gas. 
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Further, it can he shown that for the case of an axially 
elastic sphere cylinder particle the right hand side of Eqn, (5«7<1 ) 
is always negative or at most equal to zero. In the case of rigid 
spherocylinder molecules a similar result was used by Haight (l965). 
Thus for a uniform case, we can write the above equation as 

ff 

^ = i / In (~-^) f£^ vC-g.n) K dn dr dtj (3.7.2) 

f V. . 

<0 

where 


H = n Inf =5 / f Inf dr (3*7.3) 

Eqn. (3.7*2) is the statement of the H~!Rieor6m. Fbll owing the 

orgueaent given by Chapman and Cowling (1970), it can be concluded 

that H decreases monotoiiically to a constant value corresponding to 
dH 

the case = 0. Further, it is also known that this corresponds 
to maximum entropy or statistical equilibrium. Thus for equilibrium, 
we have from Eqn. (3«7.2) 

Inf^' + Inf ' = Inf^ + Inf (3.7.4) 

Comparing this v/ith Eqn, (3.6,5) it can easily be seen that Inf 
is a summational invariant or a linear combination of the suimBational 
Invariants (5.6.6). Thus we can write inf in the form 

lnf=a^+a2.c+a3{^c2+ijjo,f+i s(2B)2}+a4Fj^CmiP G^(.^) (3.7.3) 

where , ci^ , and are the quantities which may be 

related ‘fco physical parameters like number density, temperature 
and velocity of the ^s. Using the normalisation relationship^ 



CHAPITER - 4 


KJEMULAHOH OP SHOCK SQSUGPCRE HtOBREM 


4.1, latrodiiction : 

Due to the complexity of the collision terms in the Bcltsmann 
equation (especially the gain term) it has not been possible (till now) 
to solve the equation even for problems vdth simple boundary conditions 
0?o overcome this difficulty some models for the Boltzmann equation were 
proposed by making some reasonable assimiptions about the collision proc 
by Bhatuagar, Gross and Krook (1954), Gross and Jackson (1959) for the 
monatomic gas case and Brau (1969), Itorse (l964), Hanson and Jforse (196'" 
for the diatomic gas case. 5br the case of uKtnatomic gases, the model 
eqmtion was solved numerically by using an iteration scheme by 
Ijiepraann et al (1962), 

Another approach to solve a gasdynamio problem from a molecular 
point of view, consists in choosing a suitable ansata for the xmknown 
distribution function leaving a few free parameters in it. These free 
parameters are determned by using as many moments of the Boltzmann 
equation as the number of free parameters in the ansatz, The solution 
thus obtained isuociy an approximate one and depends stron^y on the fo: 
of the ansatz used and the choice of the molecifl,ar properties used in 
the moment equations, 

4.2 Ansatz for the jPistributiQn Punction : 

Mott-anith ( 1951 ) proposed a bimodal ansatz for the distribution 
fmction inside a plane shook wave in a monatomic gae considering the 
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( 5 . 6 . 2 ) and the definitions (3«6.8) of velocity and temperature 
the Eqn. ( 3 . 745 ) can be written in a standard form as 



n / m ^ 3/2 
4Tr ^2iTkT‘^ 


( tl 0 C^) 

2irkT itkT 


1/2 


exp{-, 


nic^+JTtoJ+s(2E)^ 


2kT 


} B* (tOj j3G (<{i) (3*7 46 ) 


where 

oc> 

I F(wi()da).,=l (3.7.7) 

4-,CO » ^ 

and g( iJ* ) is the distribution function for the phase angle <j) f 
of the vibrational state of the molecules. In the absence of any 
other information, the distribution g( <{> ) for <!> will be assumed to 
be uniform. Thus we shall take 

G(<!.) = ( 3 . 7 . 8 ) 

Further in writing Eqn. (3.7.5) the following approximation 
is employed: 

er£(^ W " T- ^ 

This is justifiable as the value of T^^ for Eg is 113000®K 

(see Table 2.1.1 ). Thus we have the Maxwellian distribution for . 
the case of a vibrationally relaxing diatomic gas. In the next 
chapter two solutions for the plane shock structure problem is 
discussed for a vibrationally relaxing diatomic gas. 
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physics of the problem* Ibllowing the tsro step analysis for shook 
waves in a diatomio gas given by 5!albot and Scala (l96l) within the 
frame work of fevier-Stokes theory, Ifeiglit (1965) extended the Mott- 
Smith bimodal theory to a diatomic case, Ba proposed a trimodal 
distrihution function for the shock problem. Both fifott-Smith and 
Haight used moment methods to solvo for the unknown or the free 
parameters in the ansats. Ihe aasata proposed by Haight ia capable of 
predicting only the effects of rotational excitation on the structure 
of a shock wave. In the present analysis vibrational effects are also 
included in the formulation* Ihrther in Chapter 5 the dissociation 
effects are also considered for the shock structure problem. 

A s»re general ansata was proposed by Weitssch (l96l) for "the 
distribution function in the case of strong departure from thermodynami 
eguilibrlua as an integral superposition of Maxwellians with suitable 
wei^ting functions, She integral representation may be replaced by a 
quadrature formula as given by ISarasimha and heshpande (l968) for the 
distribution function as 

N 

f. I (4.2.1) 

a»l 

where is a Maxwellian built up using ( n^ ) and is the 

product of the weighting function chosen ia the integral representation 
and the weights in the appropriate quadrature formula. It can easily 
be shown that the ifett-SiHith biaiodal aneate is a partioular case of the 
geneiml form (4*2,1 ). Iktendiag a similar representation to the 
diatomic case, we have PMghts^ trlmodal ansats as a partioular case of 
Eqn, ( 4 * 2 . 1 ), 3he choice of the parameters used in building up the 
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4 


lORMULATION OF SHOCK SmUCTOE PROBLEM 


4,1, latroduotlon i 

Due to the oomplexity of the oollision terms in the Boltzmann 
equation (especially the gain tena) it has not been possible (till now) 
to solve the eqmtion even for problems with simple boundary oonditions. 

To overcome this difficulty some models for the Boltzanana eqvation were 
proposed by making some reasonable assumptions about the collision process 
by Bhatnagar, Cross and Krook (1954), Gross and Jackson (l959) for the 
monatomic gas case and Erau (1969), Morse (1964), Hanson and Morse (1967) 
for the diatomic gas case. For the case of monatomic gases, ihe model 
equation was solved numerically by using an iteration scheme by 
liepmann et al (1962), 

Another approach to solve a gasdynamic problem from a molecular 
point of view, consists in choosing a suitable ansatz for the unknown 
distribution function leaving a few free parameters in it. These free 
parameters are determined by using as many moments of the Boltzmann 
equation as the number of free parameters in the ansatz. The solution 
thus obtained isuonly an approximate one and depends stron^y on the form 
of the ansatz used and the choice of the molecular properties used in 
the moment equations, 

4*2 Ansatz for the Distribution Function : 

Mbtt-Smith (1951) proposed a bimodal ansatz for the distribution 
function inside a plane shock wave In a monatomic gas considering the 
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MEDcwellians is arbitrary in the general fora of the distribution. 

However j if the choice ie Hade on. the basis of the physics of the 
probleia at hand, one can expect a better representation of f. 

Purther the representation of f may be hoaproved by choosing as large 
an II as possible for a given problem, fhe value of the number H is 
limited by the magnitude of computational work involved and the ability 
to guess the parameters, to build up the Maxwellians, for a given 
physical situation. 

4.5 Quadrimodal Hypothesis : 

Scala and Talbot (1963) obtained the structure of a shock wave 
in a vibrationally and rotationally relaxing diatomic gas from a 
continuum approach. For this they extended their earlier analysis 
(l96l) by defining a third temperature for the vibiatlonal degrees 
of freedom. However for the present analysis the txpeiimentally 
established fact ty »■ is used to visualise a three step 

shock. In the initial narrow region (translational shock) only the 
translational degrees of freedom are assumed to be excited fully and 
in the central aone (rotational shock) the rotational degrees of friiedom 
are expected to attain- equilibrium. In both of these regions the 
vibrational excitation is taken to be negligible. Finally, in the 
broad tail of the shock (vibrational shock) the vibrational excitations 
occur completely. With such a picture for a shock wave, it is possible 
to choose four sets of parameters to build up four Maxwellians, (in 
an actual shook these three distinct regions do not exist). These 
Maxwellians are used in (4*2,1 ) to arrive at a quadrimodal distribution 
function for the shock problm. 
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physics of the problem. Bbllowing the two step analysis for shock 

waves in a diatomio gas given by Talbot and Scala (l96l ) within the 

frame work of Navier-Stokes theory, Ifelght (1965) extended the Mott- 

Sznith bimodal theory to a diatomio case. ife proposed a trimodal 

distribution function for the shock problem. Both Mott~Sinith and 

Haight used moment methods to solve for the unknown or the free 

parameters in the ansatz* The ansatz proposed by Haight is capable of 

predicting only the effects of rotational excitation on the structure 

% 

of a shock wave. In the present analysis vibrational effects are also 
included In the formulation. ]?urther in Chapter 5 the dissociation 
effects are also considered for the shock structure problem. 


A more general ansatz was proposed by Weitzsch (l96l) for the 
distribution function in the case of strong departure from thermodynamic 
equilibrium as em integral superposition of Maxwellians with suitable 
wei^ting functions. The integral representation may be replaced by a 
quadrature formula as given by Narasimha and Deshpande (l968) for the 
distribution function as 


£ * 


N 


0=1 


'"'a 


(4.2.1 ) 


v/here P • is a Jdhxwellian built up using ( u^,, Tg, n^ ) and is the 
product of the wei^ting function chosen in the integral representation 
and the weights in the appropriate quadrature formula. It can easily 
be shown that the JMbtt-Smith bimodal ansatz is a particular case of the 
gen escal form (4* 2,1 ). Extending a similar representation to the 
diatomic case, we have Efeights’^ trimodal ansatz as a particular case of 
Eqa, (4.2.1). The choice of the parameters used in building up the 
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* 


(Uhe foiir Maxwelliaas thus chosen correspond, to the 

following four states of the gas: 

1 , uaifom free stream equilibrium condition . 

2* intermediate state existing at the end of the translational shock. 

5, intermediate state existing at the end of the rotational shock, 

4. uniform down stream equilibriiaa conditions* 

fhe parameters used in building up these four Maxwellians are obtained 
in the next section. 


4.4. Parameters in the MaxwellianB ; 


Ifexwellian f.j- Bor a given upstream Ifech nimiber M.j^ temperature 
and pressure p.j f the parameters for the Maxwellian f.j may be written 
down easily. As the rotational degrees of freedom are excited 
completely, at the upstream temperatures conflldered 29?^*, 400®, 500’’ K), 

the value of rotational temperature (f ) for the upstream gas is the eanio 
as that of liie translational temperature. However, the vibrational 
excitations of molecules at such low (upstream) touperatures are 
ne^egibie and hence the vibrational temperature , may be taken to be 
sero. 

Huts we have 


H : T » 0 

b b ^ b 


(4.4.1) 


As a consequence of zero vibaratienal temperature, the distribution 
function corresponding to vibrational degrees of freedom, in the 
Maxwellian (See Bqn. 3.7.6) becomes a Dirac delta function 
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Llaxvvellians is arbitrary in the general foira of the distribution. 

However, if the choice is made on the basis of the physics of the 
problem at hand, one can expect a better representation of f. 
lUrther the representation of f may be improved by choosing as large 
an H as possible for a given problem, Jiae value of the number N is 
limited by the magnitude of computational work involved and the ability 
to guess the parameters, to build up the Maxwellians, for a given 
physical situation. 

4.3 Quadrimodal Hypothesis ; 

Scala and Talbot (1963) obtained the structure of a shock wave 
in a vibrationally and rotationally relaxing diatomic gas from a 
continuum approach. For this they extended their earlier analysis 
(1961 ) by defining a third temperature for the vibrational degrees 
of freedom. However for the present analysis the experimentally 

established fact t„ » is used to visualise a three step 

V r t 

shock. In the initial narrow region (translational shock) only the 
translational degrees of freedom are assumed to be excited fully and 
in the central zone (rotational shock) the oxitational degrees of freedom 
are expected to attain equilibrium. In both of these regions the 
vibrational excitation is taken to be negligible. Idnally, in the 
broad tail of the shock (vibrational shock) the vibrational excitations 
occin? completely. With such a picture for a shock wave, it is possible 
to choose four, sets of parameters to bXiild up four Maxwellians. (in 
an actual shook these three distinct regions do not exist). These 
Maxwellians are used in (4.2,1 ) to arrive at a quadrimodal distribution 
function for the shock pro bl on. 
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i.e. lim [{_!?■ . ,1^''^ exp (- iSSll)] = 5(0) (4.4.2) 

Tv, ° Tvi ' 2k 

1 

further in absence of vibrational excitation, the distribution 
G-((f) ) for the phase angle (f* is also dropped. fhus we can v/rite 
f ^ f finally as 


£,=-r-C 


m 


1 47r^2TTkT. 




2TTkT. 


-3©xp {- 


mCc - i + Jjwi 

2kT, 


} 6CB) ( 4 . 4 . 3 ) 


Maxwellian f - [Hie translational temperature (T. ), the velocity 
2 tg 

(u^) and the number density (n^) of the gas at the end of the first 
aone may be fixed by using the Ranhine-Hugoniot condition for 
monatonic case (y = 5/3). However, the rotational and vibrational 
temperatures (T ,T )ofthe gas are frozen at v^ues corresponding 
to their free stream values. 


i.e. 


where 


li - 

T, " 80 
1 

^2 


(1 


25/7 

TtT 2 

1 


•) 


r'2 

(»! 



25/7. 
2 ^ 




Mr=«f/(7/5|Ti) 


fhus we have the Maxwellian f^ as 


(4.4.4) 


(4.4.5) 




f - ^ 1 ^ 

2 Air 2ir]^'r^ 



■)exp { 


mCc-i U2)^ 

■"ITtJ" 


2 

— ) 5(B) ?(«.,,) 


(4.4.6) 
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Hie foiir Maxvvelliatis = 1^4> thus chosen correspond to the 

following four states of the gas; 

1 , uniform free stream equilibrium condition , 

2. intermediate state existing at the end of the translational shock. 

3, intermediate state existing at the end of the lotational shock. 

4. uniform down stream equilibrium conditions. 

The parameters used in building up these four Maxwellians are obtained 
in the next section. 


4.4. Parameters in the Maxwelliaqs t 


Bfaxwellian f ^ - BVir a given upstream Iifiach number , temperature 
and pressure p^ , the parameters for the Maxwellian may be written 
down easily. As the rotational degrees of freedom are excited 
completely^ at the upstream temperatures consldeieed (T^ai 295 * > 400 °, 500 ° K), 
the value of rotational temperature (s ) for the upstream gas is the sam^ 
as that of the translational temperature. However, the -vibrational 
excitations of molecules at such low (upstream) temperatures are 


ne^eglble and hence the vibrational temperature 



zero. 


may be taken to be 


Thus we have 

T SB T » T. 5 T w 0 (4»4»1 ) 

t^ r^ 1 ’ v^ 


As a consequence of zero -vibra-bional temperature, the distrlbirtion 
function corresponding to vibrational degrees of freedom, in the 
Maxwellian f^ (See Eqn. 3 • 7.6 ) becomes a Dirac delta function 
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Maxwellian f„ - RETther the translational and rotational teiaperatnres 

(l , T ), the number density (n ) and the velocity (vn) of the gas 
t^ ? D 

at the end of the rotational shock may be fixed using the Eankine- 
Hugoniot conditions for a rotatlon^ly relaxing diatomic gas (y » 7/5)# 
while the vibrational temperature continues to remain at its upstream 
value ( aero ) . 

T-t _ « 


le 


n 


36 

u- 


ci + V 


q - y) 


1 ^ ^ 


^ (1 ^ 

6 m2' 


(4.4.7) 


T. = ; T , * 0* 

^3 "^3 ^ ""3 


Knowing these q.uantities the Maxwellian f^ may be written as 


^^3)^ * '^1 


21 3 

-i— :} 5 CB)ffC(i»|P (4.4.8) 


^ 4Tr 27rkT, 


2iTkT z 


2 k T, 


Maxwellian f. ~ Sie unlfoim downstream condition of the gas may be 
determined by using the iiiree familliar conservation equations 


ni = n4 U4 




" I ’’i’ = "4 (“4 * 

2 -,k„ 

“1 * = “4 * ® 


( 4 . 4 . 9 ) 


fhe values of n. ,u. predicted by these equations are exact for 
4 4 4 

perfect gases. However in the present work the more realistic values 

for n ,u. and I. for nitrogen, considering real gas effects, calculated 
4^4 4 

by Athye and Peng (1962) are used. Recently haw and ^ristov (1969) 
have given more exact values for the dovmatream parameters for normal 
shocks in a tabular form, for both nitrogen and oxygen. !Ehus we have 
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i.e. lim [{-is — exp (- — ] = 


T ->0 irk T 

VI VI 


2k T„ 


6(B) 


(4.4.2) 


Further in. absence of vibrational excitation, the distribution 
g(({) ) for the phase angle <{> is also dropped. Thus we can vorite 
f.j , finally as 


£ tl ',3/2^ > r 


1 


2vkT, 


m(c . i up^ + 
2kTi 


•> 6(0) FCu^p (4.4.3) 


Maxv/ellian f„ ~ The translational temperature (t, ), the velocity 

2 tg 

(ug) and the number density (wg) of the gas at the end of the first 
zone may be fixed by using the Rankine-Hugoniot condition for 
monatonic case (y = 5/3). However, the rotational and vibrational 
temperatures (l , T ) of the gas are frozen at values corresponding 
to their ftee stream values. 


i.e. 


where 


1 . ^ 
r. “ 80 

^2 


(1.^) ^Mr) 

Mf ^ 


n. 


1 Mj 




= 0 


“1=/^ 7/5^1,) 


(4.4.4) 


(4.4.5) 


Thus we have the Maxwellian fg as 


”2 ™ xto *^1 "'■(S-i 

^ ^-2iPr 


2TTkTj"/' -^ ' ■ 2k T 2 


(S“i ^^2)^ ^i^i 


2kTi 


.} 6 (B) F(mj j) (4.4.6) 
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^4 


m 


g^2 2TrkT4 


3/2 


, '^1 ,, 8 s 

WW**" ' »fo,^U.4.1 


o' 


in vshich all the Ihree taaperatures T. f T and 2 have the same 

\ ^4 ^4 

value, 

linally we have the ansata (4»2.1 ) in the tom 


f ^ + Wjfj + (4.4.11) 

where w^, i = ^,4 are the unknown functions ot the distance x 
varying through the shock. The hoxindary conditions on these 

ti c, 

funotionje may^ written as 


Wi s= 1 ; Wg^s w^»f K 0 for X -i. -« 


vn « w„ = w_ = 0; 4 for x ■* + 0 # 

1 23 '4 


(4.4.12) 


Thus the laeleculcs in the shock belong to one of the four states 
described by the four Maxwellians f^,i := 1,4. 

4.5 fee ro SCO pic Quantities : 

The gas dynamic variables like density, velocity and various 
temperatures, inside the shock may be written in terms of the weight 
functions w. , i = 1 ,4 and the parameters in the four Maxwellians used. 
%■ using the definitions for the macroscopic ^^in Eqns. (3*2.2. ) ahd 

( 3.6*8, ) and the expression for f given by Sqn. (4.4.12) we can 

write down the following expressions: 
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Maxwellian - ^Further the translational and rotational temperatures 
(T , T ), the number density (n ) and the velocity (u ) of the gas 

tj r^ 1) 1> 

at the end of the rotational shock may be fixed using the Rankine- 
Hugoniot conditions for a rotationally relaxing diatomic gas (y *= 7/5), 
while the vibrational temperature continues to remain at its upstream 
value ( aero ) . 


le 


36 


(1 + 

M7 


h ' V 


1 

^5 


u 


^1 


^ = T (1 + ^) 


6 

= T, 


m: 


T. = T = T„ ; T » 0. 

^^3 3 ^5 

Knowing these q.viantities the Maxwellian 


(4.4.7) 


f„ may be written as 


3 4Tr airlcTj 


( ^)exp { 

2TrkT3 


Tn(c~j^ 

2k T3 


+ J, w? 


(4.4.8) 


Maxwellian f^ ~ The uniform downstream condition of the gas may be 
determined by using the tliree familliar conservation equations 


ni = n4 U4 


2 k 2 ^^*^4 

* 5 'fl> ° ”4 (“4 * — > 

“1 ^ ^ ■■'1 = * 8 S- T. 

JL m A 4 m ^ 


(4.4.9) 


The values of 5^4jU4»f4 predicted by these equations are exact for 


perfect gases. However in the present work the more realistic values 
for and T^ for nitrogen, considering real gas effects, calculated 

by Athye and Peng (1962) are used. Recently Iiavy" and Aistov (l969) 
have given more exact values the downstream parameters for no mal 


shocks in a tabular fom, for both nitrogen and oxygen. Thus we have 
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w^ — + v/„ 

2 3 u^ 


+ w. 


u. 




)} / (n/n^ ) 


(4-5.1 ) 


T 

_J 

T. 


^ "z ^ 57 *5 557 ^ '"4 


'1 


~ «/> / (iiAi ) 

1 1 ^1 ^ 


Thus, once we knov/ the TEpriation of the weight functions vdth x we can 
calculate the various quantities, using these equations, to obtain 
the sixjck profiles. 

4.6 Moment Equations : 

As indicated earlier, after choosing a form fbr the distribution 
function for the problem, the unknown parameters (weight functions) 
are determined by satisfying appropriate number of moment of the 
Boltamann equation. The molecular property chosen for the moment 
equation should, of course, be relevant to the problem concerned. 

linm the continuity equation (4-1.4), we get a simple algebraic 
equation for w^, i 5 = 1,4 as 

W. + W„ + V/„ + W . =» 1 


(4-6.1) 
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*■4 „ '■O^-VT — — 


(2Bf 


4'“g^2 27rkT4 ^nk?/ \kT^^ 


21c®, 




4.4.10, 


in wMch all the -three terapera-fcures T. « and f have the same 

t, r, Ya 


4 *4 


value. 


ELnally we have the ansatz (4.2.1 ) in the foim 


f = w^f^ + Wgfg + Wjfj + (4.4.11 ) 

where w^, i a 1,4 are the unknown functions of the distance x 

varying through the shock. She boundary conditions on these 
t> c 

functions may^ v/ritten as 

w^ « 1 5 Wg® Wj*! « 0 for X •+ «oo 

(4.4.1 2 ) 

v/^ * '‘‘'^2 ~ ~ ^ 4 *“ ^ X .► +«* 

Thus the molecules in the shook belong to one of the four states 
described by the four Moxwelliane f^,i = 1,4. 

4.5 tlaorosoopio Quantities ; 

Ihe gas dynamic -variables like density, velocity and various 
tempera-tures, inside the shock may be written in terms of the weight 
functions w, , i => 1 ,4 and the parameters in -the four Maxwellians used, 
% using the definitions for the maoroscopic^in Eqns, ( 5.2,2. ) and 

(5.6,8. ) and the expression for f given by Eqn, (4.4.12) we can 

write down the following expressions: 
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Bie other two conservation eqmtions (luomentum and energy) along with the 

qnMrintodal ansatz, Sqn, (4*4.12), do not yield any additional independent 

equatLoas for Thus we need three loore suitable moment equations 

to complete the system of equations for ®or the relaxation flow 

situation, two of the relevant molecular properties chosen are the 

rotational energy and vibrational energy. [Che third choice closely 

2 

follows the bimodal theory and for simplicity the c moment is lised. 

.}C 

acom Sqns, (3.5*5) and (3 *5 *4)) v/e can A-adte down the steady one 
dimensional Maxv/ell transfer equation as 

^ /Q^iii£dT=~/CT!jJ+\jj -v!Jp-i!j)YC'-g.n) * ff^Kd^dx^dr (4.6.2) 
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T 

_JC 

T. 


n n n 

(w^ + — W- + — ^ w. ) 

M 2 3 4 ^ 


a. 


n 


^2 ^2 ^3 ^3 ®4 

{ W. 4* — “ W 4* W. 

^"^1 ^ 2 ^ n^ T^ T 5?^ 4 


+ ^ ‘^ ( -fj'+w, +«2 ^^+>'3^ 

+ *4 ^ ' 


Hr, *1, T, »4 T. 

{w + --2.W +-2^ w + _4._i W }/f — ) 

tw^, + Wg H W^ + w^ ;/^ 


(4.5.1 ) 


^ ‘ ^ 17 *4’ / <“^“1 ’ 


Thus, once we know iihe variation of the wei^t functions vdth x we can 
calculate the various quantities, using these equations, to obtain 
the shock profiles, 

4.6 Moment Equations ; ' 

As indicated earlier, after choosing a form fbr the distribution 
function for the problem, the unicnown parameters (weight functions) 
are determined by satisfying appropriate number of momenis of the 
Boltzmann equation. The molecular property oho sen for the moment 
equation should, of course, be relevant to the problem concerned. 

Brom the continuity equation (4*1.4), we get a simple algebraic 
equation for w^, i « 1>4 as 

w^ + Wg + ^4 “ ^ (4.6.1) 



44 



■'"'■' i T a ♦ v?e get 

SimU-arly lor ♦ = s ^1 a ’ 

d **1 1 4 „ I. . (i- Ji “^5 

; nuck T ) ^ = I «: M 1 ^ 

1*1 

= a ( 2 B)V.. - »>avo 

d! 7 ^ 


( 4 , 6 . 7 ) 


( 4 . 6 . 8 ) 


4 II WiW 3 ii.(|^( 2 Bn 

I bvi(K T^.) ^ ^ 1 

^giag the jneafi free path 

^ is noa-dimeaaioaaliaed by using 

distance ■ ^ sphere molecule). Stirther 

( eeuivaXent ragid sphere 

of «>0 upBlrea- g _ ol lae apow aquations d^y I'O 

• +ViP left band side oi 

^tten as fouotiona ol *1,. ^ 

U^tlonscan thus bs rewltten as 

three eqUat:L 03 :ii 


J J, f!i=i- ll ’< 1*3 Im 3 

,I '^i d? 2 i.j=l 

x=i 


(4*6*9) 


t=l +v,e free stream 

, c are Slvs° 

pners ooeffiolen s affined earlier (see section 

n and the ^rious ratio 

Macb number 
4,4) as follows: 


2 f Gr^) 

^ „ fj:~y +( — y w ^ 

tli 7 Hx ^ 

4 i = '^T ’ 

2 l 7 I 

Csi “ ’ 

/ 1 > 4 , 


( 4 . 6 , 10 ) 


T are given by 

^ ,he ®Uision aoelfleients \ij 
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Ttie other two conservation equations (momentum and energy) along with the 

quadrimodal anaatz, Sqn. (4.4»l2)j do not yield any additional independent 

equations for Thus we need three npre suitable moment equations 

to complete the system of equations for w^*s. Jbr the relaxation fLov/ 

situation, two of the relevant molecular properties chosen are the 

rotational energy and vibrational energy. The third choice closely 

2 

follov/s the bimodal theory and for simplicity the o moment is used. 

Eiom Eqns, ( 3 . 5*3) a^d (3«5»4), we can '.irite dom the steady one 
dimensional Maxwell transfer equation as 


^ . ffjKdiidTjdt 


(4.6.2) 


Substituting for f from Bqn. (4.4.9) in this equation we get 


4 dw. 


dT ^ “x ♦ *1 ' T bjW) 

where I^^( ) is the collision integral given by 


(4.6.3) 


('!')= /C>t»i+i|'’‘-4'i-'(')Y(“g.n)fi fij K dn dij dt 


having the symmetry property 


Iij(4') = Iji(l') 


Putting i|j a 0 ^ we get from (4.6.3) 

•X, ■■ 


(4.6 .4) 


(4.6.5) 


4 ^ . dw. -I 4 ' 

u a . J ,ir »i wj qjfcJ) 


(4.6.6) 
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lij 


I. . Co ) 
nuufX: ^ 


,2 


l2ij = 


mn uu|Xj 




( 4 - 6.11 ) 


I.. ctsC2B)2] 


ran uutXi ^ 

Ilie collision Integrals appearing on the ri^t hand side of the 
Eqns, (4.6.9) are 22 fold integrals. The evaluation of these has 
been shown in Appendix A* Out of the 22 fold integrals only 16 fold 
integrations could be performed analytically for most of the cases. 

The remaining 6 fold integrations have been carried out numerically 
using the i^nte Carlo technique^ the details of which are given in 
Chapter 6 . i9ie set of equations ( 4 . 6 , 1 ) and ( 4 - 6 . 9 ) have been solved 
numerically and shooh profiles have been calculated for four Mach 
numbers (m^ = 15 ^ 20, 25, 30 ) and three upstream pressures, (p^ = 

0.1, 0,01, 0-001 atmosphere), The results obtained are discussed in 
Chapter 7, 


4.7 Trimodal hypothesis for the Shock : 

As the relaxation time for the translational and rotational 

degrees of freedom are nearly of the same order (t^ 'V x^) it Is 

possible to further simplify the quadrimodal formulation by a merger 

of the translational and rotational relaxation sones. This reduces 

the number of ifexwellians used in distribution function by one and 

hence w© have a trimodal distribution function. The molecules inside 

the shock belong either to the upstream equilibrium state, f^ , the 

downstream Maxwellian state, f^, or an intermediate equilibrium state 

with no vibrational excitations, f,. The distribution function in thrs 

3 


case is given by 
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1 2 

Similarly for 'i> ~ 


dwi , 

I, "«» T^. ) 5;^ = j »j iij (j Ji ■«2) 

2 

and for = s ( 2 B) / 2 , we have 


(4.6.7) 


4 dw. 4 . 

I nu(h — = n w, Wj li (i s(2B)2) (4.6.8) 

X-J. 1,J-X • J 

!Hae distance x is non-^imenedonalised by using the mean free path 
Aj of the upstream gas (equivalent rigid sphere molecule). IPurther 
the quantities in the left hand side of the above equations may be 
written as functions of Ji.j f the upstream Mach number. The above 
three equations can thus be rewritten as 


4 dv 7 j -1 4 

,I/kidr°2 . ''i'J 'wj 

i“l i»3=x 


(4.6.9) 


where coefficients are given as the functions of the free stream 
Mach number and the ■variovis ratios defined earlier (see section 
4.4) as follows: 

^i.2 13 . 


cii == (^) 


7 M? '1 


^ 2 i ’ 

7 II? 

5 '^Vi 
C3. = (-A^) 


(4.6,10) 


7 Kf 


and the collision coefficients 3 ^^^ are given by 
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f = wjf^ + w^f3 + (4.7.1) 

whexe ^ ajid are the three vuiknowa weight functions to be 
determined. 

As in the earlier case, a moment method is used to deteimine the 
weight functions. Ibr the present case, in addition to one algebraic 
eqmtion provided by the three conservation equations given by 
Eqn. (4.6*1) two additional moment equations, viz., vibrational energy 
and rotational energy, are used, !I3ius the 1hree equations for the 
present formiaation are given by Bqns. (4.6.I) and (4.6,9) (k = 2,3) 
with 0 in them. The boundary conditions given by Eqn, (4.4*12) 

hold good for the weight fiJnctions and rr^. Solving these 

equations , the shock profiles have been obtained for 6 B/Each numbers 
(M^5=12,15 j 18,20,25,30) and three upstream pressure (p^ssO.1 ,0,01 ,0,001 
atisD sphere). A discussion of these results is given in CGiapter 7. 

She next chapter gives a formulation for the shock structure 
problem considering the dissociation effects. 
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I. .(o^) 

13 ^ X-^ 


L2ij 


nm uu^X]^ 




( 4 . 6 . 11 ) 



Eqns. ( 4 . 6 . 9 ) are 22 fold integrals. The evaluation of these has 
been shown in Appendix A* Out of the 22 fold integrals only 16 fold 
integrations could be performed analytically for most of the cases. 

The remaining 6 fold integrations have been carried out numerically 
using the L’bnte Carlo technique, the details of which are given in 
Chapter 6. The set of equations (4.6.1) and (4.6.9) have been solved 
numerically and shock profiles have been calculated for four Mach 
numbers (M^ = 15> 20, 25 » 30 ) and three upstream pressures, (p^ « 

0.1 , 0.01 , 0.001 atmosphere)* The results obtained are discussed in 
Chapter 7. 


4.7 Trimodal Hypothesis for the Shook : 

As the relaxation time for the translational and rotational 
degrees of freedom are nearly of the same order (t^ 'v- Tj.) it is 
possible to further simplify the quadrimodal formulation by a merger 
of the translational and rotational relaxation zones. This reduces 
the numb ec of ffaxwelllans used in distribution fxmction by one and 
hence w© have a trimodal distribution function. The molecules inside 
the shock belong either to the upstream equilibrivjm state, f ^ , the 
dovmstream Maxwellian state, f^, or an intermediate equilibrium state 
with no vibrational excitations, f„. The distribution function in this 
case is given by 
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SHJCK STRUCTURE WITH DISSQQIAOIQN EgFBOTS 

5*1* Introd-uctlon : 

We have I so far considered the shock wave structiire In a diatomic 
gas with only the vibrational degrees excited to a limited extent so 
that no molecular dissociation takes place. This postulate is valid a 
only moderately high tech numbers. As the tech number increases the 
degree of dissociation in the diatomic gas behind the shock wave become; 
increasingly larger. Thus, when feLj = 20, T^= 300®K, = 5 mm Hg. , 

the degree of dissociation in the case of nltro^n is approximately 
'lig, (5.1 .I ) gives the degree of dissociation behind a normal shock in 
nitrogen for various upstream Mach numbers and pressures ), It 
thereforej necessary to Include the effects of dissociation in the 
analysis of the shock structure at these hi^ Mach numbers. 

The analysis outlined in the Chapters 2, 5, 4- has to be modified 
in two vrays. Mrstly, because of the appearance of a second species, 
the atom, due to the onset of molecular dissociation, an additional 
Boltsmann equation has to be introduced governing 1±ie distribution func 
of the atom species. Secondly, additional collision terms appear on 
the ri^t hand side of the Boltzmann equation for the diatomic species 
representing dissociative and recomblnative collisions corresponding to 
atom - molecule collisions. 'Ehus, ihe problon reduces to the calcuiatr 
of the shock structure in a chemically reacting mixture of two species, 
atoms and molecules. 3he corresponding two Boltzmann equations may be 


written formally as; 
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£ = (4.7.1 ) 

where , yt^ and are the three \Anknown weight functions to be 
determined. 

As in the earlier case, a moment method is used to determine the 
weight functions. Ibr the present case, in addition to one algebraic 
equation provided by the three conservation equations given by 
Eqn, ( 4 . 6 . 1 ) two additional moment equations, viz,, vibrational energy 
and rotational energy, are used. Thus the three equations for the 
present formulation are given by Eqns, (4.6,1 ) and (4.6.9) (k = 2 , 3 ) 
with WgBs 0 in them. The boundary conditions given by Eqn, (4.4.12) 
hold good for the weight functions w^,Wj and ir^. Solving these 
equations, the shock pi'ofiles have been obtained for 6 Mach numbers 
(M^ss12,15j 18|20,25»30) and three upstream pressure (p^ssO.1 ,0.01 ,0.001 
atmosphere). A discussion of these results is given in Chapter 7. 

The next chapter gives a formulation for the shock structure 
problem considering the dissociation effects. 



47 


pf 

rt 


= ff.f] + [f."] - 


recorih 

( 5 . 1 ,’^ 


wliere t and I* axe idie distribution functions for the niolecules 
and atoms respectively, Hie various collision terms on the right 
hand sides of equations (5.1*1) and (5.1.2) are defined as follows: 


i) 


il) 


[f^f] ; is the net gain in number of laolecUleS' por unit time 

to 'non-reacting' collisions between two molecules 
as already obtained in Eqn. (5. 4# 5) as J, 

{f jp] : is the net gain in nwber of molecules per unit volme 
per unit time due to 'non-reacting' collisions 
between a molecule and ao atom. !Chis term can be 
obtained as a special case of (i) above. 


iii) Ef.f},.: is the number of molecules lost per unit volume 
' dis 

per unit time due to dissociative collisions 
between two molecules. 


It) tf.Fl 


■dis • 


is the number of molecules lost per unit volume 
per unit time due to dissociative collisions 
between a molecule and an atom. 

v) the number of molecules ^ned per unit volume 

per unit time due to (three body) collisions 
between two atoms in presence of a molecule 
resulting in recombination. 
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SIDOK STRUCTURE WITH DISSOGIAgEQH EPPliGTS 
5,1« Introduction : 

V/e have, so far considered the shock wave structure in a diatomic 
gas with only the vibrational degrees excited to a limited extent so 
that no molecular dissociation takes place. This postulate is valid at 
only moderately high Mach numbers. As the Mach niimber increases the 
degree of dissociation in the diatomic gas behind the shock wave becomes 
increasingly larger. Thus, when MLj = 20, T^=s 300°K, s 5 mm Hg, , 
the degree of dissociation in the case of nitrogen is approximately 355^* 
(Pig, (5,1.1 ) gives the degree of dissociation behind a noiraal shock in 
nitrogen for various upstream Mach numbers and pressures ). It is, 
therefore, necessary to include the effects of dissociation in the 
analysis of the shock structure at these high Mach numbers. 

The analysis outlined in the Chapters 2,3, 4 has to be modified 
in two ways, HLrstly, because of the appearance of a second species, 
the atom, due to the onset of molecular dissociation, an additional 
Boltzmann equation has to be introduced governing the distribution function 
of the atom species. Secondly, additional collision teims appear on 
the right hand side of the Boltzmann equation for the diatomic species 
representing dissooiative and reoombinative collisions corresponding to 
atom ~ molecule collisions, 33ius, the problem reduces to the calculation 
of the shook structure in a chemically reacting mixture of two species, 
atoms and molecules. The corresponding two Boltzmann equations may be 
written formally as: 
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vl) is the nmaher of molecules ^Ined per unit yolume 

per unit time due to (three body) collisions 
between two atoias in presence of a third, atom 
resulting in recombination, 

vii) fejj} : is the net ^lin in-* number of »tom per .unit volume per 

unit time due to elastic collisions between two 
atoms, Ihis can also be written as a special 
case of (i) above or even directly, 

viii) {f,F) : is the net gain in naaher of -atoms p*a; ► 

unit volume per unit time due to ’non reacting' 
collisions between a moleoule and an atom. !nhie 
can be obtained as a special case of (i). 

lx) number of atoms gained per unit volume per 

unit til&B due to dissociative collisions between 
two moleoules. 

x) ; is the number of atoms ^ined per unit volume per 

unit time due to dissociative collisions between 
a molecule and an atom. 

xi) {J'jl'.fl ; is the number atoms lost per unit volume per unit 

re^corob 

time due to (three body) collisions between two 
atoms in presence of a molecule 

xii) CKF.P} : is the number of atoms lost per unit volume per 

reconflb 

unit time due to (three body) collisions between 
two atoms in presence of a third atom, 

53ie fact that at a dissociative collision one molecule gets con-"'erted into 
two atoms and at a reoombinatlve collision two atoms recombine to form one 
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cs.i.i;' 


££ ={PF}+{f,F}+{ /■,£},. +{£p} -{F.P.f} ,-{FFF) ^ 1 

Pt * tlis * * •^recoml.' '• ■’ * ■'recofab v *** 


where f and S' are the distribution functions for the molecules 
and atoms respectively, The various collision terms on the right 
hand sides of equations (5,1 *1) and (5*1 ,2) are define! as follows: 

i) If»f] 5 Is the net gain in numbesr of molecoiles por raiit time duo 

to •non-reacting’ collisions between two molecules 
as already obtained in Eqn. (3«4#3) as J* 

ii) {fiJ*] : is the not gain in number of molecules per unit ■rolime 

per unit time due to 'non-reacting' collisions 
between a molecule and an atom. This term can be 
obtained as a special case of (i) above. 

iii) [f,f]^j^| is the number of molecules lost per unit volume 

per unit time due to dissociative collisions 
between two molecules. 

iv) s is the number of molecules lost per unit volume 

per unit time due to dissociative collisions 
between a molecule and an atom. 

v) f^**^»£J'Y6Ctftnfa the number ot moleovles gained per unit volume 

per unit time due to (three body) collisions 
between two atoms in presence of a molecule 
resulting in recombination. 
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ffiolecuLe is used to relate the quantities (iii) to (vl) to qualities 
(ir) to (xii). 

Ittxis we have 







(5.1.3) 

recomb “ 

4- ‘ Jreoomb 
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Dissociation and Recombination Hates ; 


fhe collisions resulting in dissociation and recombination 

represented by the following chemical formulae - 





‘^2 

Ag + A + A 

( 5 . 2.1 ) 

re comb. 



, .V- 



A, + A ^ 

2 

A + A + A 

(5.2.2) 


4 


recojab . 
dis ^ 



recomb* 


A + A 


(5.2.5) 


Sie probabilily of sufficiently strongly excited molecule 

into tw5 atoEia as indicated by forwajcd reaction (5*2*3) 
(photo Ci socto/cion) is small ( Zoldovxoh, 196&). !Qie smmu is true 

Of the photo- recombination^ repreSwutcl by reversed reaction (5.2.3). 
Hence these effects are neglected in this analysis* 
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vi) is the number of molecules mined per unit volume 
recomb 

per unit time due to (three body) collisions 
between two atoms in presence of a third atom 
resulting in recombination. 

vii) (PjF} : is the net gain . in'* nvimber of •tom per .unit volume per 

unit time due to elastic collisions between two 
atoms. ®iia can also be written as a special 
case of (i) above or even directly. 

viii) {f,i'} i is the net gain, tn nuraber pf'-atoiiis p«x - 

unit volume per vinit time due to ’non reacting' 
collisions between a molecule and an atom. !Cbis 
can be obtained as a special case of (i), 

ix) {f : is the number of atoms gained per unit volume per 

unit tia© due to dissociative collisions between 
two molecules. 

x) number of atoms gained per vmit volume per 

unit time due to dissociative collisions between 

a molecule and an atom, 

0 .^ 

xi) {]?,]?, f} : is the number atoms lost per unit volume per unit 

recomb 

time due to (three body) collisions between two 
atoms in presence of a molecule 

xii) {I’.I'.F} i is the number of atoms lost per unit volume per 

recomb 

unit time due to (three body) collisions between 
two atoms in presence of a third atom, 
fhe fact that at a dissociative collision one molecule gets con’"'ert'=!d into 
two atoms and at a recombinative oolllaioa two atoms recombine to form one 



50 


!Ehe i^te equations for the chemical reactions mentioned, above may 
be bitten as 


^ dt 


1 


=V tarb ' ■ 


dn 


^2 


^2 = 


1 

4 tjrb 


2 2 
n. + k„ n* n* 

1 * 2 \ ^ S 


2 3 

-kj n. + k^ n. 

dj Aj A A 


(5.2.4) 


(5-2.5) 


where and are the number densities of atoms and molecules 

respectively and the subscripts 1 , 2 refer to the reactions in 
Sqn, (5*2*1) and £qn. (5*2.2). Ihe reaction rate constants may be 
related by the principle of detailed balancing as - 

“'di ■'dj n* 

= __£=, {— ) . ^ CT) 

a S % 

vshere K(3?) is the equilibrium constant f and n* are n^ equilibrium 
values of the number density of atoms and molecules at the given 
tonperatiire* !I3aus, it suffices to know either the dissociation or the 
recombination rate. 


5*2*1. Galoulatlon o$, the ntsspciation Rate : 

Dissociation is assumed to take place when a molecule collides 
with another molecule or an atom with siofficient energsr, ohe frequency 
of inter — molecular collisions has already been calculated in Dqn. (3*4*5)* 
Ihe atom-molecule collision frequency has also been obtained as a special 
case of the inter-^lecule collisions in Eqn, (3»4.4). 

Hie fraction of collisions in which the total ener^r of the colliding 
parlaiers exceeds the dissociation energy is given by Zeldovlch (1966) as 
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molecule is used to relate the c[uatititlee (ill) to (vi) to qualities 
(ix) to (xii). 

Thus we have 






(5.2*5) 


The probabiliiy of sufficiently stron^y exoited molecule 
diesooiating into two atoms as indicated by forward reaction (5*2.5) 
(photo C.i;.:soGic;tion) is very sncll ( Zeldovich, ig6&). The is true 

Of the photo .rocoribinatioa^represuatcd by reversed reaction (5.2.5). 
Hence these effects are neglected in this analysis. 



57 


s 1 exp(- (5,2.i) 

o 

?i4iere is a number characterising the degree of participation of 
the internal degrees of freedom in the dissociation process and D is 
the dissociation energy, She probability of dissociation at an 
energetic collision is known as sterio factor and has to be used to 
write down "the dissociation rate. fihus, tho dissociation rate for the 
reaction Eqa, (5.2.1) may be written using Eqn. (34.4.5) and the above 
two factors as - 


0 


, n 

■ 


(5.2.7) 

Jtor the forw/ard reaction (5.2.2) the dissociation rate may be written 


as 

[f.Fldis'Pl 


X 


D 

RT 


H ■ 
■) ® 


D 

”RT 


/ fFjYC-gi.n)Kj 


dn dt dcj 


where parameters of the molecnle-^nolecule 

and molecule “ atom collisions, respectively^ Bsperimental values 
of these parameters are available in Stupochenko (l967) for nitrogen 
and oxygen^ 


5*2.2. Calculation of the Becombination Rate i 

Becombination is assumed to occur vdien a collision betv^een two 
atoms takes place in the presses of a third particle - either a 
molecule or an atom - which carries part of the bond energy. fhe 
forward reactions in Eqns. (5.2.7) and (5.2.2) correspond to diesociativ 
collistons while the reversed directions refer to recombinative co^ions. 

I II, T. KAr-POl? 

S CWj% SRARY 


i 


, 22658 { 
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(Che rate equations for the chemical reactions mentioned above may- 
be written as 


dn 


^2 1 '-“'A 0 0 

f— b ^dTb “ \ * ’''^1 "A "Aj 




dn 


r ^2, 1 “"A, ,2 1 3 

f ss-i. [— — ]« = -kj n. n* + k„ n. 

^ dt ^2 2 Mt ^2 d A. A r A 


dn 


2 2 


(5.2.4) 


(5-2.5) 


where and are the number densities of atoms and molecules 

respectively and the subscripts 1 1 2 refer to the reactions in 
Bqn. (5.2-1) and Eqn, (5«2.2), lEhe reaction rate constants may be 
related by the principle of detailed balancing as - 


K 


°2 0 
^ = { JL.}2 « K (T) 

^^2 K 

^ h. 


where K(t) is the equilibrium constant} and n-^ are n* equilibrium 

A Ag 

values of the number density of atoms and molecules at the given 
temperature. [Chus, it suffices to know either the dissociation or the 
recombination rate. 


5«2.1. Calculation of the Dissoolatlon Rate : 

Dissociation is assumed to take place vdaen a molecule collides 
v/i-th anothei? molecule or an atom with sufficient energy, 2be frequency 
of inter - molecular collisions has already been calculated in Eqn. (3.4. 3)* 
Ihe atom-ffioleoule collision frequency has also been obtained as a special 
case of the inter-molecule collisions in Bqn, (3.4.4). 

Ihe fraction of collisions in which the total enerar of the colliding 
partners exceeds the dissociation energy is given by Zeldo-vioh (1966 ) as 
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recoiabination rates for both the reaatloas are deterHiiiied*by 
Sqn, (5*2.6) using the principle of detailed balancing so that we have 


tf»P»Plrecomb ~ K(l^ 

( 5 * 2 . 9 ) 

tF.'’.<'JTeco»b = 

Thus, aTi the collision terns on the right hand sides of the Boltzjaann 
equations (5*1.1) and (5*1.2) for the moleavdes and atoms, respectively, 
are determned. 


5*3* Ihe ShocK Structure Iroblem; 


In order to determine the structure of a plane shoclc wave in a 
diesoGiating diatomic gas it is now necessary to solve the Boltzmann 
equations (5.1*0 and (5*1.2) subject to the shock boundary conditions. 
As has been done in the foregoing chapters, a moment method will be used 
to solve the Boltataam equations, Further, the trumodal ansata used 
for the non-dissociating ^s earlier, will again be used hero for the 
molecular distribution function. Ifott-Smith’s bimodal ansatz will be 
used for the atoms. Thus, v/e write 


f(c) = % + W4 f4 

^(Sl) * ^6 * '’5 ^2 


(5*3.1) 


where f are the same ifexwellian distribution functions as defined 
i 

in CSmpter 4 while Hia S', are defined by 

3 


2ifkT 


«(£i - i % 

, „ I, 



2k 

J 


(5.3*2) 


*v5triotly spealcinE this is valid only for the equilibrium situation and it is an 
approximati on for the present case 
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D 

fX-') 


exp(- 



(5.2.i) 


where is a number characterising the degree of participation of 
the internal degrees of freedom in the dissociation process and D is 
the dissooiation energy. 3!he probability of dissociation at an 
energetic collision is known as sterio factor and has to be used to 
v/rite down the dissooiation rate. Thus, the dissociation rate for the 
reaction Eqn, (5.2.1 ) may be written using Eqn. (3*4.3) and the above 
two factors as - 

= Pq ^j(^) ° « ^/f£jY(-g.n)K dn dx (5.2.7) 

R)r the forward reaction (5.2.2) the dissociation rate may be written 
as 

s 

tf.''his='’l ^ fPjYt-Sj.nJKj dn * dSi 


where (p^,s^), (p^,s^ ) are the parameters of the moleoule-flioleoule 
and molecule - atom collisions, respectively* Experimental values 
of these parameters are available in Stupochenlco (1967) for nitrogen 
and oxygen. 


5*2.2. Oalculation of the Eecombihation Rate: 


Recombination is assumed to occur when a collision betvreen two 
atoms takes place in the presence of a third particle - either a 
molecule or an atom - which carries part of the bond energy. The 
forward reactions in Eqns. (5.2.1) and (5.2.2) correspond to dissociativ: 
collisions while the reversed directions refer to recombinative collisions. 


KAKPUR 




1. 1. r. 

central USHARVi 

2265 ^ 


M:' 
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ijG which n^ are the upstream and downstream 

paac" 8 i:io 1 ; 0 p£ fox' iiie atoinic specios givBXi by 


n = 2 a. r.. * U 
o. 1 1 o- 


i 0 , 


T- and n 
1 o. 


2 a4 n. 


C5.3.3) 


^ functions of the distance through the shock. 

°a 2 

5t>r 1iie upsrream tjemperatures upto 5000*^C for nitrogen j the degree 
of dissociation is negligibly small so that for all practical purposes 
vm may assume = 0. Consequently ihe distribution function for the 
atoms becomes a Maxv^ellian 


2TTk7 mCc..i U ) 

F(c ) = 2H>3a4n4 ( exp {- — f 1 -} 

•'i n 2k T 

^2 


(5.3.4) 


with the v-f eight function w^-j temperature T and velocity TJ as the 

b Og 

three unknowns in F, As the masses of a molecifle and an atom are 


comparable the momentum and energy transfer (translational') is efficient 

at a collinion^ ?lth this physical picture^ the Temperature T and the 

^2 

velocity for the atomic species are assumed to be same as that for 

the molecular species* 

i*e* T =; T, and U = u 

°2 * ”2 

vjl 

where T^ and u are given by the Eqns. (4.5*1 ). Thus ^have altogether 
4 unknown weight functions w^, i = 1)3j4,5 which have to be determined 
by the moment method. The distribution functions must satisfy the tTO 
species conservation equations obtained by using the trimodal ansata 


of Eqn. ( 5 . 3 . 1 ) as 


4 4 


I n. u. — = X [ I I 'W ^ I ”i '^5 ^012-* " 
i=l ^ ^ i=l i = l ^ 1=1 


2 04 U4 — = -2Rj 


where 


W2 0, D / » ^i2 ^ ^'^23dis 


{Si5> 


(5 3.6] 
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Otoe recombination, rates for both the reactions are detemiined*by 
Bq.n. (5.2*6) using the principle of detailed balancing so that we have 


I P] recomb 


1 

K(T) 




(5.2.9) 


Thus, all the collision terms on the right hand sides of the Boltzsmann 
equations (5.1.1) and (5.1.2) for the molecules and atoms, respectively, 
are determined. 


5.5, Ihe Shook Structure Problem ; 

In order to determine the structure of a plane shock wave in a 
dissociating diatomic gas it is nov/ necessary to solve the Boltzmann 
equations (5.1.1) und (5.1.2) subject to the shock boundary conditions. 
As has been done in the foregoing chapters, a moment method will be used 
to solve the Boltzmann equations, Purther, the trimodal ansatz used 
for the non-dissociating gas earlier, will a^in be used here for the 
molecular distribution function. Mbtt-Smith's bimodal ansatz will be 
used for the atoms. Thus, we write 

f(£) - wj + W3 £3 + W4 £4 

(5.3.1) 

P(£l) “ ^6 ^ "5 ^2 

where f^ are the same Maxwellian distribution functions as defined 
in Chapter 4 while the are defined by 

"Cjl-iUo/ 

PjCcp-n. oxp {- — — 2 - ) (5.3.2) 

•’ i >" 21< Tq 


*Strlctly speaking this is valid only for the equilibrium situation and it is an 
^approximation for the present case. 
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'Bie integrals D have all been evaluated In Appendix G, 

‘i??o nore znoEiente of the Boltsmann equation for the molecules 
are choaen to obtain a complete set of 4 equations for the 4 unknoTOs, 
in the case of non-diesooiating molecules, we shall choose the 
rotational energy and vibrational energy the Ebltamann equation (5,1.1 )* 
Ihe restating moment equations nay be written as j - 



1=1 i2 


+ 


i=l 


n^Uik 


dw^ 4 

idr= S 


' J J 

X [ II W^w- R.. cisC2B)2) 
ij = l 1 2 11 2 


I (73(26)^)] 

i=l ''i2 2 


" .1 

1=1 i2 


(s ^ 


where 4 ) are the collision integrals already defined in 

Chapter 4, and 




(5.3,8) 


^2 = G. 

fhese integrals are also 22 fold, as in the earlier case and 16 fold 
integration are done analytically in Appendix 0 and the remaining 
integrations are carried out numerically using Monte Carlo technique 
given in Chapter 6* 
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in which (n ,U ,!!? ) and n are the upstream and downstream 

O,’ O, Oj 

parameters for the atomic species gi-ven by 

n = 2 a, n, , U = u, , T = and n = 2 n- C5.3,33 

Oj 1 1 1’ 1 O2 


and U , are functions of the distance through the shock. 

2 2 

Eor the upstream temperatures upto 5000®C for nitrogen, the degree 
of dissociation is negligibly small so that for all practical purposes 
we may assume - 0. Consequently the distribution function for the 
atoBis becomes a Maxwellian 


2TrlcT tn(c,-i U 

FCc ) = 2w5a4n4 ( exp {- 'i—) (5*3.4) 

-1 ^ ^ m 2k T 

°2 

with the v/eight function w^, temperature 3?^ and velocity as the 

three unknowns in P, As the masses of a molecule and an atom are 

comparable the momentum and* energy transfer ( translational') is efficient 

at a collision^ With this physical picture, the Temperature T and the 

■ °2 

velocity U for the atomic species are assumed to be same as that for 

the molecular species* 

i*e, T 55 T. and U = u 

Oj t Oj 

Vs/«. 

where and u are given by the Eqns^ (4.5*1 ). Thus^lxive altogether 
4 unknown weight functions w^, i = 1,3f4>5 which have to be determined 
by the moment method, !Phe distribution functions must satisfy the tv/o 
species conservation equations obtained by using the trimodal ansatz 
of Eqn. (5.3.1 ) as 


4 dw. 44 4 

x=l 1-1 3-1 . i»l 


dWr 


2 «4 1^4>4 djr “ 


where 


W2 = 0. D. . = / » V2 “ 


n^/n 


\6-l-S) 


C5.3*,6) 
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Stoq the foregoiiig discussion, y^a obtain a sot of 4 ordjnaryj 
coupledj non—lijiear^ differential equations yducli aust be solved 
sisjultaneously fox' the weight functions subject to the boundary conditions 


v?j - 1 ; W5 = W4 = \75 = t) for X -> -« 
/ "*4 ® for X -+ «■ 


( 5 . 5 ‘ 9 ) 


53ae Sqns. (5.3*5) and (5*3*'7) are solved nmerically for these boundary 
conditions and shock profiles are obtained for a range of Mioh niHibers and 
three upstream pressures. The results thus obtained are disGUSBod in 
Ghaptof 7. 
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The integrals ^ have all been evaluated in Appendix C.' 

Two more momenta of the Boltzmann equation for the molecules 
are chosen to obtain a complete set of 4 equations for the 4 unknowns. 

As in the case of non-dissociating molecules, we shall choose the 
rotational enej^gy and vibrational energy the Boltzmann equation (5.1.1 )• 
The resulting moment equations may be i,vritten as : - 


VlW- -t 10^ 1-7 

I hiU.kT ll ,1 "V's^o./l 

1=1 1 1>3 = 1 •’ 1=1 i 2 

+ X [.n_ 1. 

1 ^ ^ 1 X” X 

4 dv;. 4 1 -7 19 

.1 "i“lk M , 'Vj * .1 WjWsI^ cfC2Ii}2) 

i=l 1 J ^ 1=1 12 

X ill_ WiWj (“S(2B)^) 

.1 'Vs P^o. Cjs(2B)2)] 




where I. . ( <? ) are the collision integrals already defined in 
ij 

Chapter 4, and 

R..C$) = / dT 

B.o. 2C^) = / 4>[fi,F2ldis ' ^5.3.8) 

W2 = 

These integrals are also 22 fold, as in the earlier case and 16 fold 
integration are done analytically in Appendix C and the remaining 
integrations are carried out numerically using Monte Carlo technique 
given in Chapter 6d • . ^ ^ 



OHAMSl - 6 
EVi^IiUATIOH 0^* U? TBPtR AT.R 


6,1 Jjitroductloa : 

31ie collision intsgrals in the moment equations encountei-ed in 
the preceding tm chapters are 22 fold end in most of the cases 16 
fold integrations were carried out analytically. It ms not possible 
to do the remining 6 fold integration analytically and henoe a numerical 
approach was resorted to. In any of the standard techniques of multiple 
quadrature, the number of points required to evaluate an integral goes 
up in geometrical progression as the order of the manifold integral, 

'SxLs exponential increase never occurs in case of Monte Carlo evaluation. 
Ihus for the numerical integration of manifold integrals this method is 
Tory economical, Ihe Jtonte Carlo technique gives a reasonable estimate 
of the integral and errors in the integration can also be estimated during 
the process of integration, 

^P-ts Oorlo Integration ; 

A typical integral which has to be evaluated may be written 
formally from Appendix A or 0 as 

I = / /dB/’^de dB, (6.2.1) 

0 0 0 0 0 0 

fbr a Monte Carlo evaluation of this integral the limits of integration 
are changed to 0 to 1 for all the variables by using linear transformations 
for then* Qius we have the integral 
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From, the foregping discussion we obtain a set of 4 ordinary, 
coupled, non-linear, differential equations which must be solved 
simultaneously for the weight functions subject to the boundary conditions 


- I ; W 5 = tW 4 = Vir 5 “ for x -> -«> 

v/j = W 3 = 0 ; W 4 ® for x -+• « 


(5.3.9) 


G}he Eqns. (5.3.5) and (5»3.7) are solved numerically for these boundary 
conditions and shock profiles are obtained for a range of Mach numbers and 
three upstream pressures. The results thus obtained are diacu^od in 
Chapter 7 . , 
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1 


dx/^dy/^dz/^du/^dv/^ M(x,y ,z ,u,v,w) dw 
0 0 0 0 0 0 


( 6 . 2 , 2 ) 


where 


H(x,y,z,u,v,w} 



max 


F (BjBj 


(6.2.5) 


The method consists of evaluating the integrand H at U randomly 
chosen points i = 1 ,1T lying with-iii the domain 

of region of integration and finding average of these ¥ values. 0x13 


average value itself is an estimate 


i,e* 



I HCx-.yj 


I „ of the integral, 
e.¥ 


,Z^>Ui.Vi,W.) 


( 6 , 2 ^ 4 ) 


Sbr this calculation the points choaen as random 

numbers belonging to 6 independent series of unxfomnly distributed random 
numbers between 0 to 1 . For the numerical computational purposes the 
random digits available in the Rand Corporation feble have been used in 
this thesis. 


6.5 Error Estimation t 


Ihe variance of H is defined by 

dT = h 2 - I"' (6.5.1 ) 

0 

where 

dx = dx dy dz du dv dw, H = / 

0 

and the standard deviation i® given by 

= a/t^ 


(6.3.5) 



CHAPTER 


6 


NUMERIC A1 EVAIUATIOH OF IKTEGRAIiS 


6,1 Introduction : 

The collision, integrals in the moment eqmtions encountered in 
the preceding two chapters are 22 fold and in most of the cases 16 
fold Integrations were carried out analytically. It was not possible 
to do the remaining 6 fold integration analytically and hence a numerical 
approach was resorted to. In any of the standard techniques of multiple 
quadrature, the number of points required to evaluate an integral goes 
up in geometrical progression as the order of the manifold integral. 

This exponential increase never occurs in case of Monte Carlo evaluation. 
Thus for the numerical integration of manifold integrals this method is 
very economical. The Monte Carlo technique gives a reasonable estimate 
of the integral and errors in the integration can also be estimated during 
the process of integration. 

6,2. Monte Carlo Integration ; 

A typical integral which has to be evaluated may be written 
formally from Appendix A or 0 as 

2tt 2Tf 

I = / d<|>/.d.j.J de/ deJ dfi/ dB, (6.2,1) 

0 0 0 0 ^0 0 i i i 

Ibr a Monte Carlo evaluation of this integral the limits of integration 
are changed to 0 to 1 for all the variables by using linear transformations 
for then* SOius we have the integral 
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The following table given by Kahn {l965) gives the probability that the 
deviations of I from I will esceed + A S . 

eu Q 

^ Hx^bability 

0.6745 0.5000 

1.0000 0.3173 

2.0000 0,0455 

3.0000 0.0027 

4.0000 0.0001 

It is seen from this table that deviations exceeding + 3^ will be 

frequent whereas deviations esceeding + 2 3^ or + 5 3^ are not so 
coEsaon, Ihe quantity is used as a measure of the sampling error 

of the estimate !Ehus it is possible to limit the error to a small 

value by taking a large number of sample points II, 

In actual computations of i estimates of and I^ are 
substituted in the expression for (Eqn. 6,3,1 ) 

l-e. » C<r5«,. = - 4 (6.3.4) 

where 

— j N 

" H ^ ^ (6,3,5) 

ial 

In otherwords, since the value of is also estimated by using 

the Honte Carlo technique, strictly speaking, the errors involved 

((j2 „ should also be considered before using the value of 

eiN eh 

instead of , However, this is not done inview of the fkct that, 
a 10-15^ error in estimating the variance itself affects the overall 
error criterion negligibly. Ibr example - let 10^ error iie predicted in 
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I = /^ dx/^dy/^dz/^ du/^dv/ H(x,y,z,u,v,w) dv; (6*2«2) 

0 0 0 0 0 0 

where 

^K^»y » 2 ,u, v,w) “ 417 j*!* (6»2«3) 

Ihe method coneists of evaluating the integrand H at N randomly 
chosen points (x^,y^,z^,u^,rj^,vv^), i » 1 ,ll lying vd.th-in the domain 
of region of integration and finding average of these IT values. This 
average value itself is an estimate of the integral, 

1 

i,e. I Igj,, J (6,2.4) 

i^l 

jBbr this calculation the points chosen as random 

numbers belonging to 6 independent series of uniformly distributed random 
numbers between 0 to 1. libr the numerical oomputational purposes the 
random digits available in the Rand Corporation Table have been used in 
this thesis. 

6.5 Error Estimation t 

The variance of H is defined by 

a2 ^ jl dx = (6.3.1) 

0 

where 

dx « dx dy dz du dv dw, = /^ dr (6.3»2) 

O', 

and the standard deviation is given by 

= a/*^ (6.3.3) 
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the estimation of Integral I, then a + 15^ error in this error 
estimation allows the overall error in I to take a value between 
8.5^ and 11.5^* 


In pratical calculations the choice of Bf depends on the error 
that can be allowed. Lordi and Mates (l970) evaluated 6 fold 
integrals encountered in their analysis, using Monte Carlo technique 
and obtained satisfactory estimates for the muLti-dimensional integrals. 
They us^ 400-500 trials for different oases depending on the error 
criterion. Ihree typical plots of the estimates and versus 

17 are given in Fig, (6.3.1) for the case of three collision integrals 
for f Rotational aiergy and Vibrational Energy. It is found that 
600 trials are required to limit the overall error to _+ 10^ in the 
last two oases* However, for the raoment integral as many as 1000 
trials were necessary to restrict the eriror to Jr 10^* In a few oases 
of the 0^ jaioment integrals^ the errors were as high as 26% for 1000 
trials . 

Ibr evaluating a typical six fold integral (f) or 
given in Appendix A or G it takes about 3 minutes on the IM 7044, 
with 600 Monte Carlo trials. Fov the cases when i the integrals 

are 7 fold and the same procedure is adopted for evaluating this 
integral. !Ihe comput.-r time taken for these cases is about 3^ minutes. 

As the integrations are carried out approximately the shock profiles 
calculated may also be expected to be approxunate to the same extent. 

In the next, chapter a discussion on the results and conclusions are 


given 



58 


She following fable given by Kahn (1965) gives the probability that the 

deviations of I „ from I vd,ll exceed + X S . 

elf — 0 


X 

Probability 

0.6745 

0.5000 

1.0000 

0.3173 

2.0000 

0.0455 

3.0000 

0,0027 

4.0000 

0.0001 


It is seen from this table that deviations exceeding j+ will be 
frequent whereas deviations exceeding _+ 2 or + 3 are not so 
co3Hmon. The quantity is used as a measure of the sampling error 
of the estimate Thus it is possible to limit the error to a sniall 

value by taking a large number of sample points if. 

In actual computations of $ estimates of and are 

substituted in the expression for (Eqn. 6.5«l) 


i.e. 

where 




eN 


“ 2 * 1 ^ 
1=1 


(6.3.4) 


(6.3.5) 


In othexwords, since the value of Is also estimated by using 
the Monte Carlo technique, strictly speaking, the errors involved 
(o^ - e|j^j) should also be considered before using the value of 
instead of , Hjwever, this is not done inview of the feet that, 
a 10-15^ error in estimating the variance Itself affects the overall 
error oriterion negligibly. ibr example - let 10^ error l>a predicted in 



CHAH’ER .7 

DISCUSSIONS AND GOHGDUSIONS 


7»1* Shock Structure in a Vibraidonally Relaxing Gas 

An attempt was made in Chapter 4 to investigate the effects of 
vihrational relaxation on the structure of a shock wave and two 
formulations were presented for diatomic gases. Both the trimodal 
and the quadrimodal theories are extensions of the Mott Staith himodal 
approach and the trimodal theory of Haight for a rotationally relaxing 
diatomic gas. 'lire various shock profiles for nitrogen predicted by 
the trimodai theory of Chapter 4 ere shown in Pigs. 7.1.1 to 16 and 
the results given by the quadrirodai approach are presented in Jigs. 
7,2.1 to 9. 

7*1*1. Shock Profiles (^Crlmodal ICheory): 55ie trimodai calculations 
are carried out for the upstream Mach numbers ranging from 12 to 30 ai 
for the upstream pressures of 0.1, 0.01 and 0.001 atmosphere. Ihe 
density profiles for these Iiflach numbers are given in HLge, 7*1*1 to 6 
It is seen from these profiles that for a given upstream Mach number 
the shock transition takes place faster with the decrease of the 
upstream pressure, whereas for a given upstream pressvire the shock 
becomes thinner with the increase of the upstream Mach number. 

!I!he corresponding temperature profiles are given in Jigs* 7.1* 
to 12. The dependence of the width of the aone of temperature varia" 
on the free— stream conditions, is same as in the case of density. In 
all the cases, the translational temperature initially rises at a fas' 
rate to a maximum and then decreases gradually to its down-stream 
e<juiHbnum value, thus exhibiting a pionounced temperature peak. Tht 
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the estimation of Integral I, then a jf 15^ error in this error 
estimation allows the overall error in I to take a value between 
8.5^ and 11.5^ 

In pratical calculations the choice of N depends on the error 

that can be allowed 4 Ijordi and Mates (1970) evaluated 6 fold 

integrals encountered in their analysis, using Monte Carlo technique 

and obtained satisfactory estimates for the multi-dimensional integrals. 

They used 400-500 trials for different cases depending on the error 

criterion. Three typical plots of the estimates and versus 

N are given in lig, (6.3.I) for the case of three collision integrals 
2 

Tor o' t Rotational ISnergy and Vibrational Energy. It is fomd that 

600 trials are required to limit the overall error to + 10^ in the 

2 

last two cases. Ifowever, for the o^ moment integral as many as 1000 
trials were necessary to restrict the error to _+ 10^. In a few oases 
of the c moment integrals, the errors were as hi^ as 26^ for 1000 

Jim 

trials. 

Ibr evaluating a typical six fold integral 1^^ (f) or D^^{^)or 
given in Appendix A or 0 it takes about 3 minutes on the IBI 7044, 
v/ith 600 Monte Carlo trials, Ibr the cases when i V j, The integrals 
are 7 fold and the same procedure is adopted for evaluating this 
integral. The ooiaputv'r time taken for these oases is about 3s minutes. 

As the integrations are carried out approximately the shock profiles 
calculated may also be expected to be approximate to the same extent. 

In the neact, chapter a discussion on the results and conclusions are 


given 
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the rotational profile lags behind the translational one* it shows a peak 
for Maoh numbers above 20* Shis overshoot is very much less pronounced 
and it occurs slightly later when compared to that of translational 
temperature. however the vibrational temperature monotonically increases 
xo its down-stream equillbriimi value and further it is found to lag behind 
the other two temperatures. Ibis lag in the vibrational tempera tuire is 
significant for Mach numbers above 20* 

Hhe above pattern of the three temperature distributions may be 
explained as follows: 2he translational degrees being the fastest of 
the three, the cxjrres pond lag temperature attains its equilibrium value 
(monatomio) first with- In a short distance ( 'u This is also true 

of the rotational degrees of freedom which is slower than the translational 
but faster than the vibrational degrees of freedom. Hnally the slowest 
of the three namely the vibrational degree absorb energy from the faster 
ones and relax to the equilibrium state gradually. The degree of lag of 
vibrational temperature at lower Mach numbers is small mainly due to the 
classical nature of the vibrational mode assumed. 

In Tigs. 7.1.13 and 14 are shown the density and temperature 
profiles with the upstream temperature as a parameter. It is seen from 
these profiles that for a given , an increase in the upstream temperature 
reduces the extent of the shock due to lii^er absolute temperatures inside 

the shook. 

A typical variation of the wei^t functions w^,w^ and w^ is given 
in ELg. 7.1.15. The functions and w^ show a monotonie behaviour 
within the sho ck, whereas as expected lias a maximum inside the 
shook with aero values at the two ends, Ibr > 20, w^ is found to 

This fact is reflected in the non-existence of a peak 


take small values. 



CHAPfEa -7 

Discussiosrs AND COITGDUSIONS 


7»1» Shock Structure la a VlbratLonally Relaxing Gas 

An attempt was made in Chapter 4 to investigate the-effeots of the 
vibrational relaxation on the structure of a shock wave and two 
formulations v/ere presented for diatomic gases. Both the trimodal 
and the quadrimodal theories are extensions of the Mott Sbiith bimodal 
approach and the trimodal theory of Haight for a rotationally relaxing 
diatomic gas. fflhe various shock profiles for nitrogen predicted by 
the trimodal theory of Chapter 4 are shown in Pigs, 7*1 *1 to 16 and 
the results given by the quadrimodal approach are presented in Pigs, 

7.2.1 to 9. 

7.1.1. Shock Profiles (itimodal Theory): Bie trimodal calculations 
are carried out for the upstream Mach nmbers ranging from 12 to 30 and 
for the upstream pressures of 0.1, 0.01 and 0.001 atmosphere. The 
density profiles for these Mach numbers are given in -liligs, 7*1.1 to 6. 

It is seen from these profiles that for a given upstream Mach number 
the shook transition takes place faster with the decrease of the 
upstream pressure, whereas for a given upstream pressure the shock 
becomes thinner with the increase of the upstream Mach number. 

The corresponding temperature profiles are given in iigs. 7.1.7 
to 12. The dependence of the width of the zone of temperature variations 
on the free-str earn conditions, is same as in the case of density. In 
all the cases, the translational temperature initially rises at a fast 
rate to a maximum and then decreases gradually to its down-stream 
equilibrium value, thus exhibiting a pronounced temperature peak. Though 
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in the rotational temperature and a faster approach of the vibrational 
temperature to the rotational temperature. Ehie is perhaps attributable 
to the lower population of the molecules characterised by f^ inside the 
shock, .as the dov/n-etream condition (f^) is attained faster due to the 
classical behaviour of the vibrational degrees of freedom assumed. Ihls 
is not so at higher iilach numbers as this assumption regarding the vibratioi 
degrees of freedom is valid to a better extent. 

A com|sa,rision of 'the, theoretical- profiles with the e^cperimental el 
profiles given by Hobben and lalbot (l966) is made in Elg. 7*1»16 for 

12.9. ihe density profile predicted by the trimodal theory oorapare 
very well with the experimental profiles over more than half the distance 

throu^ the shock and differ slightly over rest of the distance. But 

.t 

for a shift in the calculated rotational temperature prof ile, also : compares 
very well vdth the experimental profile. 

7.1.2 Shock Profiles (Quadrimodal Theory): 'Che denslly profiles for th* 
quadrimodal case are given in Tigs. 7.2.1 to 3 and the temperature profilee 
in Pigs, 7.2.4 to 6. The calculations were done for 11 ^- 15, 20 and 25 and 
for the upstream pressures indicated earlier. (Bor M^«=30j the shock profile 
ai*e obtained at p^== 0,1 atmosphere and are shown in Pig. 7.2, 7*) Tlie 
profiles In this case are qualitatively same as in the trimodal case. IThe 
translational temperature overshoots are some what hi^er in this case 
than for the earlier theory. 

In iig, 7.2.9 a comparison of the various profiles given by both 
the theories is made (while comparing p w 0.5 is matched for both the 
cases where p - ( p - p^) / ( P4 - Pi )* density profiles match olaselj 

upto a certain distance from the upstream whereas, temperature profile 

compare fairly closely but for the over shoot in 1^. Pig. 7.2,8 shows a 
typical toon of the weight functions w^. Ehe flmetions nnd 
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the rotational profile lags behind the translational on.e>it shove a pealc 
for Mach numbers above 20. Ihis overshoot is very much less pronounced 
and it occurs slightly later when compared to that of translational 
ttsiiperature . hovfever the vibrational temperature monotoxiically increases 
to its down-stream equilibrium value and further it is found to lag behind 
the other two temperatures. OMs lag in the vibrational temperature is 
significant for Mach numbers above 20. 

She above pattern of the three temperature distributions may be 
explained as follows: Ihe translational degrees being the fastest of 
the three, the corresponding temperature attains its equilibrium value 
(monatomic) first with- in a short distance ( ''< Xj ). IHiis is also true 
of the rotational degrees of freedom which is slower than the translational 
but faster than the vibrational degrees of freedom. Finally the slowest 
of the three namely the vibrational degree absorb energy from the faster 
ones and relax to the equilibrium state gradually. The degree of lag of 
vibrational temperature at lower Mach numbers is small mainly due to the 
classical nature of the vibrational mode assumed. 

In ligs. 7.1.13 and 14 are shown the density and temperature 
profiles with the upstream temperature as a parameter. It is seen from 
these profiles that for a given , an increase in the upstream temperature 
reduces the extent of the shock due to higher absolute temperatures inside 
the shock, 

A typical variation of the weight functions and w^ is given 

in Ilg. 7. 1.1-5. The functions and w^ show a monotonio behaviour 
within the shock, whereas Wjf as expooted has a. maxiiaum inside the 
shook with zero values at the two ends. Ibr >- 20, w^ is foxmd to 
take small values. Ihis fact is reflected in the non-existenoe.of a peak 
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have monotocic behaviour as in the earlier ease and and possess 

a mscd-iHUEi vsdthin the shock} having zero values at the two ends. 

7,1.3. Shock Thiclmess: The following three shock thicknesses are 

calculated for the vihrationally relaxing diatoinio gas using the profiles 

given by both the theories} density thickness ( ), vibrational tanpera- 

,ure thicloiess ( Ay ) and the rotational temperature thickness ( Aj^ )* iu 

all the cases the maximum slopes of the corresponding profiles are used. 

These thickness indicate the extent of the shock when the actual profiles 

are replaced by linear profiles with the corresponding maximum slopes. 

The variation of the reciprocal shock thicknesses vifitb the upstream Mach 

number is shown in ligs, 7.4.1 to 3* for nondiraensionalising the shock 

thickness a length L* adopted by Itobben and Talbot (l966) is used instead 

of the free-stream • It is seen from the :e.gures that the density 

thickness is lower than the two temperature thicknesses at higher Mach 

the 

numbers « -This -i a expected due to the presence o:^ rotational and vibrational 

relaxation zones' at higher Mach numbers. The quadiimodal approach predict 

lov/er values for the shook thicknesses. The experimental values of 

and A at M,= 12.9 given by Robben and Talbot (l966) are also shown 
R * 

for oomparision. The theoretical Aj^ compares well with the experimental 
value at ?.I^=12.9. finally} for a given the shock thicknesses are found 
to decrease with the upstream pressure. 

7,2. Shock Structure in a Bjssociating Diatomic Gas c 

The analysis given in Chapter 4 ms further extended to include 
the effects of dissociation in Chapter 5. shook profiles thus 

obtained are presented in figs, 7*3.1 to 9. 

7.2.1, Stock Profiles with Dissociation Effects; In Elgs, 7.3.1 to 3 
the degree of diseooiatioh, total density and noleoular deasity profUes 
are given for !l, 15. 20 and 25 respectively. Bie profiles ora also given 
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in the rotational temperature and a faster approach of the vibrational 
temperature to the rotational temperature, 5his is perhaps attributable 
to the lower population of the molecules charaoterised by inside the 
shock, .as the down-stream condition (f^) is attained faster due to the 
classical behaviour of the vibrational degrees of freedom assumed, Kiis 
is not so at higher tlach numbers as this assumption regarding the vibrational 
degrees of freedom is valid to a better extent, 

A comparision of 'the.thieoretical ... profiles with the experimental shook 
profiles given by Eobben and Talbot (l966) is made in 5ig, 7«1»16 for 
M.jsa 12,9, The density profile predicted by the trimodal theory compare 
very well with the experimental profiles over more than half the distance 

throu^ the shock and differ sligJitly over rest of the distance. But 

it 

for a shift in the calculated rotational temperature profile, also: compares 
very well vd.th the experimental profile, 

7.1.2 Shock Profiles ( CJuadrimodal Theory): The densily profiles for the 
quadrimodal case are given in Pigs, 7.2»1 to 3 and the temperature profiles 
in ELgs. 7.2,4 to 6, The calcvilations were done for 15, 20 and 25 and 
for the upstream pressvires indicated earlier, (Bor M.j=30, the shock profiles 
are obtained ai p.|= 0,1 atmosphere and are shown in Big, 7,2.7.) Eie 
profiles in this case are qualitatively same as in the trimodal case. The 
translational temperature overshoots are some what higjaer in this case 
than for the earlier theory. 

In Big, 7,2,9 a comparison of the various profiles given by both 
the theories is made (while comparing p s* 0.5 is matched for both the 
cases where p “ (p - p^) / (p^ •• )• 2he density profiles match closely 

upto a certain distance from the upstream whereas., the. temperature profiles 
compare fairly closely but for the over shoot in T.^.* Big. 7*2,8 shows a 
typical vaiiation of the vreight functions w^. The functions w.j and 
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for M^sa 25 wi1ii as parameter in 7*5 »4» I!he eorresponding 
teGnperature profiles are given in Figs. 7.3.5 to 8. In all these 
profiles the total density is found to increase monotonioally whei^as 
the nolecular density has a maxLiiium inside the shock. !Ehe difference 
of these densities is related to the degree of dissociation by ct = 1 - Pi^/p 
Ihe study made by Itorranee (l96l) indicates the possibility of maxima 
for both p and « . However in the present case inspite of a peak in 

the profiles for ^ the profiles for degree of dissociation have a 
mono tonic behaviour, 

Eor a low a ('V 0.2) the three temperatures and are 

found to equilibrate in a narrow region v/here the dissociation is 

practically absent. This region is followed by a very broad clearly 

distinguishable tail in which the degree of dissociation attains its 

the 

equilibrium value. In this broad regioa|/chan.ges in p , ’’t- and 
are very gradual whereas the molecular density and the degree of dissocia- 
tion change significantly (see lig, T.3.1, 2, 5 and 6), 

For a high ts('u ,9) the two regions lose their identity and 
merge together. However there exists a broad sono in which the throe 
temperatures attain nearly their equilibrium value. Erls is followed by 
a tail whose width is also of the same order. In this region the only 
quantities which vary significantly are the two densities and the degree 
of dissociation, (see Bigs, 7,5.3 and 7 ). 

Finally for a moderately high a ( '\i 0.6) there exists a broad 
region in which the three temperatures and the total density attain 
aliaost their equilibrium value. This is followed by a very broad region 
in which the mol8c\flar density and the degree of dissociation attain 
their equilibrium values as in the case of low p » In otherwords the 
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have mottotonic behaviour as in the earlier case and and possess 

a maximum within the shock, having zero values at the two ends. 

7.1,3. Shook Thickness! The following three shock thicknesses are 

calculated for the vibrationally relaxing diatomic gas using the profiles 

given by both the theories: density thickness ( Ap ), vibrational terapera— ' 

ture thickness ( Ay ) and the rotational temperature thickness ( ). In 

all the cases the maximum slopes of the corresponding profiles are used. 

These thickness indicate the extent of the shock when the actual profiles 

are replaced by linear profiles with the corresponding maxi mum slopes, 

The variation of the reciprxical shock thicknesses with the upstream Mach 

number is shown in Tigs. 7,4.1 to 3, Por nondimensionalising the shock 

thickness a length L* adopted by Robben and Talbot (l966) is used instead 

of the free-atream * It is seen from the figures that the density 

thickness is lower than the two tcaaperatvire thicknesses at higher Mach 

the 

numbers .'iMs -is expected due to the presence o:^_ rotational and vibrational 
relaxation 2 sones at higher Mach numbers. The quadrdmodal approach predicts 
lower values for the shook thicknesses. The experimental values of Ap 
and Aj^ at M^= 12,9 given by Robben and Talbot (l966) are also shown 
for comparision. The theoretical An compares well with the experimental 
value at M^=12.9» linaily, for a given the shock thicknesses are found 
to decrease with the upstream pressure, 

7,2, Shock S tructure in a Dissociating Diatomio Gas t 

The analysis given in Chapter 4 v/as further extended to include 
the effects of dissociation in Chapter 5* Uho shock profiles thus 
obtained are presented in Mgs. 7,3,1 to 9. 

7*2*1, Shock Profiles with Dissociation Siffeotss In Mgs, 7,3,1 to 3 
the degree of dissociation, total density and molecular density profiles 
are given for 15, 20 and 25 respectively, lobe profiles are also given 
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merger of the two zones vri.ll haTe just begun, (see Elge. 7.3,4 and 8), 
ihese figures also indicate the effects of the upstream pressure 
on the various profiles. As found in the case of shock structure without 
dissociatiouj in this case also^ with the decrease of p.| the extent of the 
shock decreases, 

7*2.1, Shock Hmckness with Slssociation Effects: She definition of 
shock thickness should have some bearing to the kind of relaxation 
phenomena predomlnent, Ibr the shock in dissociating gas the degree of 
dissociation of the gas has to be incorporated in the definition of a 
thickness, Mor the present case In the absence of any other definition 
in the literature, a shock thickness A is defined as a distance from 
the begining of the shock transition to a point inside the shock whore 
the degree of dissociation is 90 % of the down-stream equilibrium value, 
(a similar definition has been adopted in the case of ionising shock 
waves in mono tonic gases). Ihe position where shock transition starts 
is of course not very well defined and to this extent there is certain 
aESount of arbitrariness in the value of A chosen. 

5br a given upstream pressure the equilibrium degree of dissoci- 
ation increases with Ifae Mach number. Rxrther with the increase of 
equilibrium degree of dissociation the shock thickness decreases, liis 
fact is reflected in Hg- 7*4.4 which shows the variation of a and 
1/a with M.J for nitrogoi 


7 
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for M^=i 25 with as pairameter ia lig. 7»5»4» ttie oorrespondiag 
temperature profiles are given in Pigs. 7«3*5 to 8. In all these 
piofiles the total density is found to increase monotonically whereas 
the molecular density has a maximum inside the shock. The difference 
of these densities is related to the degree of dissociation by a = i « p^/p 
The study made by Itorranoe (l96l) indicates the possibility of maxima 
for both p and a , However in the present case inspite of a peak in 
the profiles for ^ the profiles for degree of dissociation have a 
mono tonic behaviour. 

Ibr a lev/ a (a- 0.2) the three temperatures and. Ty are 

found to equilibrate in a narrow region where the dissociation is 

practically absent. This region is followed by a very broad clearly 

distinguishable tail in which the degree of dissociation attains its 

the 

equilibrium value. In this broad regior^changes in p , and 

are very gradual whereas the molecular density and the degree of dissooia.-- 
tion change significantly (see llg. T«3.1j 2, 5 and 6). 

Ibr a high ctC'v .9) the two regions lose their identity and 
marge together. However there exists a broad zone in which the three 
temperatures attain nearly their equilibrium value. This is followed by 
a tail whose width is also of the same order. In this region the only 
quantities which vary significantly are the two densities and the degree 
of dissociation, (see Pigs, 7.3.3 and 7). 

Finally for a moderately high a( 0.6) there exists a broad 
region in which the three temperatures and the total density attain 
almost their eqmlibrium value. This is followed by a very broad region 
in which the molecular density and the degree of dissociation attain 
their equilibrium values as in the case of low a • ^ otherwords the 



7 *3 . Opnc lug ion s ; 

For the comparison purposes hardly any experiaGntal shock 
profile data for nitrogen in the Mach number range considered is 
available. Thus no defini:e conclusion can be made as to tho validity 
of the theory, Rowever the two foronlations-quadriiTiodal and triKiodal- 
given in Chapter 4 are capable of predicting the shook structure in a 
vibrationally relaxing diatomic gas. Out of these two, the quadrimodal 
theory is better than the trimodal one especially for the Hach number 
range 15 to 20, Further the analysis presented in Chapter 5 includes the 
dissociation effects in the gas and hence is valid for higher Hach number 
upto 50 . for > 30 , the ionisation and photo-recombination of the 
atoms and photo-dissociation of the moleoulos become appreciable, These 
effects are not included in the present formilation. 

As particular cases of the quadrimodal theory, for vihrationally 
relaxing diatomic gas given in Chg,pter 4j it is possible to recover, both 
the trimodal theory proposed by Haight for rotationally relaxing diatomic 
gases (w^=0) and the Mott Smith bimodal theory (w^«w^bO) for inert gases. 

Further to enable to carryout meaningful calculations for shook 
structure, the parameters b, a and s for axially elastic sphoroc^^'linder 
models coo^sponding to real diatomic gas molecules should be available. 
In the present case the values indicated in the Table 2,1.1 wore used. 

The shock structure analysis given in Chapter 5 for the 
dissociating gas may further be extended for a mixture of atoms caid 
molecules of the same -gas. In such a case the upaoreom a.tcmie number 
density is nonaero and hence the full bimodal ansata for the distribution 
function for atoms has to be retained. Correspondingly, appjxspriate 
number of relevant moments of the two Boltsmann equations have to bo 
consndered in addition to the equations of Chapter 
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merger of the tv/o zones vri.ll have just begun, (see ELgs. 7 •3*4 and 8), 
Shese figures also indicate the effects of the upstream pressure p.| 
on the various profiles. As found in the case of shock structure vri-thout 
dissociation, in this case also# with the decrease of the extent of the 
shock decreases# 

7.2.1. Shock Ihiokness vdth Dissociation Effects: fhe definition of 
shock thickness should have some bearing to the kind of relaxation 
phenomena predondnont# Ebr the shock in dissociating gas the degree of 
dissociation of the gas has to be incorporated in the definition of a 
thickness. Ibr the present case in the absence of any other definition 
in the literature, a shock thickness A is defined as a distance from 
the begining of the shock transition to a point inside the shock v'lbere 
the degree of dissociation is 90 ^ of the down-stream equilibrium value. 
(a similar definition has been adopted in the case of ionising shock 
waves in monotonic gases). The position where shock transition starts 
is of course not very well defined and to this extent there is certain 
amount of arbitrariness in the value of a chosen. 

3br a given upstream pressure the equilibrium degree of dissoci- 
ation increases with -the Mach number. Further with the increase of 
equilibiium degree of dissociation the shock thickness decreases. Ihis 
fact is reflected in Elg. 7»4.4 which, shows the variation of a an*! 

1/i with M.J for nitrogen 
7 .. 
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Por the comparison purposes hardly any ejtperimental shock 
profile data for nitrogen in the Mach number range considered is 
available. Thus no definite conclusion can be made as to the validity 
of the theory. However the two formlations~quadrin!odal and triiviodal- 
given in Chapter 4 are capable of predicting the shock structure in a 
vibrationally relaxing diatomic gas. Out of these tvra, tho qixadritiiodal 
theory is better than the trimodal one especially for the Uach number 
range 15 to 20, Further the analysis presented in Chapter 5 includes the 
dissociation effects in the gas and hence is valid for higher Mach numbers 
upto 30, For > 30, the ionisation and photo-recombination of the 
atoms and photo-dissociation of the molecules become appreciable. These 
effects are not included in the present fo rraulation. 

As particular cases of the quadrimodal theory, for vibrationally 
relaxing diatomic gas given in Chapter 4> it is possible to recover, botli 
the trimodal theory proposed by Haight for rotationally relaxing diatomic 
gases (vf^=0) and the Mott Smith bimodal theory (w^>sw^«a0) for inert gases. 
Further to enable to carryout meaningful calculations for shook 
structure, the parameters b, a and s for axially elastic sphorooylinder 
models corresponding to real diatomic gas molecules should be available. 

In the present case -^he values indicated in the Table 2,1.1 were used. 

The slio ok structure analysis given in Chapter 5 iot 
dissociating gas may farther be extended for a mixture of atoms and 
molecules of the same igas. In such a case the upstream atomic number 
density is nonzero and hence the full bimodal ans at z for the distribution 
funotiott for atoms has to be retained. Corrospondiitigly, appi'opr’iate 
nianber of relevant moments of the two Boltzmann equations have to be 
considered in addition to the equations of Chapter 5 
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APPMDXX - A 


pOIilSIOU INTEGRALS 1 . 


33ie collision integrals in the moment equations are given by 


IijC« = 


- Y(**g»n)K dn dx dx, (A,1,l) 


where 


g,n = {(Cj^ + ijjj^ X - Ej V sin 

- (c + wxr-Bv sin <j) e) }, n 
^ ** %• 

K dn elemental area on the interaction surface and y(^) = ^ 
for X ■< 0, “ 1 fot- X > 0. 

Siam the dynamics of binary collision (Chapter 2)we have for the impulse 
I transmitted at a collision from molecule 1 to molecule 2 


I = - 


(A*1,3) 


where 


^ = 1 + cos 4 >)^Ce>fn)^+ Cb+Bj cos (A»1.4) 

2 X 

further the post collision variables are given by the following 


equations 


c' =C +--Lji;C' = 

*'1 -1 m ~ ' m " 


> I 

u), - to, + -— 
-1 *1 Ji 


(b + Sj cos Cej n) 


(a.1.5) 


(ii = a ^ (b + B cos 4) Ce n) 
- » 
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Tae collision integral given by Eqn. {A,1 .i) Is a 22 fold 

2 2 

Integral, Without making a particiiLar choice for f (c ^ 1/2 J. w or 

X 1 4^ 

2 

1/2 s( 2B) ) it is possible to carry out 11 integratiom analytically, 
She dependence of the integrand on e, e. , n , w and w, is limited to 
the following combination of them; 


e*ii , e^»n , exe,'n, «,e(a) .1, m,*e, ( w.exn 

-rw .^X ^ ^ 

and Uj • 

IlhuB the vectors oj and iji, may be expressed in terms of e j ®i 
as follows (see Jig. for facilitating the integrations in w 

t^^space. In such a case we have 


Si 

and 


exn 


exn 


m = m{e cos 0 + -r cos ijj + e x 

•* ” exn 


exn 


sin i))) sin 6} 


e,xn e^xn 

SI = .Oj (ei cos 01 * c-,^^ cos * Si » sin *i5sln ep 


( A , 1 * b , 


= oi sin 6 , sin 

dm = sin 8 das d6 diji = dm^ doi^^ d^ 


and 


= Wj sin diOj^ dSj^ '^’^1 


Further representing e and ^ spherical co-ordinates 


we get 


nxci 

e 5» n cos 8 + {— — oos E 2 

* " 


fgl 


sin e} sin 8 


« n cos 8^ + (s cos ^ sin Ej) si^ ^1 



APPEUDIX - A 


collision IKTEGEAIS 

ft u--- — XJ-'in— r -1 u_. IP . i J 

The collision integrals in the moment equations are given by 

Y(-g*n)K dn dr dt, (A,1.1) 

IJ IJi •r A 

where 

g,n « {(Cj + toj^ X V sin 

- (c+wxr-Sv sin i e) }♦ n 
K dn elemental area on the interaction suirfaoe and y(x) « 0 
for X ■< 0, * 1 for X > 0, 

Prom the dynamics of binary collision (Chapter 2)we have for the impplse 
I transmitted at a collision from molecule 1 to molecule 2 

g.n 

1 = ^ ' (a.1,3) 

C/m 

where 

C ® 1 + ^^Cb+B cos (j))^(exn)^+(l>+Bj cos (A*1.4) 

Purther the post collision variables are given by the following 
equations 



“l (b + Bj cos (e^ X n) 

(0 * 6) « (b + B cos ij)) C® ^ tv) 

Jp 


(A.1.5) 
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and 


de - sin 0 d$ de 


de^ = sin dB^ de^ 


CAa.7) 


where s , t; and n form an orthogonal set. Consequently we have the 
folloY/ing relationships useful for the integration 


and. 


e^n = cos B ; 


evgj-n = • sin B sin Bj cos e 


(o.exn = ta sin 6 cosijj sinp 


toj»ejXn = sin 9^ cos if>j sin Bj^ 


Ca.1.8) 


Using these relationships we can rewrite Eqns. (a. 1.2) and (A.1.4) 
as 

g.n = v»n + cos 4i) ^ Bj - (b + B cos (fi) u_(_ cos i}* 

X sin B - V (B^ sin cos Bj -- B sin (?) cos B) (a.1.9) 


? = 1 + 




mb 

IT, 


{(l+¥^ cos 1 ^ 2 )^ 3 in^P|+Cl+B cos (f))^ sin^ B) (A.I.IO) 


V * " c 


(A.1.11 ) 


Now the integrations with r^^psot to ? and e-j^ may be 
Carried out^ as the dependence on them is only throu^ K» !Phus 
we have 


/ K de de^ (27r)^ K 


(A. 1 . 12 ) 
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The collision integral given by Eqn. (a, 1.1) is a 22 fold 

2 2 

integral. Without making a particular choice for f (c , 1/2 J. w or 

X i 4, 

2 

1/2 s(2B) ) it is possible to carry out 11 integrations analytically. 
The dependence of the integrand on e, , n , w and is limited to 
the following combination of them: 


e^n , e ‘n , exe,»n, m.eCw.,'1, w,.,), m.exn 

~ ~ ..1 — ■>. ~ ~ ~ ..(I ~i. HI ~ 

and (Un e.xn . 

M# JL "V X ^ 

Thus the vectors ^ and uij may be expressed interms of a > ©i and 
as follows (sea Big. for facilitating the integrations in w and 

space. In such a case we have 


exn 


exn 


0 ) * ti){e cos $ + (-r COS \I) + 6 X 

•' " 'exn 


exn 


e,xn 
## 1 ** 


sin tf') sin 6} 




CA.l.oj 


wi = (0, {e. cos 61 + (1 cos ill, + e, X sin 4'i)sin 61} 

1 -^1 1 Isixsil -1 Tii^ ^ ^ 

* (u sin 0 , sin 0j 

dm = (0^ sin 0 dm d6 dtjt = d(ij_j^ dw^j d^ 

and 

2 

dmj^ * sin d0j 


Further representing e and ^ in spherical co-ordinates 


we get 


e s? n cos $ + . cos e + n {r— — sin c) sin p 

" pxei 


Sj « n cos pj + (s cos Ej + t sin aj) sin Pj 
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where 


B = B/b , = B^b 


and 


K - 4a^ [1+7 ^Ci+B cos <f.)5Ccos a3sin^ 6 + (l+B^ cos (^j)5C-cos $^) 


. J ^2/s -t ^ 2 _rf- 


sin^gj}+(^)^(l+R cos cos sin^B^dfcos 3) <S(-cos 3)3 


In liie 11117 egrand . n appeajcs tixarougli the texia y*n « Hepresenting n 
intenns r and i we get 


CA.l* 


V yxi V V>ci 

n = — cos a + {- — — cos 6 + — x .p. * . .Z — 5} sin 

vxi 


V vxi 


ce 


(A.t.U) 


vn = V cos a 


and 


dn - sin ct da d3 


Ih.e dependence of the integrand on 6 if any is only through {|} and 
we can write^ at this state 


27T 

/ {i|} d6 = 2-^ {f} 

0 

ihus we have the collision integral given by Eqn. (a,1.i) as* 


(A. 1.15} 


I .(I) = -(27r)^ / TjT y'{-g*n)f££^l K sin a da dr d7, (A,1,16) 


13) 


where 
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and 


de = sin 3 d3 de 
dgj^ = sin 3j dBj^ de^ 


(A, 1.7) 


where s , t and n form an oriihogonal set.: Consequently we have the 
following relationships useful for the integration 
e»n = cos 3 ; £i»n « - cos 3j^» 


exei«n « •» sin 3 sin 3j cos e (A. 1.8) 

w.exn * (0 sin 8 cosip sin3 
^ 

and 

Using these relationships we can rewrite Eqns. {a. 1.2) and (A,1.4) 
as 

g*n = *^14- 'f' 3j - (b + B cos Wj, cos li) 

X sin 3 V (Pij sin 4»j cos 3^ B sin cos 3) (A,1.9) 


? s 1 + ^ {(1+T?^ cos sin^pj+Cl+B cos sin^ 3) (A.i.io) 

where 

V * 0. Q (A.1 .11 ) 

« wt « 

•i 

Mow the Integrations with rfiSpeQ’^ lio e and may be 
carried out, as the dependence on them is only throu^ K. Thus 
we have .. 

♦ 

/ K de dEj = C2 it)^ K . • (A.1. 12) 
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and 


dtj^ = b dCj^ du'j_j_ dto^j^ sin 6^ dgj d'H’j d<^j (A,1.17) 


f C-^) {- 

P 3^2 2TikT^ 2TrkT^ SrkT,, ^ 2kT+ 


2TfkT- >kTv 
P P 


J, (j) 


1 


2kT, 


- ^2£2i} FCwj^) 
2kT 


(A.1.18) 


in nju-ch ) aatisfiee the fell owing condition 

*«3 

/ PCw^j) d4»JjJ » 1. 


(A.1.19) 


157 naing this (a. 1.19) two trivial integrations with respect to 
and inay be carried out straight away. 

£he integrations in the velocity space c , C-j W be carried 
out by using the standard transformations 


c = u. i + V - G 


^4.- 

tl 


1 •* 


^ 5? .+ T, - 
ti tz 




(A,1.20) 


£r “.1 i - I - 2 


such that 


(A.1.21 ) 


dCdC^ =dVdG 

^ iOBpeotion it can easily be seen that integrand is dependant 
on T only through the exponential tam, Thus perfiopilng the 
inte'grations in V spaee we get the eoUision integral given by 

(A,1,17) as 
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v/here 

B = B/b , = B^/b 

and 

K = 4a^ [1+2 ^(1+B cos (())6(cos 6)sin^ 3 + (1+^j^ cos ^j^)6(-cos Bj^) 

sin^B]^}+(^)^(l+¥ cos Cl+¥j^ cos 4)psin^B sin^3j5(cos 3) 6 (-cos $)] 

CAa.l3) 

In the integrand, n appears through the term v.n , Representing n 

inteims v and i we get 
** 

V vxi V vxi 

n = — cos a + {— H-I— cos 6 + ~ x S rt .. sin 6} sin a 
" V vxi V vxi 

• A* fv * ' «V ' 

(a.i.u) 

V‘n = V cos a 
and 

dn = sin a da d6 

A# 

The dependence of the integrand on 6 if any is only through {^} and 
we can write, at this state 

2v 

/ {^} d6 = 2ir {^} (A«1,15) 

0 

Thus we have the collision integral given by Eqn. (A.I.I) as. 

!• .($) « -(27r)^ / TiT y(“g.n3fj£^i K sinada dT dT*, (A,1.16) 


where 
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m 


I. .($) = —{— 

( 8 ^ 2)2 2 ^MTt,^Tt.) 




J- 


^1 

_> {.., .A .. } 


27rkT^ 2trkT_ 
H 3 


fW YC-g.n) exp {- 


mG 


■nkX^^ 


J. wf J-i w 


‘V. 


2 k CTt.-^Tt.) 
1 1 


1 “'U s( 2 B) s( 2 B 0 ‘ 


2 kT_ 2 kT^ 2 kT^ 2 kT- . 

»i ^3 


■) K sin a da dr’ dr. 


where 

dr * = dG {J3jjjdw^_ di|) dw^_L 

dTg = sin P sin p| d7 34> dy^ dp dP^ 
anil the ibll owing result is used 


{A,1.22) 


(a. 1.25) 




Tt ’^t 

2 T?k H ^ 3 ^ 3/2 


, f i 3 w 2 'i r"^. 


(A .1 , 24 ) 


'i ''3 J 

Rirther transforming ) and (^ix. s cartesian 

co-ordinates, after suitably non-^mensionalising tbcm we have 

m. dw, di^ - 


TT ■‘"y 


2k T^* 


(A. 1 . 25 ) 


w^L = —j^ 


^ dwj^y 


ndkere 




/■^ 


y y 2 k Tj.. 


u!_l sin 


(A b 26 ; 


“lx v3r ^ 


/ Ji 


= / — — (D^ , sin 4*, 

* rt*. rri 1 i- 
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cIt = bd c 01 d(o dw-- di|i sin 3 d3 d^ d <>) 

X 4h X* 


« b dCj^ ^“14. ‘^“Hi 


(A.1.17) 


« = !2_ (_JL_)V2 [_£l_) (J2„,V2 ,. 

P .„2 4i:kTt ^ '-2itl<T^ '' \kTv ^ 

rtTT T.p Tp Vp 


fflCC - Up i)‘ 

2kTt 


. . S(2B)2, p 


2kT^ 2kT, 


•} PConv) 


(A.1.18) 


In v/Jiich P(w^ ) satisfies the fell owing condition 

.CO 

J PCoiH') dW|j * 1. 


(a.1.19) 


By using this (A. 1.19) two trivial integrations with respect to “f 4 
and oij II mny be carried out straight away. 

The integrations in the velocity space c » laay be carried 
out by using the standard transformations 


s-^ — u a. 

"1 j - 


o. ir 

n 

^4.- 

tl 


*p V ^ 

fJ- 

•V 

T + 
ti 

5^4.- 

■fca 


n 



'TV T“ 

\T 

, 

“57. 


(A.1.20) 


such that 


d c d c = d V d a 


(A.1.21 ) 


inspection it can easily be seen that integrand is dependent 
on V only through the exponential temif Thvis perfopning the 
integrations in V spsce we get the collision integral given by 
Eqn. (A,1.17) as 
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It be easily seeb that the Ihtegraad depends on and ttough 

only the exponential terms wMoh may be integrated by using the following 
well" Imown integral 


-2 


»(a. 


/ e ^ dS'y = 


(A-1.27) 


finally we use the following transformation for w and w 

'X lx 


n\ T. 


n = {b+B cos 


Z S' (b+Bj^ cos 4ij) 


fi _ 


/2k T 


Sjf sin 3 - fb+3j cos 0^) 


■] 


j; ' “lx 


^2k T. 


u sin 0 

Ji 


1 - v ■Tp”ix 


(Aa.28) 


(Dt„ sin B 


such that 




where the Ibcobian of transformation is giyen by 


2k T. 


r. 

.1 


2k 


J S' {(b+B cos — j — ~ sin'^ 6 + (b+Bj cos — i- sin^ Bj) (A,!, 

As t£® dependence of the integrand on Z is only through an exponential 
term, the 2 integration can be performed using the standard result 

^2 

'O 

(a, 1.30) 

(2Bf 


/ e dZ = FS 


Sius we tore for any f (c^, ^/2 (^Z or ^ 

integral giyen by the 11 fold integral as follows 


29) 


) the collision 
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_ Ti • n • 

= -(2Tr)2 — i-l_ {, m 


(8^2)2 2rk(Tt.+T^.) 


,3/2 ^ 


J 1 

.} { ~ .. A ,., } 


2TrkTr 2nkT-. 
i J 


— 3^/^ C- ■ ) /T|T Y(-g.n) exp {- 

TrkTy. 


mG 


2k (T^-.+T^J 
1 D 


^ ^ s (2B) s(2B0 \ ~ • j j t j“ 

2kT 2kTT, 2kTy 2kT.. 6 

i n 

where 


(A. 1.22) 


dt' =■ dG ujj_dtD_^^ d(|> (0j4. dij)2 

dtg = sin B sin 3j dB dF^^ 6^ dy, dB dSj 


(a. 1.23) 


and the following result is used 


/ exp {- 


T +T 


Tt. T^. 


V^} dV « { 


T. T. 
2irk ^ 


m T* +T, 


i}3/2 


ti*t. 


(A.1.24) 


Further transforming ) and (wjj^ >'l'j ) to s cartesian 

co-ordinates , after suitably non-dimensionalising them we have 


2k 


u, d(0j_ di|) “ 


1 — 
— dll) , dio 
X y 


where 


2k Ty . __ __ 

"lu do>i.v di}i^ = - — ■■*■■ d(0j^ dwjy 


/“Tf _ /-JJ- 

Wv = / T , — cos ij/, u)„ » / — (I) X sin ’1' 

V 2kT^^ -v y y' 2k 

"lx V iTiT “1^“' *1’ “ly V *' 

j M 


(A.1.25) 


( A -1. 26; 
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r IS ,3/2 

" ^3/2 2TrkCTt.+Tt.) 


, 8s ^1/2 


f 8s \ l/2 v2 
3 


/ T|T YC-g.n) e3cp {- 


mG 


2kCTt.+Tt,) 


-n2 


s(2B)^ 
2k Ty. 


5C2Bi): 


2k T. 


} ^ sin a da d6 dn dt\ 


(A.1,31 ) 


V. 


where 


ri = n//T 


g,n = V cos a » n + X 
X =--v(3^ sin cos - D sin<fi cos P) 

®iere are yet 11 fold gt integrations which are to be carried out* 
!Phese can be done only after a particular choice is made for f » 

A.2,1 ^ * C ^ : 


(A.U32) 


Using the Eqns. (A. 1*5) we can write 

*x 


{c^} . 2 V 1 . 21^ i 

^ ^ ^ m2 


(a. 2.1 ) 


witere 


* x.n = 


cos 




sin 6 sin a 


V 


(A. 2.2) 


Eons. (A.I.II) and U.1.14) are used in writing down the Bqns. (A.2.2). 

As meiitioned the dependence on 3 is only throu^ 't ^ sJtd the 
integration with respect to 5 is carried out as follows (Use Eqn. (A. 1.15)) 
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It can be easily seen that the integrand depends on and through 

only the exponential terms which may be integrated by using the following 
well known integral 

-2 

(A.1-27) 


« '•Wy __ 

/ e dm = /a 


Ilnally we use the following transformation for and’ , 




n = (b+B cos 


— tUjj sin B - (b+Bj cos 


Z *> (b+Bj cos 


/2k T-. 

i/— !i 

V Jl 


Wjj sin - (b+B cos 


/2k T, 


- “lx ^1 


/ 2k T 


(A.1.28) 
lx ® 


such that 


-2 . -2 


b) 


. 2.,2 


+ J(t) +Z ) ; dm^ ~ ^ dn dZ 


where the Jocobian of transformation is given by 

2k 2k Tj., 

J = {(b+B cos ^ sin^ B + (b+Bj cos •■» j sin^ Bj) (A, 1 

As the dependence of the integrand on Z is only through an exponential 
term, the Z integration can be performed using the standard result 

» -Z^/J 

/ e dZ = /ffJ (a.1.30) 

mmOO 

2 

Thus we have for any f (c^, 1/2 J. wj or ) the collision 

integral given by the 11 fold integral as follows 


♦ 29) 



(a. 2.4) 


Substituting from Bans. (a, 2.1) and (A.2.2) in (A,2 o) we get 


21 21^ 3v?-v^ 

ill ,/ 


COS 2a} 


A.2.1 Integration with respect to n •■ 

Hb perform the integration with respect to rf ^ we substitute 
the Talue of X from Eqns. {A.1.5 ), (a. 1.9) and (a.I.II) in Eqn. (a, 2,4) 
and in turn substitute this value for in Eqn. (a, 1.31 ). fhis 

reduces the p integration to the following 


/ YCn-A}{(g^A^+g2A^)-(2Ag^+3A2g2^'^'^Cgi+3Ag2}ri^-g2b^}e"’^ dp (A.2.5) 


v^ere 


(A.2.6) 


A * V cos a + X 

Xhe above integral can also be written in the form 

__ -2 _ 

j 3/2£ { ^A+3g2A^)^ +(gj+3g2A)Ti^*g2n^}e ^ (A.2.7) 

A 


Si 


=•>2 


^ cos 0< 


3v?»v^ 


'^2 " 


25^ 


{- 


cos 2a} 


V 




and 


v&ere 
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ninj 


m 


I. . fft") a ^ ^ { - 

ir*-' ^S/2 2TTk(Tt.+Tt.) 

1 j 


f T|T Y(~g.n) Qxo {■ 


^3/2 , 8s ^1/2 , 8s U/2 ^2 


mG^ :r2 5(2B)^ 

i»ri ^ 

2k Tv. 


2kCTt.+Tt,) 


where 


s(2Bi) 

2k Tv7 

3 


2 


} F sin a da dG dn dT. 

m 5 


(A.1.31 ) 


n » T\/\/J 

g*n = V cos a n + X (A.U 32 ) 

X --v(Sj^ sin (|)j^ cos "• B sin^ cos 8) 

There are yet 11 fold gf integrations which are to be carried out. 

These can be done only after a particular choice is made for f • 

A.2.1 $ » C ^ : 


X 


Using the Eqns. (A,1*5) we can write 


{cb = 2 V I + 21^ ^ 

^ » ».2 

^ m 


(a.2.1 ) 


where 


n^ = i.n = — cos 

A y 




sin 6 sin a 


^ ""lx - S 


(A, 2. 2) . 


Eqns. (A.1,11) and (a. 1.14) are used in writing down the Eqns. (A.2.2). 

' f '2 1 

As mentioned the dependence on 5 is only through 'C > and the 

'X 

integration With res peot to d is carried gut as follows (Use Eqn. (A. 1.15)) 
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g = , 2? ^ ; g2 *= {(UZ) + m-l) C— 

^ V 

V 

A = a/j ; V = vt/J ; Z = Czr)^ J "n 

V 

The integration of (A.2.7) gives us a function of A, , g^jJ 

^ + 2g2^^) 

1 1 _ ■5 

- 2 ^ §2 ' 2^Sl ^ 2i A 


(A.2.8) 


+ g2 


erf 7^} 


(A.2,9) 


A^2*2 Integration with respect to a : 

We have from Sqns. (a.2,S) and (a, 2.8) 

A = V cos ct + X (A.2.10) 

This trail sformation is used to change the a integration to integration 
with respect to a . In such a ease we have 


sin a da = -dA/v 


V 2 = X -S' V and = X - v 


(A.2.11 ) 


vdiere and are the lower and upper limits (corresponding to 

as 0 and a * IT } for the new variable A, of integration. 

'Phjs means the following integral has to be evaluated 

3/2 _ 

/ F C A ; g,,g2) dA 

V 

It stovaa be Mted at this stage that i, and appealing in F , 
th^elves are funetlons of A, given by Bins. (A. 2 . 8 ). B>«s after 
substituting &r i, and we get the above integral as 
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2Tr {c^} = / {Cj} 

Substituting from Eqns* (A.2.1) and (A,2.2) in. (A,2.3) we get 


(A.2.5) 


{c2} 

X 


21 ''x ^ 21^ 

i- - c“ » ♦ -2- <— r 

m 


3 v?>v 2 


cos 2a} 


(A.2.4) 


A, 2.1 Integration with respect to n - 

To perform the integration with respect to n , we substitute 

the value of I from Eqns. (A,1,3), (a. 1.9) and (A.1.11) in Eqn, (a,2.4) 

2 

and in turn substitute this value for {c^ } in Eqn, (a. 1 .31 )• This 
reduces the rf integration to the following 


f YCn-A}{Cg5^A^+g2A^>-(2Agj^+3A2g2}n+(gi+3Ag2}ri^-g2n^>e"'^ dn (A. 2.5) 


where 


Si 


=.2 


^ cos (X 


v^+v^ 3v^»v^ 

1 j- X X - , 

>g2 * + - — “ — cos 2a} 


yp- 


v'- 


(A. 2. 6) 


and 

A * V cos a + X 

The above integral can also be written in the form 


-2 


J^/^£’{(gj_^+g2^)-(2giA+3g2A^)in +Cgj+3g2A)^-g2^)®*'’^ dn 
A 


(A. 2. 7) 


where 
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3/2 v, 5 -2 

£ " £ - '£ 


I ^ I I ^^ 0 * 0. erf A} ^ dA 


V £=1 


r3/2 


= -r- 0 Uv2;v, 5,C,50 - H(vj;v,C,C,X)}. 


where 


/n _ 3 f3£-l -2 


„ _ 2a 3x( 5ai) _ _ 3(3ai) 2K ^ i-e ^ c3C-nx 

„ j r T ” r—'TT' * t: — ■*■ '"^T" ' ‘ IT"' 


72 ?2 


2 2 y2 ? ^2 


2 


2C3?*1)X . ^ ^ 5? - 1 


^2 y2 


;2 72 


Qo 


1-g ^ (Sg-.!))^ . Q _ g X X(5ai) 

9 — o 9 ■* '1 r 9 —'7 

2c Sv'^ ? ^ C" V 


^ 3S - 1 £ , ^ C3ai)X 


2^2 y2 C 


25^ 25^ ^ 


^3 


5g.l 

2 -2 


0, = ^ 1 ^ 

2?2 72 


Qc = 0 


K = 

o 5 


^-3^ . 3(30]^ 3(1-5) L 


2,2 v2 


2 , 


2 5 




2CX 3(35-l)X 


72 


g . SiK-D 25 ^ J 

* ' 2^2 72 ■ ? c2 ?V 


^2 


K,.. 2 W-m R . 35-1 

,2 ;t 2 3 ^ 2:72 


a v‘ 


5 ^v‘ 


(A. 2 . 12 ) 


CA. 2 a 3 ) 
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gj = , 2? ^ ; g2 » ^ {(1-0 + (301D C^)^} 

V „ (Ai 2 .8 ) 

A = a/j ; 7 » v/j ; C « C::—)*^ ; n == n/^ 

V 

The Integration of (A.2.7) gives us a function of A, 

P (A,g]^,g2)<.T.- « {^ (gj + 3g2 A + 2gj^ + 2g2A^) 

1 ^ 1 
- J (g2 Si ^ g2 ~ l(8l + 3g2 F + g^ 

+ g2 ®rf X} (a. 2. 9) 

A,2»2 Integration with respect to o : 

V/e have from Eqns, (A,2.6) and (A,2.8) 

A ss V cos a + X (A.24IO) 


This transformation is used to change the a integration to integration 
with respect to A . In such a case we have 


sin a da =5 -dA/v 

V 2 = X + 7 and Vj^ = T - V 


{a. 2.11 ) 


vfhere Vg and are the lower and upper limits (corresponding to 

a =5 0 and a * ir ) for the new variable A, of integration. 

This means the following integral has to be evaluated 
,3/2 

* / P (A ; gi»g2) dA 


It shomd be n»ted at this stage that i, and ^ appeariag in F , 
themselves are functions of A, given by Eqns. (A. 2. 8), Thus after 
substituting for and we get the above Integral as 
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and 




S 

I < 

t-l S,+ l 


P. 

. X. p 




(Vp) ) 


(A.2.14) 


Ihe iategrals 

\(x) 


X ^ 

/ e"^ dx and (h)^{x} 


/^ x” erf X dx 


(A.2.15> 


appearing in {A.2.I4) have been listed in Appendix B 


Integrations in the G space ; 


Bqns* (a. 1*11) and (A,1,20) give the Allowing transformation for G 


V = G - 


U, , X 

13 


ifoere 


.= u. ~ u 

i 


Id 


(A.2.16) 


{A.2.17) 


Bepresenting G in spherical co-ordinates we have for the elemental 

w 

voligae in G space as 


dG 


G G* sin Yj^ dVj dG 
Bro® Bq^n, (A,2*16) v^e can get the following results easily 


(A.2.18) 


* G ^ - 2G U- . cos 7i ’*■ 
T ^ 5 cos Y “ ^ 

X 3,3 


(a.2.19) 



60 . 


3/2 V 2 5 

~I II (^0 - Q, - R. erf A} dA 

V Vj It-l “ ^ ‘ 


j3/2 

V 

v/here 


n ^ 

'’o = r- 




.35. 




1 ^ 


+ 1-5} 


25X 3X(55-1) p ^ 3C3€-1) 25 , lyg ^ (35^1) 

5 “ ^ ^2 ’ 3 - 2 2 ^ " 5 ^2 ^ ^2 


(A. 2 . 12 ) 



2(35-l}X 


35 ~ 1 
^2 ^ 



lzl+ . o 

2^2 2^ 


5 X X(35.1) 

5 ^2 y2 


Q, = 35^ 
^ 2 c 2 72 


i + M + (3g-l)X^ 
^ 2^2 2^2 ^ 



35-1 

2 —2 
C V'^ 


Q4 


35 - 1 
2^2 ^ 


Qs " 0 




3(3g-l)l2 ^ 3(1,5) , i 
2^2 2^2 ^ 


CA..2.a3) 


= , 3(;35-l)X 


c2 


R - 3(35-1) 35-1 ^ 


2^2 ^ 


C2^ 
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first of these two, gives a transformation, for integrating with 
respect to Yj 


sin YidYjS= 


GO,. 


(A.2.20) 


where G » G/ /J and ^ = u / /j 

-i- J •A.J 

further ttie (A*2.8) and (A. 2. 19) give e 


function of G , 

and 

u. . 

ID 



G 

- 2 2 

T - U. 7 

2 


€= ( — - ) 
2 V u . . 


(A,2,21 ) 

as the following 


(A.2.22) 


^■us -aie volxime ele&eat d G given by Eqn. (A.2.18) reduce to 


dG » G — — dG d'f'j^ (A.2»25) 

' n. . 

and the region of integration is 
i 0 to 2v 

V iCg-uT-) to{G + u. . ) 

13 *•’» 

G { 0 to « 

!Hie integration with respect to gives a factor of 2ir an 

ragion of integration (OBOT) in the 5, opeee is aho« in ?1«. A.2.1. 

25 iue we Jhav^ evuiuste foXlowjJ^g iiitogrsl 



CHj - HjJG 




dv 


(A42.24) 
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and 


rh 


I'p=HCv ,v,5.C,X) = I {: _L-. , 0^ A (v ) - Rp®, (vj) (A,2.u) 

Jl=l A+1 p p 


The integrals 

X 2 

Ar.(x) = / x" dx and (k)„(x) = P x” erf x dx (A.2.15) 

0 0 

appearing in (A.2.14) have been listed in Appendix B 
Integ 3 ?ations in the G space : 

Eqns, (A.I.II) and (A.i,20) give the following transformation for ff 


V = G - 




where 


Uj r u. 


id 


d 


(a,2.16) 


(A.2.17) 


Hepresenting 6 in spherical co-ordinates we have for the elemental 
volume in G space as 


dG « sin d\l?^ dG 


(A.2.t8) 


From Eqn, (A, 2.16) we can get the follovdng results easily 


v^ » G 2 - 2G u. , cos Yi + u. ? 

id ^ IJ 


(A. 2. 19) 
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a = n J/2k(T^ + ) 

i 3 

“ ^‘2 = HCv 2 ; v^£,c,X) 

(a. 2. 25) 

Hi = Hi C7 ,G;c,x) = HCvi;7,5,C,X) 

It can be seen from Sqa. (a.2.11), {A. 2.I3), (A,2.22) and (a. 2.25) that 

H2£-7,-G;C,X) = H2C-v,G;c,X) = KiC7,-G;c,X) = Hi(v,G;C,X) (A.2.26) 

fhe integral (A. 2. 24) is eqxd, valent to the difference of the contributions 
due to the region AfE (Sector of circle with infinite radius and centre 
Dj Ixjunded by the lines G = 0 and = G — iT . ) and the region 
ABC (sector of circle with infinite radius and centre bounded by 
lines G = 0 arod tT = G + ). !5bis is possible only when the integrand 

is well behaved in the domain ADS (which contains ABO as a sub domain 
in it), Such that the integral with « limit is not an improper 
integral, %■ careful study of the integrand in (A 2.24), that is, 
the functions given by Bkjn. (A. 2.1 3), it can be seen that the Integrand 
is bounded due to the exponential damping terms and hence it is well 
behaved, in the domain ADS. 

Bans we can write the integral (A. 2. 24) as 

S^iT- . “7^2 G+u. . ^ 

r ^ f dv - dG f (H 2 -Hi)Go' dv 

•Wj 0 “11°. 
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The firet of these two, gives a transformation for integrating with 
respect to 

sin YidY-i» . (A.2.20) 

11 Ouij 

where G » Q/ /j and u!". =:u..//j (a,2.21 ) 

Xj 

Further the Eqns» (A.2»8) and (A, 2. 19) give g as the following 
function of G , and uT. 

G V u. ? 2 

-JO-) (A.2,22) 

Th\is the -volume element d G given by Eqn. (A,2.18) reduce to 

dG » G -- — " dG d<f<, (A.2.23) 

U. . 

iD 

and the region of integration is 
’I'l : 0 to 2n 

V : ( G - ) to ( G + \^^ ) 

G { 0 to « . 

The integration with respect to gives a factor of Zn and -the 

region of integration (OBDB) in the G, space is showti in Slg« A. 2.1. 
Thus we have to evaluate the following integrali 

/ ^^^2 ” ^1^^’ ® (A^2.24) 

0 G-u. , 

, -ID 
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M-D&lly f by effecting the chajige of order of integrations and using the 
syHHjetrtc properties (A. 2. 26) we get the integral 


r d7 H2 d6 

-CO v-u. . 


(A. 2. 27) 


where 


H2=HCv2;v,e,c:>x) = 


I 

£=1 £+1 


- Qji ^£^^2) - ^£0/^2) > 


are functions of v , C , 5 ^ y and 5 is a function of &, v , 

given respectively by Sjns, (a. 2,15) and (A.2.22). Ihe integration with 

_ 

respect G- can be carried out, as is an integrable function ^G, 

throng tfie functions which are in turn dependent on g * 

However the analytiGal integration, with respect to v has been not 
possible. iSius we get an integral with respect to v as 


/ H (v; a,C,X) dv 


(A. 2. 28) 


where 

c TT 

H2(v,a,5,X) = I ^ ^ 


£=1 £+1 






where 


(A.2.29) 


^ Pj G e dG 

2 

%(v,a,?,X) =7^2 G dG 

_-2_ 

== Rj G dG 

V2 - T + V 


(A. 2 . 30 ) 


Gi Y and V + ^ 



&3 


where 

a = m J/2k(T^ + ) 

i j 

Kj « ^2 (7 .G;c,X) = H(v 2 ; V»g,C,X) 

Hi = Hi Cv,G;c,^ * H(Vijv,?,s,X) 

It oan be seen from Eqn. (A,2»11 ), (A.2.I5)) (A,2.22) and (A. 2, 25) that 


H2(-v,-G;c/X) » H2(-v,G;c,X) = Hi(v,-G;?,X) = Hi(v,G;?,X) 


(A. 2. 26 ) 


The integral (A. 2. 24) is equivalent to the difference of the contributions 
due to the region ADE (Sector of circle with infinite radius and centre 
D, bounded by the lines G c 0 and V « S'— u^^ ) and the region 

ABO (sector of circle with infinite radius and centre B, bounded by 
lines G s 0 and v = G + uT . ). This is possible only when the integrand 

Xj] 

is well behaved in the domain ADE (which contains ABC as a sub domain 
in it), such that the integral with « Ircait is not an improper 
integral, ISy careful study of the integrand in (A 2.24), that is, 
the functions given by Eqn. (A. 2 .I 5 ), it can be seen that the integrand 
is bounded due to the exponential damping terms and hence it is well 
behaved, in the domain ADE. 

Thus we can write the integral (A.2.24) as 

,« _ ^ 2+u* . 

/_ dG / •’ (H2-Hi)Ge" dv - £ dG / CH2-Hi)Ge*“'’ iv 

0 u- . 0 


— 2 
wttG 


11 


IJ 
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Bie integrals (A. 2 . 50 ) are listed below 


p = 


; Px = 

{hp ^ r } 

0 

2 

^ i 

2,2 ^ 


3Ci X 

3 X £3 _ 

(3+2X)53 

^2 = 

? 

—2 2 ’ ^3 

2,^v2 ^ 


-25 jX 


^3 

P4 = 

2 —2 

5 v^ 

; P5 = 

52 72 


_ ?2 C, _ £l X X £3 

^ ■ 25^ 272;2 ' 1 £ ?l2 




„ !i. 

Qa = 

2t} ^ 

c 


,,' 3 , 


Q 4 == 


> 

% = 

SjX 

c 

; 


2^1 X 


R2 = 


^ ^3 

tl . 

-2^3 X 

; Rg 

^4 = 

2 —2 
, V 


£2 , - _ ~^3 

2;X ^ 

Q 5-0 

3^3 2^2 ^1 

2?2 ^ ^ 

(3+2X^)?3 2C^ e2 


C2 ^ 


CA.2,31) 


w^iera 

6 , . 5 Ge-“®*dG= £j 2 - £ii 

52 = C1-£)G e'“®^ dG = (£02 - Epi) - CEl2 ' 'iP 

^1 
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finally, by effecting the change of order of integrations and using the 
syiometrio properties (a. 2.26) we get the integral 



v+uj . -s*2 

dv ]_ ^ H2 Ge"®^ dG 

V-Ui- 


(A.2.27) 


where 


H2=H(v2JV,5,C,X) 


£=1 




U1 



Jj, are functions of v , C , C , X and t is a function of G, v , 

given respectively by Eqns. (a. 2.13) and (A.2.22)t The integration with 

- - 

respect G can be carried out, as is an integrable function ^G, 

through the functions which are in turn dependent on £ . 

However the analytical integration, with respect to v has been not 
possible. Ihua we get an integral wtlii respect to 7 as 

CO 

/ H* (vj a,c,x) dv (A. 2. 28) 

^00 

where 

i; rT ^+1 

^ ^11 ^^2 

H2Cv,a,5,)0 » I — - A^(v2) - ®, (^ 2 ) (A.2.29) 

where 

Q - . 2 

= / ^ p, G e'“® a5 

5.Cv.o,c,X) = /°2 Q G e'“° dG 

=1 U.2.30) 

— . G'? ^ r ^ ^ 

Rs,(v,a,?,X) = / R. G e““^ dG 

■Oj 

V2 = T + v* 

Gj = V - and G2 = V + u^j . 
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5j = r {35-i)G K„2- w 


- 2(72 * s?.) ^Aj (G^ * i-AjCG^ (A.2. 

op 


32 ' 


5op = Z h (Gp and (? = -r G 


Shios T?7e havey -fee collision integral for ? ^ c as a seven fold 
ijitegr^l given by 




-2i1j11. 
^ J 

A 


in 


TT^ 2k(t^ +T ) 


,3/2 , 8s , 1/2, 8s , 1/2 r ,3 
kT kT ^ J 

Kiy, ^ 


(v;a,?,X} exp {- }K dv dx 


2kT, 


21cT, 


ij 

6 


J 


(a,2,35) 


2 

Ihie expression I. , (c ) teikes o/o form when i = j. However it is 

i J ^ 

more coavinient to carry out the integration in G space, separately 
for this case ( i = j) » 


iioitt Spa, (a,2,16)j (A.2.ig) and (A,2,22) we have for i « j 


u. . *= 0 ; V = G ; V = G cos Yi 

ij -V A ^ 



2 

cos Yj, 


(A.2,35) 


V 

Bius in E(in, (A.2.12) the functions dependent on K 

c^ be integrated with respect to Y p to obtain 


Pjl CvX?) ' r Y| dYj 

Cv*X,5) = f Qj^ sin Yi ^Yl 

*^1 (v,X,f) R. sin Yi dY^ 

0 


(A.2.34) 
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Kie integrals (A. 2 . 30 ) are listed below 


P- = 


/fT ^ 


o 2 




Hi X 3X ^3 „ C3+2X)5- 2e, 

p = . ...... ^ , n _ ^ 

" C -2 ,2 ’ ^3 - “ 


V c' 


2c 


-2e3X 


^4 ® 2—2 


J Pp K 


5 7 

V 


Qo“— ^ 


?! X X 

2 ’ ^1 = — ' 


2c^ 2v^C 
CjCl+l?) Cj ^2 


Q_ s ' " '- + ■■■ ■■^■' j O 7 ^ 

A ^ o ^ r ^ 7 ' 3 




252 v 2 


2C' 


i;2 ^ 


(A.2,31) 


Q 4 




R 2 * 


«4“ 


^3 


^ 5 " “ 

2cV 


5jX 


_ 3£j x2 ^ 2£2 

c 

; Ri 

” 2^2 72 ■" 2?2 

2Cl X 

; R 3 

(3+2X^)€3 


2^2 ^ 

- 213 ? 

» % 

^3 

^272 

C 2 ^ 


C ^2 


where 


r* *“** Z 

-2 . ^ -aG 


h ‘ I 5 '^‘= " hi • hi 


G ***"2 

§2 “ / Cl-e)G e dG * (5^2 *' ^61^ " ^^12 " 

.' ®1 , \ 
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so tbat liie collision integral be written as 


o 2 

2. 

T. , fc 7 = 

IX '■ 4 




1 

'~2 


V. 


______ s(2E]2 + s(2B.)2 _ „ _ 

H2Cv,a,CjX) V ejq) {-a — ) K dv dt, 


2k T 


V. 

1 


(a.2.35) 


where 


^ ^ I ^1^2 - =/-\ 

H2 Cv,«,c,X) J “ Qjt 


(A. 2. 36) 


1^=1 £+1 

are the functions ij_ given by (A.2.31 ) 

in wkLeh replaced by and respectively 


isiieir0 


= /^ e sin dYj = 2/3 


0 


^2 = r ip 


0 


(a.2.37) 


€3 “ /^ (35-1) sin Yp dYp = 0 . 


0 


^Ilius the collision integral for f ^ reduced to a seven fold 

integral both for the cases i ^ 3 i » D* 

A.3 1/2 fe>4-^ 


s, _i^ (j,)2 (i,(e*n)^) +(b+Bj^ cos 4’i) 

■'IT ^ «• w 


2J 


(1 - (en,n)^}} + I (Cb+Bj cos (^p) wp 


exn 


(A 3 fl 


. ( 1 ) + B cos ^)w . Cp’^^ 
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?3 = / (3C-1)G dG * 3(5^2 “ hO “ (^o2 - hO 


hv = -zr-rr + S2 ) (G') + ^AiCGM (A.2.32^ 

-^p atTI axp 


4^1? . v2 53 ' P 

=! — A, (C) and G* = *'*a^ G 

^op 1 '> p-* p 


a 


IHius we have, the collision integral for | = as a seven fold 
integral given by 


ID X 4 


m ^3/2 ^ 8s -,1/2, 8s ^1/2 b ( ^3 

^15-’ 


'J 


1^2 (v;a,C,50 exp {- |i|^ - 


13 

6 


{a.2.33) 


aiiis expression I. . (c ) takes o/o form when i = j. However it is 

ij X 

more convinient to carry out the integration in G space, separately 
for this case ( i = j) * 

Prom Eclu, (^.2*16), (A, 2. 19) and (Ai2.22) we have for i » j 


u. . <» 0 ; V = G ; V « G cos Yi 

13 ^ ^ X J 

/x.2 2 

5 = (— ) = cos Yj. 

V 


(A.2.33) 


Thus in Eqn. (A.2.12) the functions P. jQ-t®. vdiich are dependent on 5 

Xt jC 

can be integrated with respect to Yj to obtain ‘ 

Pjl Yj dYj 


(v>X,?) = /^ sin Yi ^Yi 


(A.2.34) 


(v,X,C) = f sin Yj dYi_ 
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uslJig the transformation for e j n , w and we can rewrite 

this as 

2 

{f} = <^j^)^sin^ej+Cb+B cos d)^sin^B] - In (A. 3. 2) 

wbare I,n,^,p 3_are given by Egns. (A.I.3 ), (a,1«8), (a. 1.9), (A.I.IO), 
(A, 1,28) and {A.I.32). In this case {|} is independent of 5 , so 

2^ - — 

/ {{} AS - 2^W) = li 'i 

0 

Int egration with respeot to rT 

Substituting this C?! ia (a. 1.31) we find that the Allowing 

integral in n appears in I )• 

— 2 _ . 

/ vCn-'A) {pj^CA-n)^ * ^ 2 dn (a. 3. 3 ) 


where 


hi = 


{(b+B^ cos sin^ + Cb + B cos 4)^ sin^ 6 > 


2^2 


>^2 = 


B 


and 


(A.5.4) 


A is given by Eqn- (A. 2. 6) 

Saie Integra (A.^53) can be re^itten vdth the proper limits for n 
and aft^ carryliig out the n integration we get the 

A ^ ^ " * '^“1 ■ "" 2 ) ^ " '*^2 ' “ i ’ “ 

A 

» j5/2 [(3Vj-2v2)A+2yjA^] j - Pi A ] 
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so that the collision integral can be written as 


0 2 
*2n-i 


X rJ-s _ ' i . ra ,3/2 , 8s ^ ^2 f ,5/2 

’ii -4- ’ 'isr’ ' 


^ _ s( 2B)2 + s(2B,)^ _ _ „ 

K2Cv,a,C,X) V exp {-a ^ - — * — ) K dv dx^ 


2k T 


V* 

1 


(A.2.35) 


where 


H2 Cv,a,C,X) I 


c F ^+1 

5 ^^^2 


(Ji-s-se) 


l-l A+l 

are the functions Pjj, Qg) Bjj, given by (A.2,3l ) 
in which gj , ^2 ?3 replaced by and respectively 


where 


= /^ 5 sin y. dy, * 2/3 
0 

^2 * r Cl-C) sin Yj dYi = 4/3 


X3 ® r C3C-1) sin Yx dYj * 0 . 
0 


(a.2.37) 


Thus the collision integral for | = is reduced to a seven fold 

integral both for the cases i / 3 and i » j . 

A.5 1=1/2 J-, 


{$} = {(b+B cos (|>)^ Cl-»(e.Ti)^) +(b+B, cos 

^ ^ J. ^ 

2 

(1 - (cj^.n) )} + I {(b+Bj cos 4>i) . exn 


(A 3- r") 


- (b + B cos ({>)u > e^xn} 





ee 


* - j [(% - 2u 2)^ erf A) 


* ^ ^ (A ; pj, y^} . 

lategFatloa wi1it respect to a 

'j&e sasi0 traaafora^tion used for- iiie earlier case, given by Eqti, (A. 2.1 
is adopted :for the ct integration, to obtain the integral as 

5/2 V2 _ _ 3/2 

-TT” / — :r" {S(v2;yi,y2)“®Cviivi. ‘ 

¥ V. V 


wi^r® 


i/F 




16 




2 ^ 2 ) - 4y^ vg] erf [C^ya * 

^1 3, -P 

- — vie 
8 P-* 


la^' 


in G space : 


Hers the integrations are done by adopting the transformations 

used in the earlier case given by Eqns, (a. 2. 19) and (A. 2. 23)* ^hus 
we l^ve find the integral 


7/2 


G+Ui 


G e dv 


(A.3.8) 


r r I CS2 - %) 

^ij 

Bo^ ttese integrations can be done analytically , wi-th out eff S 
the of order of integra-tions, unlihe in t/he case 

Inters ting with respect -to v we get from (a.3«8) 


5yi)Vp 
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Further using the transformation for e , n , oj and we can rewrite 
this as 



(A. 5. 2) 


v/here I, n , 3 , are given by Eq.ns. (A,1.3 ), (A.US), (a. 1*9), (A.I.IO), 
(A. 1.28) and (A. I.32). In this case is independent of 6 , so 

2Tr 

/ d6 ® 2 ti TfT j 


Integration with respect to n 

Substituting this {fT in Eqn. (A.I.31) we find that the following 
integral in n appears in I. .( ^ ). 

2 

/ vCh-A) {yj^(A-ri)^ + h (A-n)^} dn" {A.3,3) 

where 

2 

hj “ — f {(b+Bj^ cos sin + (b + B cos (|>)^ sin^ 3} 

2c 

V2 “ 7 ^>=1^ (a.3*4) 

A is given by Egn. (A. 2, 6) 

The integral (4,3.3) can be rewritten with the . proper liniits for n 
and after carrying out the n integration we get the form 

—2 

Z ^ ^ ~ ^^ 2 ^ ^ ■*■ (^2 ■* 

= jV2 [(3pj.2p2)A.2Mjit3] * 1 - vj 5^] 
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jn 


y* . 
3-3 


/ Cs, 







dG 


(a.3.9) 


= Crib " 3uj s^i 

4 ^ |1S ~ 3i‘j)G^ - 20 (3^3^ - 2ii2)5'^ - 12y G^] (A.3.10) 

^ ** ^ " 54yi)G^- i2y G^] e“^P] > 

" P * P 

« G 4- X^ ^ ^2 ~ ^ ^ ^2 

CA.3,n) 

2s = Xj; 5-^ = -5-. Xj. 


Int^gSESffeleii with rfcspect to 9- is carri«i out using the transfoimations 
(4*5,1l} ^after rearranging the teisjs we get from (A. 5 . 9 ) the required 



i-. 


J/2 


u. 


{SCX^,*a,C} - ? (X. ; a,0) 


13 


(a«3.32) 





- %P ! «,?) = % [5fe-33lX4(«.Xp) 

c 

- \ XjCo.X )) + 3fe-l) Ix^CS’.Xp) - Xp xCtXp)}] 

* 2M C'>2f“'Xp) - Xp ^ga.Xp)} 

♦ 20CC- 3 ) {Xp <^3(0, Xp) - i>4(“'Xp)) + 


12ft-l) {X„ ♦., Ca,X^) * 240 ' 
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^ 1 

* ® ^ 2u2)/^ +■ Uj erf A) 


“ y 1 (A- ; pj, )i2^ • 

Integration vath reaioecl; to a 

i'he same transformation used for' the earlier oasej given by Eqn. (A, 2*11) 
is adopted for the o integration) to obtain the integral as 


V2 _ .3/2 

V Vj V 


{SCv2;Mi,y2)"S(vj;uj,P2)> ' (A-5-6) 


where 


+ ^ [4^2 - 3yp - 4(3pj - 2 y 2 ) - dy^ vg] erf Vp + ~ [(Ayg - 5yj)v^ 

V, , -vf (A«3»7) 


Integration in G- space : 


3 

r''pi ^ ’ 


Here again the integrations are done by adopting the transformations 
used in the earlier case given by Eqnst (a. 2,18) and (A*2.23)* Thus 
we have tfo find the integral 


r a. 

“ij 0 G-n-jj 


~ -aG'^ _ 


(S2 - Sj) G e " dv 


(A. 3 . 8 ) 


Both these integrations can be done analytioally, wi1ii out effecting 

9 

the change of order of integrations, unlilce in 'the case of ? = C^' - 
Integrating with respect to v we get from (a.3«8) 



{ 2 j(a,Xp) - - ( 27 - 7 ?) 

+ &(C-l){X^ 2 ^Ca,X^)-Z 2 (a-X )}]> 

■* -’ P 
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ae iBw^igeale 


]^C«#x) K / g” dG 


'^if 






{A, 3 . 14 ) 


€U2d 


XCa.x) 





erf 6 


dG 


Bsf0 eral'oateS and given in Appendix B, 

tefore arriving at (a.5»12)j the quantities and P2 appearing 
in (Aj»J*tO} are esgvreased in teisis of C using Bqns. (A.I.10) and (A. 3 * 4 )* 
Sslsf res^ilt CA.3,12) in Bon. (-^. 1 . 31 ) we get the collision integral 
toT ‘^xlB case as 


4 , 2. 


*13 “ 


T3 1^/2 ^ 85 j.1/2 ^ 8s 1 ^1/2 b 

.4 ' '■2k(T..-Ht4 kTy. ^ 

1 J 


‘ij 


/ ($, - S4 c]a' { - 


.oa'>2 _ 

’ _ ■} dr^ 


2kT, 




2k T 


(A.3.15 


V. 


a six fold integral which has to he evaluated numerically, 
ia ttj® c ^ ca®0y ?&en x - h reduces to a o/O form. !ro 

.tfesoQiiisloia integrals for such a case we follow^ 


as foil owed for the case. 
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— / (S^ - s 

Ui. 0 


2 + - S4) G e 


<iG 


(a.5.9) 


where 


=p = : vi-va) - (scsui - 2 p2)g 3 * spj s^j 

+ ^ [15 (4y2 - 3yi)Gp - 20 (Sy^ - 2y2)G^ - 12y G^] (A.3.10) 

HZ 

erf Gp + ^ [(40y2 - 24yj) + (40y2 - 54yi)Gp^ - 12yj G^] 

Gi = G + Xi ; G2 = G + X2 

(A.3ai) 

Gj = - G + X2; 64= - G + Xj, 


Integration vdth respect to G is carried out using the transformations 
(A«3.11) and after rearranging -toe terms we get from (A, 3. 9) the required 
integral as 

ii— {S(X, ;a,c) - S (X2 ; a,0> (A.3a2) 

where 


Sp = s(Xp : cc) = t5(;-3){x4(:r,Xp) 

" X3(h,Xp)} + 3(s-l) {x^Ca.Xp) - Xp xC«,Xp)}] 

+ ^ [15 C?+3) {i|)2Ca,Xp} - Xp i(/i(a,Xp)} 


+ 20 (?- 3) {Xp (ct, Xp) - ^4(a,Xp) } + 
i2(c-i) «p <.5 (?,Xp) - 4,5cs-.y)i * ^ [ 



'-Cl j? 
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of , the integration with 

resp^t to yi gtTG a £&sv^t 


r 


sin Y, dv^ = 2. 


te ixiilisic-n integral 




* ^ 3 

Z1 ( > .2 


(^) r sCv3,n,,U3) 


V {«iv2 - 


sC2n-3" - 5(220' 


?M 


-} K dv dT 




6 


^eh iM ohyiously a seven fold integral. It has not been possible to 
cari^i^ toe v integration analytically for 1 = j case. 

A.4 I * 1/2 : 

,{A^ t #5 j obtain for Ihis case 


{f} = - I V (Ij sin cos 6j - B sin ^ cos B) 
toto is lfidep®£^&£it of 5 in this case we get 

TfT - {.',) 

with gesveot -o n : 


(A.4.1 ) 


for I frcE Iton, (A.l.8),(^.1.9)j (A.1.10)sCA. 1.28) and 
(a*! *52) gst toe follo’sing integral in n 


J I V (A-y)*' e ■ dn 


H 


¥ 


^ fl? sin t cos B - B-, sin ^ cos B) 
^ ■* 


(A.4.2) 


f ~ X 


(A.4.3) 
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4(2?+3) {Zj(a,Xp) - Xp ZgCa^Xp) - (27-70 

(Xp ■" ^^^■^5{XpZ4(a,Xp)-Z5(a-Xp)}]} 


The integrals 

Xj^(a>x) = / g” dG 

-CO 

2 9 (A«3*14-) 

(c.x) = r a" e-= dG 

-lOO 

and 

ZCa,x) = /'“ g” erf G dG 

have been evaluated and given in Appendix B. 

Further before arriving at (A,3,12), the quantities wj and M2 appearing 
in (a, 3. 10) are expressed in teims of C using Eqnst (A.i.io) and (A.3»4). 
Using the result (A,3.12) in Eqn. (A, 1,31 ) we get the collision integral 
for this case as 


ro 


hj <2 •’l ’ ‘2«Tt +T ) 


j3/2 {_8s_jl/2 ^ 8s jl/2 b 


kX 


kX 


13 


/ (Sj - S2) exp 




2kX 


dt. 


2k T. 


(a.3.15) 


V. 

3 


This is a six fold integral which has to be evaluated nvunerically, 
2 

As in the c„ case, when i = j, this reduces to a o/d form, 3k) 
obtain the c ollision integrals for such a case we follow^ the same 


procedure as followed for the case. 
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Integrating (A.4.2) we get 


- r/F 


)= V C 1 + 2 A^) 


- ri.^T2 


2 (1+2A'^) erf A - 


Integration with respect to a : 

Tbl3 is done bjr using new variable A, as indicated earliecr 
A integration we get 


V 


Jv 

V 


/ ^ I. 2CA^.<iA = ^ " V(Vj}} 


V 


wfeera 


V = 


V(v ) = {•— C3v + 2v^} - r (1 + ® ^ 

t> 12 p P ^ P 


i (3Vp * zvp erf v^) 


Integration in 6 space ; 

Here ag at n the procedure used in the earlier case is adopted 
the required integral as 


,5/2 — 00 G+Uj ^ „ 

^ ^ (V 2 - vp G e dv 

ti-' ' 0 G-u - - 

2.3 13 

Integration with respect v gives us the above integral as 


ff.I r cv, - V2 . Vj - V4)e e-“« da 

u. . n 

13 


} {A. 4 . 4 ) 


After the 


(A.4*5) 


(A.4.6) 


to obtain 


(1*4.7) 


(A.4.8) 


where 
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As the integrand^s independent of > "tiie integration with 
respect to give a factor 

fit 

j sin Yi dYi = 2. 

0 

Thus we get the collision integral 


.4n? 


hi'^1 = -4^ CsTT’ 

ir '^x ““ 

sC2C3^ + s(2D^3^ 

V exp {-av^ K dv dr. 


which is obviously a seven fold integral. It has not been possible to 
carryout the v integration analytically for i » j ease. 


A.4 5 = 1/2 3 (2B)^ : 

Using! Eqns* (A.1 ,5 ) we can obtain for this case 


{$} = - I V (B^ sin cfij^ cos ■. B sin cos B) (a,4.1 ) 

As this is independent of in this case we get 

W ={0 

Integration with respect to n ; 

Substituting for I from Eqn, (A,1 .8),(A.1.9), (A. 1, 10), (A,1 .28) and 
(a. 1.32) we get the following integral in n 

—2 

J X” ^ (A-n)^ e”’^ dn (A.4.2) 

A 

where v ® ^ ‘f’ B - Bj sin ^ cos B) 



(a.4.3) 
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V = V(G ) = {£ (3G: - o'*) . i_ 

p ? 24 p pj 4g 


C5G + 


®p same as given by Eqn, (-^,3,11), 

Oai-iylDg out the 5 ihtegration and rearranging terme given the 
Integral (a, 4.8) as 


C^.4.9) 


,5/2 - 

J V 


u. . 
^3 


{VCa,X2) - VCa,X )} 


where 


(a.4.10) 


Vp = V(a,Xp} = [(XjCa.Xp) -Jx^ X2(?,Xj,)} 


* fXj Co.X,] - X4(d.X 1)] 


- 3g (S {i?2(o,xp - ♦iCd.X^)} * 2 {♦iCa.y 




r 


-Xp -fjCd.Xp)}] - Ij- [3{Ziff,X^) - X Z„Ca,X )} 


+ 12(2j(a,Xp3-Xp ZjCa.Xp)} ♦ 4 (25(0, X^) + Z^(o,Xp)}]} 


^^'‘ihaliy we can write down the collision integral for this case as 


2 -2nj n. 


I.. ( ^(25)" „ JliJlr 

7 i ^ -f A ^ 


m 


, 3/2 . 8s ,1/2 , 8s 1/2 


" 2““" ,;4 “ TkCTt.+Tt.}^ "kt- ^ ^kT 


‘'1 -j 




/ gS /2 u (v^.vp exp {. - 5 — > dr, 


(a. 4. 12) 
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Integrating (a. 4.2) we get 

V lj2(A, V (1+2^) - ~ (1+2/?) erf A - Ae'^ } (A.4.4) 

Integration with arespect to a ; 

This is done by using new variable A, as indicated earlier. After the 
A integration we get 



where 


/ ^ '(. 2(^|dA = ^ {V(V2) - V(v,)} 
1/ « -1^ 


V 

P 


V(v ) = (3v 

p 12 p 


- S' ''P 


(A.4.5) 


i (3Vp + 2v^) erf v^} (A.4.6) 

Integration in & space ; 

Here again the procedure used in the earlier case is adopted to obtain 
the required integral as 


f5/2 “ CO G+Uj • 

J t .tn [ ^3 


“ /^IGX - V,) G e"“^ dv 

/T TT ^ A 


(a.4.7) 


“ij 0 


Integration with respect v gives us the above integral as 


1^/2 -00 

/ (Vj - V. + V3 - V4)G df 

u. . 0 

13 


(A.4.8) 


where 
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^IMs iSj s fold integral, to be evalmted nmerically. 

Ftjr the case v/hen i = vfe can obtain the expression for the collision 
integral as 


I. . = 

xx 2 


-/-n? 


^2 I j2 ^ j 

t ^4]c T. ^ ^ 


kX 




V ex^> {-a - 


_ sC2B)^ ^ sC2I}^)2 _ 


K dv dT, 


?k T, 


{a.4.'I5) 


■IMS is a seven fold integral, to be evaluated numerically. 


Sins we have all the collision integral appearing in the 
moment e<5uations given in Chapter 4.. 
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''p ' ''"’p’ " 4 " 


2G^D e P -^(3+12(^ + Ach er£ G„} 
p 48 P P P 

and G are same as given by Ban. (",3.11 ). 

P 

Carrying out the G integration and rearranging terms gives the 
integral (A, 4.8) as 




{VCa,X2) - V(a,Xp} 


where 


vir 


Vj, » VCcXj,) . [(XjCcXp) -3Xp X2(»,Xp)) 

+ (Xj - Xp X4(?.Xp)>] 

- Ib ('!'2t“>Xp) - Xp ♦jCS'.Xp)) . 2 {*,ff.Xp) 

-Xp ♦j&Xp)}] - -ij [S{Zi(S-.Xp) - Xp Z„C?,Xp)) 

* lZ{23(a,Xp]-Xp Z2C“.Xp)) ♦ * (Xsff.Xp) * Xp Z^Ca.Xp))] 


Finally we can write down the collision integral for this case as 


.2 •-2n4 n- 

1 = 


13 2 


s(2B)^v ^ ““1 "j m 


/ 2kCTt.^. 


\5/2 f 8s . 1/2 . 8s 1/2 t 

)V V. ^ 


/ j5/2 7 {VV,} =XP t- ■‘7 


2k T„ 2k ' 6 

Vy 


(^.4.9) 

(a. 4.10) 


2 

ij 

(A.4-12) 



MPMDIX B 


Some of the importaat integrals used in Appendix A have been listed 
here. 


3.1 Integrals = /x^ e”^ dx 


- erf X 



A - (erf X - X e'^ )/2 

2 

= - I (1 + x^) 

2 , 

h, = tc 3 erf X - e~^ (3x + 2x^) > 

4 

f ^ 2 -x^ 

= - f {% + 23i ■¥ 2) e 

2 

t “*5 erf X - (I5x + 10x^+ 4x^)e“^ > 

o o 

B.2 Integrals #^(x) = / x^ erf x dx : 

jlX 


2 

{§)^ = X erf X + e /2 

2 

=: i {(2x^-1) erf X + xe**^ ) 

@2 = - { (2x^ erf X + (x^+1 ) e ^ > 

-x^ 

(H) ^ ( 4 x^- 3 ) erf X + (2x +3x)e ^ ^ 

3 16 

- ^ 

@ = 2x^ erf X + (x^+ 2x^+ 2 ) e } 

4 10 

* j_{( 8 x^_ 15 ) erf X + (4x^+ 1 Qx^+ 15 x) e 


5 
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This is, a six fold integral, to he evaluated numerically. 

Ebr the case when i = j, v/e can obtain Ihe expression for the collision 
integral as 


>>i'i2 -4nj 

I ,. (^=^1 - — ( 


11 


sC2Dj''i _ “"“3^ ^ E_)3/2 V V(V2) 


2 ^4 Mk 


kT, 


_ _ s(2B)^ + s(2Dj)^ _ _ _ 

V exj) {-a K dv dt 

2k T„ 


(A.4.13) 


This is a seven fold integral, to be evaluated numerically. 


Thus we have all the collision integral appearing in the 
moment equations given in Chapter 4 .. 
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B*3 Integra I & 

X„C«,x) = /" g" dG 

7= /irTa 

X^= (1+2 - x^) x^/2a 

X^ = (3 + 2 Cl ) X X^/2 Ct 

= {3 + 12a 4aSc’^) Xq/4«^ 

^ — (15 + 20ci X -f- 4^^x"^)x Xq/4 

Xg - {15 + 90ot x^+ 60 + Q x /ea^ 

0 

B.4 Integrals T)^Cct,x} = /” G^ g-«(G-x) -G^ 

“ J^/Cl+a) exp {-ciX^/(l+a)| 

% = ax 

^2 ~ Cl+a+2a^x^} i^f^/2(a+l)^ 

- {3a (l+a) + 2a^ x^} x if/^/2(l+a)^ 

^4 = {3(a+l)^ + 12 Ca+i}a^ i^^3/4 (a+l)"^ 


^5 = {15o:(I+a)^ + ZOa"* (l+a)x^ + 4a^ x^} x iJIq/ 4 (a+l)"^ 



APPHJDIX B 


Some of the important integrals used in Appendix A have been listed 
here# 


B.1 



B.2 



_ 2 

Integrals A^(x) = e”^ dx : 

= erf X 

= -rrx2/2 2 

_ (erf X - X e“^ 0/2 

1 2\ -X 

= - i (1 + X ) e 

2 , 

== i{ 3 erf X - e"^ (3x + 2x^) > 

2 

1/4 ,-.2 \ —x 

=: -;^(x4-2x + 2 ;e 

2 

* ^ {15 erf X - (I5x + 10x^+ 4x^)e~^ > 

Integrals <B)j^(x) = / x^ erf x dx j 
= X erf X + e /2 


^ {(2x^-1 ) erf X + xe”^ } 

2 

®o “ T {(2x^ erf x + (x^+1 ) e~^ }• 

2 

(g)^ = { |g- ( 4 x^- 3 ) erf X + (2x^+3x)e‘* ) 

2 

®4 “ erf X + (x^+ 2x^+ 2) e"^ } 

®5 * ^ + ( 4x^+1 0x^+1 5x) e''^ 
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-V 


B,5 Integrals 2(a^x) = /“ erf r 


dG 


^ putting y = )?■ fouomng reaiate are listed. 


S “ 


% 


eri y 


■3. = /a /n'(ci +T" ! + X erf y) > 




erty+ 


2^6 


^ (2a X + 3)x 
2a 


2 t'o(1 + a. 
erf y + 


> Xo 


{gx^g^ ^ (5 + 2x^)a+ 2} ^ 

4 (a 4- 

r / ^ 4 2 1 

L V4 a"x + 12a x +5) erf y/4 


] Xo 




+ f 8x^4' 4- (l2x\ 12x) ot^+ (ax^+ 26 x) a^ 

2 

4* (2x^4- 12x) -a + 5x } e~^ A + 1 >] 

I { + 20a x^+ 15> X erf y/4a^ 

4- f (l5a^+ 20a 8} (1 + a)^ + 2x^ (a + 1 ) a 

{JG^ 4- 5Ca^ -f 35a + 9) + 4x^a^ ^0<i^+ 10a^4 5 


erf y + 3e ^ ) Xq /(“ + 1 ) 


5 .7/2 


+■! )}e ^ ]X(j« 


wiiers 
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B,5 Integrals 

XnC“,x) * / dG 

*«co 

f= /tt/oc 

X1 = ^Xo 

p 

Xg = (1 + 2CI X ) XQ/2a 

X^ « (3 + 2ctx^) X Xy2a 

X^ « {3 + 12a x^+ 4a^x'^) X^Aa^ 

X^ = (15 + 20a x^+ 4a^x‘^)x X^A 

Xg = (15 + 90a x^+ 60 a^x'^ + 8 oPx^) X^/8a^ 

2 p2 

B.4 Integrals if)yjC«»x) - f 

»«> 

2 

= /x/ (l+a) ejcp {-ax /(1+a)} 

= ax ij;y (l+a) 

1^2 (l+a+2a^x^) i|jQ/2(a+l)^ 

413 = {3a (l+a) + 2a^ x^} x t|>^/2(l+a)^ 

11^4 = {3(a+l)^ + 12Ca+l)a^ + 40^^ x"^} \l>^/4 {a+l)^ 

- {15a(l+a)^ + 20a^ Cl+a)x^ + 4a^ x^) x i(»Q/4{a+l)^ 
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A « 


15 .3 


4- 


4*^ ^■ 


f o^ (^ + l) X^ + C(^ (os + -)) 


2 4 

X 


3 '* ^\3 ^ 45 2 45 2 2 ^ 

^ a 1^1 4- 3} rt (a + 1 ) 


15 „£ ^4 


a'*^ (cs + 1 )^ 4- ~ 3 + ax^ (a + 


Ji 


45 2 


X’f15a('x + 1) 3ci(a+l) x} 


STTTT) - -C 4r 4--^ (a+ 1 ) 


1*^ 4 -O ^ 

^ X a (1 -*-a) }y+{15sa+10x 

a(a -f 1 )} (1 + y^)/a (a + l") - { -^ cc + ^ ^ 

4-2 

«(a +1)! (3y X 2y'^) + ( 3x (y + 2y + 2) 

1 ^ a (a ' 4 - t ) - ^ ( 1 ^ 10y^+ Ay^) y * 





B.5 Integrals Z(a,x) = /“ g” etf G 


dG 


By putting y = ^ -j+cT ^ following results are listed. 


Z = 


erf y 


Xo 


= {e~^ /2 )^ar{a~rT~) * x erf y) Xq 

„ , (l+2)ix^) (2a + l) X -y^ , 

Z » { .i L erf y + — e '>' } 

2a 2/a(1 + (xf 

_ r (2a x^+ 3)x „ 

Z„ = [-: erf y + 

5 2a 


{6x^a^ + (346x^)a^ + (5 + 2x^)a+ 2> „-y^ t 

e j Xo 

4(a + 1 a^'^^ 

( (4 a^x^ + 1 2a x^ + 3 ) erf y/4a^ 


Z. = 


+ { ax^+ a"^ + (l2x^+ 12x) a^+ (ax\ 26x) 

+ (2x\ 19x) a + 5x } e""^ /{ 4a^/^(a + 1 

t { 4a^x'^ + 20a x^+ 15) x erf y/4a^ 

+ { (15«^+ 20a + 8) (l + a)^ + 2 x^ (a + 1 ) a 

(30a^ + 50a^ + 35a + 9) + 4x^a^ (5o^+ 10a\ 10a^+ 

2 

(a erf y + Be~^ ) /(® + 1 


5a +1 )} 


v/here 
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APPMDK C 


C.t Golllsion Integrals D. . and li i 
- - 0^.1 


!Rroia Chapter 5 v/e hare 


B. . 
13 


a-t 


/ f f y(- g,a) K d n dx. dx 

i J *4^ 4^ X 


(c,l *1 ) 


Poll owing the procedure given in Appendix A, 16 integrations were 
carried, out for the case i ^ to get the six fold integral 


D. - 
IJ 








{R CS,X.)-R3(a,X^)>exp { 


5 (25) 

2k T,, 


4^'v■ 2xk{T. +T..) 

^JL ^3 

2 s C22j^) ^ 

} sin 0 sin d6 dpj^ dB dEj^ d^ 


2k 7 


V 


3 


(G*1.2) 


where 

R^Ca^X) = ^ J" tX4 - J' X3) ’ ^ f-6x Z3, 

+ 6Z2 - 4x ^3 + 4 Z4) - ^ C-2x 

, 2ii ^ - 2i>2 + 2^3) (0,1.?) 

ffla .iuaatltles and J appearing in Ikjn. (0.2.2) have 

already been defined in Appendix A and the integrals x„ , 

listed in Appendix B- 
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+ — a'' (a + 1 ) (“ + “I ) * 

+ a^d + a)3 {a + ^ ) 

+ J|. a^ x^ (a + 1 f + ^ « + ccx^ (a + 0 

15 

8 * 

B = { -^ a^x + 1 5 a^(a + 1 ) x^+ 5 a^(a + 1 )^ } 

4 

/a(a + 1 ) - t ^ a2 4 ^ x^ a2 (a+ 1 ) 

+ x^ a^(l + a)^} y + {15x a+ lOx^ 

a(a + 1)} (1 + y^)v^ (a 7T) - { ot + x^ 

'3» 4 2 

a(a + 1 )) (3y + 2y^) + { 3x (y + 2y + 2) 

/o.(a +T) - -g (15 + 10y^+ 4y^) y . 





5br i 5= 3 ¥/e have 


.i = }V2 I j2 


4TTk T* 


^sd i- [ 


s(2:i)^ + s(2B,)2 


'ffheT‘B 


sin 6 sin dg dB^ dB dE^ do 


(C.1.4) 


»sd C“.*) ' r Cxj - X Xj) - X P, . 2x Zj*Zj.xZ„)*l (tj.xtj) (c.,.5) 


Here again all the quantities used are defined in either Appendix A 
or B* Tile collision integral in the atoiaic species conservation 
equation laay be obtained as a special case of B. . by effecting following 

i J 

substitution as indicated In Chapter 2 or 3* 


B, = 0, = 0, to = 0, B = 0 


(G.1,6) 


C*2 Collision Integrals ) sind ^ ^ * 

'ie can write the expression for integrals H. s{ 2B) ) straight 
sway using results of section C*1 as follov/s. 


®5r i 3^ 3 


1 '? m 

H-- ds(2E)^) - { SL 

^ 4 *V 2iTk(Tt +T^ ) 

i i 


j3/2 J_ 


/ CR 3 Ca,X,) - RJ5,X,) 4sC2rO^} exp U 

^ ^ ^ ^ 2k Tv 


s{2h r 

• " • ) sin B sin Bj dB dBj dB dBj d'^i 


CC.2,1) 
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APPETOIX a 


C.1 Collision Integrals P. , and D j 

■ 


Prom Chapter 5 we have 




S / f.f . Y (- g.n) K d n dt, dr 

X 3 m 


(C.1 .1 ) 


Following the procedure given in Appendix A, 16 integrations were 
carried out for the case i j, to get the six fold integral 


Hi nj 

Dij - - — < 


m 


>3/2 j^/2 {R (S,xj-Rs(«,X2)}exp { 

A/S' 2TTk(T^.+T^.) iiij 


'i "'3 


sC2B) 


2 s(2B^r 


2k T. 


'v. 2k Tv. 
1 ''j 


-} sin 6 sin f3j^ dB dB dBj^ d(!) 


( 0 . 1 . 2 ) 


whei’e 


/n 1 

RsCa,X 3 = { j- (X4 - X X3) - C-6x 


- - 2i}j2 + 2\|)^) 


24 *0 


(c.1 . 3 ) 


The quantities a, X.j and J appearing in Eqn. (C.2.2) have 
already been defined in Appendix A and the integrals j 
listed in Appendix B. 
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:^r 


(|sC2n^= - J ? 


A 2 s( 2 B) + s(2H,)^ 


^sd s(2S)^} eXT {[- 


sill B SiTi dB dB dB^ d<t) 


(g, 2.2) 


Idr the integrals H (2 the method indicated in Appendix A 

ij] » 

is folXcHved and .ve can write for the final six fold integral for the 
case 1 3 as 


«ij =/ (jJi <) dt 


1 T ,.,2^ 


^ _ _i 1 ^ m ^ 3/2 1 j j3/2 

4/7 27rk(T^ +T^ ) “13 ^ 

i j 

— _ sC 2B)2 

{Pj, (a,Xj) - (a.X^;)} exp {- — 


s(2E,3" 

} Sin $ sin dp dp^ dB dE, d(^ d(J)^ 


(c.2.3) 


^ (X4 “ 5C X3) - [^3 - (3-x) Z2 


CC*2.4) 


Cx + i) Zi + f ZJ M I Ml - f) - X ^(,1) 


Bar i = 3 we have 
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Bor i = 0 we have 


D-i = ~nf E_)V2 j j2 

4Trk T.. 

sC2r5)^ + s(2B,)2 

P-sd ®3cp {- [ ]} 

2k T,. 

sin e sin dB dBj^ dD dBj^ d(Jij^ 


(C.1.4) 


where 


^sd (X 3 “ X X 2 ) “ (^3 ” 2x 2^2^^!'' (c.1*5) 


Here again all the quantities used are defined in either Appendix A 
or B, The collision integral in the atomic species conservation 

equation may be obtained as a special case of by effecting following 
substitution as indicated in Chapter 2 or 5* 


= 0, (*>j = 0, « = 0, = 0 


(C,1 .6) 


Integrals ^ ' 


2 

We can write the expression for integrals s( 2B) ) straight- 

away using results of section C.1 as follows. 

Ibr i j 

tA . - 

l 3/2 1 


1 2 

R^. (|s(2B)^) = 1 { 


tn 


/ (RsC5,X,) R_(a,X 2 ) {^s(2B)^} exp { 


4 VIT a-iikCTt +T^ ) 
i 3 
:i 




u 


13 


2-' ^2 


s(2Br 

2k T„ 


( 0 . 2 , 1 ) 


; :,s(2B ) ■ : . 

- sin B sin Bj^ ds dBj^ dB dBj^ d<^ d^j 


V, 

3 






t 


kJ" 
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a - n? n, 

^ ^ ‘Uk T 


/ 0^ 


1 


s(2r)"- + sC2tp' 

r,x) ^.xp (- — 1 sin g sin 6, 

2k T„ ' - 




(C.2.5) 


tk-a integrals from S, .(§) the same substitiition 

Oii: 13 

is aqn* (C.1.6) is usei. 

fes i^ve all the reauire4 collision integrals for the 


luring case* 
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and for i = j 


R ds(2B)^)= - n? | j2 

i 4TrkTti 

„ 5 o s(2B)^ + s(2Bv)2 

R-. (a,X) {ls(2B)^} exp {[ 


2k T, 


■]> 


sin 3 sin d 6 d 3 j^ dD dBj d(J) di})j (c, 2 . 2 ) 

2 

Ibr the integrals H. . (i Wj_) the method indicated in Appendix A 

IL Q i 

is followed and we can write for the final six fold integral for the 
case i ^ j as 




where 


{- 


m 


-} 


"i 


4/jr 2Trk (T^ +T^ ) 
i j 


3/2 j j3/2 ■ 


u 


iJ 


_ sf2D)2 

{Rj. (a,X^) . R (a,X )} exp {- 

s(2B,)2 

— } sin 6 sin Bj d3 dBj^ d3 dBj^ d({) d<J)j^ 


R,.(a,X) = I CX4 - X X3) - [Z3 - (3-x3 Z2 


- (X + j) Zj + f ZJ + [ I + (1 - f) h 


(G.2,3) 


(C,2,*4) 


Star 1 = 3 we have 





FLAST !<: bPH[ KOC\ L INL't'A 



ijP^ACE. 


RPPPF PCN^AliCN 
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Rll 4 ».f) = - nj vTT I 

*• 4Trl; 

9 ^ 

s(2r0 + s( 2 Ej )2 

> P sin 3, 

2k Ty ' ^ 

i 

dg dBj dB dBj^ d({) dijtj. 


(C.2.5) 


!Do obtain the integrals ^Qj, 2 (^) from the same substitution 

suggested in Eqn, (c.1,6) is used. 

®b.us we have all the required, collision integrals for the 
dissociating case. 




104 



FIG, 2.2.1 ^.SURFACE REPRESENTATION,' 
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FIG.5,1.1_EQUIUBRIUIv| DEGREE OF DISSOCIATION FOR; 
NITROGEN BEHIND A NORMAL SHOCK WAVE 
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fig' 7,4:2. -ROTATIONAL TEMPERATURE THICKNESS FOR 

Ov ■■V'V/'C :,- : nitrogen . ' ' . , 
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